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FORMATION  OF  POROUS  SILICON  CARBIDE  AND  ITS  SUITABILITY  AS  A 
CHEMICAL  AND  TEMPERATURE  DETECTOR 


Tilghman  Lee  Rittenhouse,  Ph.  D. 

Department  of  Chemistry 
University  of  Illinois  at  Urbana-Champaign,  2004 
Paul  W.  Bohn,  Advisor 

The  need  to  sense  chemical  mixtures  in  a  variety  of  hostile  environments  (such  as  high 
temperature  caustic  gases)  continues  to  grow.  However,  silicon  electrical  devices  are  limited  to  relatively 
low  temperatures  (<  250  °C).  For  this  reason,  wide  bandgap  materials  such  as  silicon  carbide  have 
received  increased  attention.  Current  SiC  sensors  such  as  Schottky  diodes  composed  of  catalytic  metals 
show  deficiencies  such  as  unacceptable  drift  in  the  signal.  Alternative  sensor  structures,  such  as  porous 
semiconductors,  may  provide  improved  sensor  performance. 

A  novel  electroless  method  of  producing  porous  silicon  carbide  (PSiC)  is  presented.  Unlike 
anodic  methods  of  producing  PSiC,  the  electroless  process  does  not  require  electrical  contact  during 
etching.  Rather,  platinum  metal  deposited  on  the  wafer  before  etching  serves  as  a  catalyst  for  the 
reduction  of  a  chemical  oxidant,  which  combined  with  UV  illumination  injects  holes  into  the  valence 
band,  the  holes  subsequently  participating  in  the  oxidation  and  dissolution  of  the  substrate.  The  etchant  is 
composed  of  HF  and  K2S2O8  in  water.  Various  porous  morphologies  are  presented  as  a  function  of 
etchant  concentration,  time  of  etching,  and  SiC  polytype.  Wafer  quality  is  of  the  utmost  concern  when 
utilizing  the  electroless  wet  etchant,  since  defects  such  as  stacking  faults,  dislocations,  and  micropipes 
have  a  large  impact  on  the  resulting  porous  structure. 

Results  of  imaging  and  spectroscopic  characterization  are  presented  and  compared  to  PSiC 
produced  via  anodic  etching  of  the  same  wafer  material.  In  general,  it  is  found  that  electrolessly  etched 
PSiC  has  photoluminescent,  cathodoluminescent,  and  Raman  scattering  properties  relatively  unchanged 
from  bulk  SiC.  However,  the  spectroscopic  properties  of  the  anodically  etched  PSiC  were  dramatically 
different  than  those  found  for  bulk  SiC.  Cathodoluminescence  studies  of  anodically  etched  porous  layers 
yield  luminescence  peaks  in  the  ultraviolet  region  (energies  higher  than  the  bandgap  energy  of  bulk  SiC). 
Additional  spectroscopic  results  (photoluminescence,  Raman  emission)  have  been  collected  from  the 


same  porous  films  in  order  to  determine  the  origin  of  the  UV  cathodoluminescence  emission.  The 
photoluminescence  from  these  samples  showed  a  wide  range  of  peak  luminescence  and  can  be  linked  to 
the  degree  of  etching  in  the  porous  film.  Raman  scattering  spectra  show  the  development  of  Frdhlich 
modes  as  etching  proceeds.  Therefore,  although  quantum  confinement  is  one  possible  cause  of  the  UV 
cathodoluminescence,  other  results  indicate  that  defect  structures  or  surface  state  species  are  the  most 
likely  origin  of  the  UV  emission. 

Anodically  etched  PSiC  films  are  then  investigated  for  their  suitability  as  chemical  and 
temperature  detectors  particularly  suited  for  high  temperature,  caustic  environments.  Studies  utilizing 
interferometry  as  well  as  conductivity  have  been  conducted.  PSiC  layers  were  unresponsive  in 
interferometry  for  a  wide  variety  of  adsorbate  gases,  and  this  method  of  analyte  detection  proved 
unsuccessful. 

Ohmic  contacts  to  PSiC  were  investigated  in  order  to  study  changes  in  the  bulk  conductivity  of 
the  semiconducting  skeleton  when  exposed  to  various  chemical  and  temperature  conditions. 
Conductometric  investigations  show  that  PSiC  has  a  measurable  response  to  high  concentrations  of 
hydrogen  gas  (~6%  decrease  in  resistance  when  exposed  to  20%  H2)  comparable  to  currently  used 
platinum  alloy  resistive  films.  However,  similar  to  the  resistive  films,  strict  temperature  control  of  the 
PSiC  sensor  is  necessary.  Current  hydrogen  sensing  arrays  may  be  improved  by  the  introduction  of  a 
PSiC  sensor  for  the  specific  role  of  high  concentration  sensing.  Eventually,  it  may  be  possible  to  form  an 
array  of  PSiC  sensors  capable  of  distinguishing  the  concentrations  of  individual  components  within  a  gas 
flow  as  well  as  the  temperature  of  the  gas  flow. 
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CHAPTER  1 


INTRODUCTION 


1.1  Motivation 

Today,  the  need  for  high  temperature  sensing  of  hydrogen  within  gas  mixtures  is  finding 
increased  attention.  Within  the  aerospace  industry  specifically,  there  is  a  need  for  sensors  that  can  be 
used  in  fuel  leak  detection,  fire  safety  monitoring,  and  engine  emission  monitoring.  For  example,  the 
ability  to  monitor  the  exhaust  of  a  jet  or  rocket  engine  would  allow  for  much  greater  safety  and  engine 
performance.  However,  these  exhausts  can  be  expected  to  be  at  high  temperature  and  contain  caustic 
species  in  addition  to  the  hydrogen  to  be  detected.  A  solid  state  device  that  can  monitor  these  exhaust 
gases  in  situ  does  not  exist. 

The  current  widespread  use  of  silicon  in  electronics  (primarily  due  to  the  facility  with  which  an 
insulating  oxide  can  be  formed)  combined  with  the  ability  to  obtain  highly  pure  silicon  inexpensively  has 
made  it  the  most  studied  semiconductor  to  date.  However,  silicon  carbide  has  several  intrinsic  properties 
that  make  it  favorable  to  silicon  for  some  applications,  especially  those  that  are  at  elevated  temperature  or 
in  caustic  environments,1'3  as  shown  in  Table  1.1.  Si  is  typically  limited  to  approximately  250  °C  in 
electronic  applications  and  about  600  °C  in  mechanical  devices;  Si  microelectromechanical  systems 
operating  in  excess  of  250  °C  would  require  the  use  of  bulky  cooling  schemes.4  However,  due  to  its 
stability  at  high  temperature  and  in  harsh  chemical  environments,  SiC  can  operate  in  conditions  in  which 
Si  fails.  Sensing  applications  that  play  to  silicon  carbide’s  strengths  include  rocket  engines,  nuclear 
reactors,  materials  growth  chambers,  automotive  systems,  and  other  high  temperature,  caustic 
environments.  For  these  reasons,  SiC  was  chosen  as  the  sensor  material  for  this  work. 

SiC-based  distributed  control  electronics  are  envisioned  to  eliminate  90%  of  the  wiring  and 
connectors  needed  in  conventional  sheltered-electronic  aircraft  control  systems.5  For  the  U.S.  Air  Force, 
SiC  sensors  implemented  on  an  F-16  would  result  in  800  pounds  of  reduced  weight.  In  addition,  the 
occurrence  of  wiring  and  connector  problems  (the  most  frequent  cause  of  propulsion  maintenance  action 
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Table  1.1.  Properties  of  Si  and  SiC 


Melting  Point  (Si)/  Thermal  Breakdown  Saturation  Electron 

Sublimation  Point  (SiCd  Conductivity  at  Electric  Field  Drift  Velocity  Bandgap 


Material 

298  K  fW/m*Kl 

(V/m) 

at  298  K  (m/s) 

(eV) 

Si 

1410 

150 

3x  107 

0.9  x  10s 

1.1 

4H/6H-SiC 

2700 

370-490 

2-4  x  108 

2.0  x  10s 

3-3.3 

and  downtime  in  aircraft  today)  would  be  greatly  reduced.  For  commercial  jetliners,  the  increased  fuel 
efficiency,  weight  savings,  and  reduced  pollution  would  carry  large  economic  and  environment  payoffs. 

It  has  been  estimated  that  the  economic  savings  value  of  high  temperature  SiC-based  electronics  would 
result  in  savings  of  millions  of  dollars  per  aircraft  over  the  course  of  an  aircraft’s  multi-decade  operational 
lifetime. 

1.2  Current  Crystalline  Semiconductor  Sensor  Studies 

The  literature  details  numerous  efforts  utilizing  crystalline  wafers  of  silicon  carbide  to  detect 
hydrogen  gas  flow  at  various  elevated  temperatures.  Most  of  these  studies  have  revolved  around 
conductometric  methods  of  detection.  SiC  sensors  for  hydrogen  gas  have  taken  the  form  of  either 
Schottky  diodes,  metal-insulator-SiC  capacitors,  or  occasionally  field  effect  transistor  (FET)  geometries. 

The  Schottky  diodes  utilize  a  metal  gate  material  (usually  the  catalytic  metals  Pt  or  Pd)  bonded 
directly  to  the  SiC  substrate  with  only  a  very  thin  (~2-5  nm  thick)  oxide  separating  the  metal  from  the 
semiconductor.  Hydrogen  or  hydrogen-containing  species  dissociatively  adsorb  to  the  metal  and 
hydrogen  atoms  migrate  to  the  metal-SiC  interface,  creating  a  polarized  layer  on  the  semiconductor 
surface.  The  work  function  of  the  metal  decreases,  accompanied  by  a  lowering  of  the  Schottky  barrier 
height  at  the  metal/SiC  interface.  Furthermore,  in  a  Schottky  diode  structure,  it  is  possible  the  thinness  of 
the  oxide  allows  for  tunneling  of  electrons  through  the  junction.  The  result  is  an  exponential  change  in 
the  current  or  a  quadratic  change  in  the  capacitance  of  the  device.  The  change  in  conductivity  is 
measured  and  can  be  correlated  to  surface  concentrations  and  to  the  levels  of  the  sampled  gas  in  the 
ambient.  The  advantage  of  a  Schottky  diode  sensing  structure  in  gas  sensing  applications  is  its  high 
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sensitivity.  However,  numerous  studies  have  shown  that  interfacial  reactions  between  the  gate  metal  and 
the  SiC  change  the  diode  properties,  and  the  metal-SiC  structure  must  some  how  be  stabilized  in  order  to 
eliminate  drift  in  sensor  properties  for  long  term  operation  at  elevated  temperatures. 

Metal-insulator-SiC  (MISiC)  devices  utilize  a  thicker  insulating  layer  between  the  metal  gate  and 
the  SiC  substrate.  Often  the  insulating  layer  is  thermally-  or  ozone-grown  oxide  Si02  (MOSiC),  but  other 
insulating  layers  such  as  AIN6’7  have  been  studied.  Recent  studies  have  also  investigated  the  use  of 
reactive  insulators  (such  as  Sn02)  as  the  barrier  between  the  gate  metal  and  the  SiC  substrate.8 
Regardless  of  the  insulating  material  used,  the  adsorbing  gas  changes  the  space  charge  region  under  the 
metal  clusters  at  the  metal-insulator  interface.  This  thicker  insulating  junction  does  not  allow  for 
tunneling,  but  rather  relies  solely  upon  a  field  effect  for  its  sensing  capabilities.  The  presence  of  the 
insulating  layer  aids  in  reducing  reactions  between  the  gate  metal  and  the  underlying  semiconductor 
substrate.  However,  interfacial  reactions  have  still  been  shown  to  have  detrimental  effects  on  device 
performance.  Also  as  a  result  of  the  relatively  thick  insulating  layer,  MISiC  devices  have  lower 
sensitivity  than  Schottky  diode  devices. 

Field  effect  transistor  devices  again  rely  upon  a  fairly  thick  insulating  layer  at  the  gate,  covered 
by  a  catalyzing  metal  such  as  Pt.  However,  now  instead  of  a  capacitor  configuration,  additional 
electrodes  are  configured  as  a  source  and  a  drain.  A  voltage  applied  between  the  source  and  drain 
contacts  causes  a  current  to  flow  between  the  contacts,  through  a  channel  region.  The  concentration  of 
mobile  carriers  in  the  channel  region  is  responsive  to  the  ambient  gas  environment  surrounding  the  gate 
electrode.  Unfortunately,  these  field  effect  transistor  devices  suffer  from  many  of  the  same  problems 
encountered  with  MISiC  devices,  namely  interdiffusion  reactions  and  reduced  sensitivity  as  compared  to 
Schottky  diodes. 

Therefore,  despite  rapid  improvement  over  the  past  decade,  each  of  these  devices  has  encountered 
technological  hurdles  that  still  need  to  be  overcome  before  an  operable  device  can  be  manufactured  and 
installed.  These  shortcomings  will  be  discussed  and  illustrated  further  in  Chapter  4  of  this  dissertation. 
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1.3  Current  Porous  Semiconductor  Sensor  Studies 


Although  first  discovered  in  1956,9  porous  semiconductors  have  received  particularly  increased 
attention  in  the  past  decade.  Porous  semiconductors  are  attractive  candidates  for  several  potential 
applications;  their  inherent  high  surface  area  makes  porous  semiconductors  particularly  attractive  for 
sensing  applications.  However,  few  reports  detail  the  investigation  of  porous  silicon  (PSi)  sensors,  and 
even  fewer  studies  of  porous  silicon  carbide  (PSiC)  sensing  have  been  reported.  It  appears  that  the 
widespread  use  of  a  porous  semiconductor  is  most  likely  years  away. 

One  vein  of  research  with  PSi  sensors  has  utilized  Fabry-Perot  interference  fringes.10"13  When 
light  is  reflected  from  a  porous  film,  the  reflected  spectrum  will  show  fringes  whose  spacing  is  dependent 
upon  the  refractive  index  of  the  porous  film  (Fig.  1.1). 11  This  detection  method  has  been  successful  for 
measuring  biomolecules  as  well  as  vapors  of  organic  solvents  such  as  ethanol  in  the  hundreds  of  ppb.  It 
has  been  proposed  that  capillary  condensation  is  the  primary  cause  of  refractive  index  changes  in  the  PSi 
film.12  However,  this  method  has  not  been  shown  to  be  effective  in  measuring  gases  that  have  refractive 
indices  very  near  that  of  air  or  inert  gases  (such  as  H2).  In  addition,  signal  from  these  devices  can  drift  as 
silicon  in  the  porous  film  oxidizes  (Fig.  1.2).12'14 

Similar  to  the  Fabry-Perot  interference  fringe  work  described  above,  Snow  et.  al.  has  made  PSi 
devices  that  act  as  Bragg  mirrors.15  Snow  et.  al.  theorize  that  capillary  condensation  of  organic  vapors 
(saturated  vapors  of  acetone  and  chlorobenzene)  lead  to  the  relatively  large  (75-125  nm)  shifts  in 
reflection  fringes  they  witness. 

PSi  devices  have  also  been  studied  via  conductometric  studies.  In  general,  conductometric 
devices  can  be  expected  to  be  more  straightforward  and  easier  to  apply  in  many  sensing  applications  than 
optical  interference  measurements.  Schechter  et.  al.  conducted  some  of  the  earliest  conductometric 
studies.16  Schechter  et.  al.  report  large  changes  in  conductivity  through  a  Schottky  contact  when  a  PSi 
sensor  is  exposed  to  a  vapor  with  a  large  dipole  moment  (such  as  methanol  or  acetone).  After 
investigating  a  variety  of  adsorbate  gases,  Schechter  et.  al.  concluded  that  the  increase  in  conductivity  was 
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Figure  1 . 1 . 1 1  Reflection  spectra  of  two  porous  silicon  thin  films  with  different  thicknesses,  in  air  (solid 
curves)  and  after  exposure  to  1.7%  ethanol  in  air  (dashed  curves).  The  Fabry-Perot  interference  fringes 
shift  as  a  result  of  the  change  in  the  dielectric  constant  of  the  porous  layer. 
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Figure  1.2. 14  Reflection  spectra  from  a  porous  silicon  film  during  the  growth  of  native  oxide  over  the 
course  of  ten  days.  The  Fabry-Perot  fringe  shifts  to  the  blue  upon  oxidation  of  the  PSi  film,  representing 
a  change  in  baseline  for  the  device  over  time. 
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dependent  upon  the  adsorbate’s  dipole  moment.  They  further  proposed  that  the  increases  in  conductivity 
are  due  to  an  increase  in  the  density  of  charge  carriers  within  the  porous  nanostructure. 

The  results17  of  Kim  et.  al.  are  similar  in  nature  to  those  reported  by  Schechter.  Kim  et.  al.  report 
ethanol,  methanol,  and  acetone  vapor  detection  near  room  temperature  with  porous  silicon  devices. 

Kim’s  FET  devices  are  rectifying  in  nature,  at  least  partially  due  to  interfaces  between  n+  silicon  with  p- 
type  porous  silicon  regions.  With  5  V  applied,  Kim  et.  al.  report  a  seven-fold  increase  in  current  flow 
when  the  device  is  exposed  to  0.5%  ethanol  (Fig.  1.3). 17  They  propose  that  capillary  condensation  of  the 
analyte  gases  within  the  porous  layer,  accompanied  by  a  reduction  of  the  built-in  potential  of  the  p-n 
junction,  is  responsible  for  the  change  in  conductivity.  Their  initial  studies  also  show  a  strong 
dependence  of  the  device’s  response  to  processing  of  the  PSi  surface.  It  was  found  that  oxidized 
(hydrophilic)  porous  silicon  layers  are  relatively  insensitive  to  organic  vapors  as  compared  to  reduced 
(hydrophobic)  porous  silicon  layers. 

Angelucci  et.  al.  have  researched  porous  silicon  membranes  that  are  permeated  with 
semiconducting  oxides  such  as  SnC>2.  Their  devices  have  shown  80-100%  increases  in  current  when 
exposed  to  30  ppm  of  benzene  and  CO.18  Furthermore,  the  porous  devices  have  shown  a  sensitivity  to 
humidity. 

A  thin  (20  nm),  continuous  layer  of  gold  was  used  as  the  top  electrical  contact  to  porous  silicon  in 
a  device  fabricated  by  Foucaran  et.  al.19  The  Schottky  contacts  showed  a  response  to  acetone  and  ethanol. 
Their  results  indicate  that  the  thin  gold  layer  on  the  porous  silicon  surface  does  not  affect  the  adsoiption 
of  polar  molecules.  Like  Angelucci’s  device,  Foucaran’s  experiments  indicate  an  irreversible  response  of 
the  sensor  to  humidity. 

The  investigation  of  porous  SiC  devices  has  been  more  limited.  Shields  et.  al.20  investigated  a 
PSiC  sensor  in  the  presence  of  hydrocarbons  at  200-400  °C.  This  work  was  able  to  detect  0.5%  methane 
and  propane  and  rudimentarily  discriminate  between  the  two  gases.  Their  contacts  were  made  of 
chromium,  and  the  I/V  characteristics  indicate  fairly  Schottky  behavior  of  the  contacts  to  the  PSiC. 
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Figure  1.3. 17  I/V  curves  measured  for  ethanol  in  air  exposed  to  a  bare  porous  silicon  film. 

The  sensing  studies  referenced  above  all  indicate  that  PSiC  sensors  show  promise  as  a  suitable 
device  for  engine  exhaust  monitoring.  If  the  high  surface  area  of  a  PSiC  film  can  be  properly  harnessed, 
porous  semiconductor  devices  have  the  ability  to  provide  high  sensitivity  at  high  temperatures  and  caustic 
environments  where  current  sensors  cannot  operate.  However,  many  unresolved  issues  remain.  For 
example,  all  of  the  conductometric  studies  reported  above  utilized  Schottky  contacts  to  the  porous 
network.  Therefore,  it  is  difficult  to  distinguish  whether  the  changes  in  conductivity  were  due  to 
interactions  between  the  adsorbate  gas  and  the  porous  skeleton,  or  due  to  reduction  of  the  Schottky  barrier 
height  at  the  metal-semiconductor  interface. 

1.4  Outline  of  Work 

Initially,  this  work  investigated  electroless  etching  of  SiC  to  form  PSiC.  After  discovering  a 
suitable  electroless  etching  technique  for  PSiC,  this  research  progressed  to  measure  the  spectroscopic 
properties  (cathodoluminescence,  photoluminescence,  and  Raman  scattering)  of  the  porous  material. 
Particularly,  cathodoluminescence  and  photoluminescence  were  conducted  to  look  for  light  emission  from 
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the  PSiC  that  was  blue  of  the  bandgap  emission.  Such  emission  could  be  indicative  of  nanocrystallite 
formation.  Raman  studies  were  conducted  to  further  investigate  quantum  confinement  effects  in  the 
nanoparticles. 

After  characterizing  the  PSiC  formed  with  the  electroless  etchant,  this  work  then  proceeded  to  a 
study  of  the  formation  of  PSiC  by  anodic  etching  under  conditions  similar  to  those  reported  in  the 
literature.  These  samples  served  to  enable  a  comparison  of  the  spectroscopic  properties  between 
electrolessly  etched  and  anodically  etched  PSiC,  because  the  anodically  etched  pieces  came  from  the 
same  wafer  as  the  PSiC  etched  with  the  electroless  etchant.  One  highlight  of  this  section  of  the  research 
was  the  discovery  of  ultraviolet  emission  in  cathodoluminescence.  Anodic  etching  was  also  found  to 
produce  porous  layers  of  a  more  reproducible  and  controllable  thickness. 

This  work  next  proceeded  to  investigate  the  sensing  capabilities  of  PSiC.  The  anodic  etching 
process  was  fine  tuned  until  PSiC  films  could  be  produced  that  yielded  Fabry-Perot  interference  fringes  in 
their  reflectance  spectra.  Shifting  of  Fabry-Perot  fringes  due  to  changes  in  the  dielectric  constant  of  the 
film  was  sought  to  provide  sensing  capabilities  complementary  to  conductometric  measurements.  The 
porous  films  were  investigated  for  their  optical  sensing  properties  with  several  gaseous  species  at  a  range 
of  elevated  temperatures. 

Finally,  conductivity  measurements  were  made  on  PSiC  layers  in  a  variety  of  gas  atmospheres  at 
elevated  temperatures.  One  major  hurdle  during  this  phase  of  the  work  was  the  formation  of  ohmic 
contacts  to  the  PSiC.  Results  from  the  conductivity  test  series  were  evaluated  for  their  ability  to  produce 
an  array  of  PSiC  sensors  able  to  identify  and  quantify  a  combination  of  gases  within  a  range  of  elevated 
temperatures. 
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CHAPTER  2 


ETCHING  AND  MORPHOLOGY 


2.1  Introduction 

Due  to  their  unique  morphological  and  optical  properties,  porous  semiconductors  are  candidates 
for  several  varied  applications.  For  example,  porous  semiconductors  have  been  investigated  as  compliant 
substrates  for  epitaxial  growth.1'6  Some  research  indicates  that  epitaxial  layers  grown  on  top  of  porous 
substrates  may  lead  to  fewer  defects  as  well  as  strain  relaxation.2  In  addition,  because  of  their  inherent 
high  surface-to-volume  ratios,  porous  semiconductors  are  attractive  candidates  for  chemical711  or 
physical  (temperature,  pressure,  flow)12  sensing.  The  suitability  of  porous  semiconductors  as  epitaxial 
substrates  or  sensors  is  highly  dependent  upon  the  porous  nanostructure. 

The  literature  documents  several  early  variations  of  forming  porous  semiconductors  that 
eventually  lead  to  the  aqueous  anodic  methods  commonly  used  today.  PSi  has  been  formed  with 
anhydrous  etchants  such  as  HF/acetonitrile  mixtures. 13-14  Also,  stain  etching  of  Si  to  form  thin  PSiC 
layers  utilizes  A1  films  in  the  presence  of  an  aqueous  HF/HN03  etching  solution.15  Electroless  etching  of 
GaN  has  been  demonstrated  using  K2S20g/K0H  etching  solutions.16  Finally,  as  a  forbearer  to  current 
PSiC  production,  Shor  et.  al.  conducted  a  series  of  experiments  investigating  the  use  of 
photoelectrochemical  etching  of  SiC  in  the  hope  of  finding  faster,  smoother,  and  more  selective  bulk 
etching  processes.17"20 

Semiconductors,  such  as  silicon  and  gallium  nitride,  can  easily  be  made  porous  via  wet  etching 
techniques  (both  electrochemical  and  more  recently  metal-assisted  electroless  etching  techniques21"24), 
and  a  variety  of  porous  morphologies  can  be  formed.  In  silicon,  the  morphology  is  highly  dependent 
upon  the  etching  conditions.  In  gallium  nitride,  the  etching  conditions  have  significant  impact  on  the 
resulting  morphology,  but  the  porous  structure  is  also  highly  dependent  upon  the  particular  properties  of 
the  wafer  (doping,  defect  density,  etc.),  since  the  material  quality  of  available  gallium  nitride  does  not 
approach  that  of  commercially  available  single  crystal  silicon.  Morphologies  of  porous  materials  range 
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from  sponge-like  networks  of  pores  to  very  regular  arrays  of  long  columnar  nanocrystals/pores  oriented 
perpendicular  to  the  semiconductor  free  surface.  Regardless  of  the  morphology,  though,  all  porous 
semiconductors  are  characterized  by  a  high  surface  area. 

Currently  available  anodic  etching  techniques  produce  PSiC  of  varying  morphology  and 
properties.25  45  The  results  reported  in  the  literature  detail  specific  changes  in  morphological  and 
spectroscopic  properties  upon  formation  of  a  porous  nanostructure  in  SiC;  however,  wide  discrepancies 
often  exist  in  the  properties  measured  by  independent  researchers.  As  a  result,  definitive  conclusions 
about  the  causes  of  spectroscopic  phenomena  in  PSiC  are  difficult  to  determine.  In  addition,  the  need  for 
electrical  contacts  in  the  anodization  process  is  a  limitation  on  the  large-scale  production  of  PSiC.  The 
development  of  an  electroless  wet  etching  chemistry  could  permit  large  areas  of  homogeneous  PSiC  to  be 
produced,  and  the  electroless  technique  is  better  suited  for  small  or  odd  shaped  pieces.  Electroless 
etching  also  offers  the  advantage  of  masking  areas  to  remain  unetched.  Finally,  electroless  etching  has 
shown  the  ability  to  produce  PSi  with  tunable  optical  properties.24 

The  electroless  etching  of  Si  and  GaN  is  made  possible  by  using  ultrathin  (d  <  1 0  nm)  films  of  Pt 
metal  to  catalyze  the  reduction  of  a  chemical  oxidant  in  the  presence  of  HF,  which  combined  with  UV 
illumination  injects  holes  into  the  valence  band,  the  holes  subsequently  participating  in  the  oxidation  and 
dissolution  of  the  substrate.  FECF  is  used  as  the  oxidant  for  Si  and  GaN.  However,  although  H2O2 
produces  PSiC  anodically,  preliminary  research  showed  that  SiC  etches  very  inefficiently  in  HF:H202 
etchant  without  anodization  (even  in  the  presence  of  Pt  catalyst).  The  more  active  oxidant  potassium 
persulfate,  K2S208,  is  required  for  the  efficient  electroless  production  of  PSiC.  K2S208  in  the  presence  of 
Pt  and  illumination  is  competent  to  inject  holes  into  the  valence  band  of  SiC.  The  holes  may  then  travel 
laterally  and  participate  in  etching  reactions  in  a  manner  similar  to  that  observed  in  PGaN.16  One 
hypothesis  is  that  SiC  is  oxidized  in  an  8e"  process  analogous  to  the  electroless  etching  of  PSi  and  PGaN, 
the  products  subsequently  reacting  with  HF  to  produce  a  soluble  product  according  to: 

Cathode  (Pt): 

4S2082'  +  hv  ->  8S042'+  8h+ 
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Anode  (SiC): 


SiC  +  4H20  +  8h+  ->  Si02  +  8H+  +  C02 
Si02  +  6HF  ->  H2SiF6  +  2H20 

Overall: 

4S2082'  +  SiC  +  2H20  +  6HF  ->  8S042'  +  8H+  +  C02  +  H2SiF6 
It  is  also  possible  that  the  ultraviolet  irradiation  utilized  during  etching  is  creating  a  short-lived  diffusible 
oxidant  species  (a  free  radical  form  of  SOx,  for  example)  that  participates  in  the  etching  process;  no 
mechanistic  studies  have  been  investigated  in  this  work  to  verify  the  exact  reaction  mechanism. 

Silicon  carbide  forms  over  250  different  polytypes.  Polytype  refers  to  the  stacking  sequences 
within  the  crystalline  material.  The  polytypes  of  silicon  carbide  generally  fall  into  three  primary 
categories:  cubic,  hexagonal,  and  rhombohedral.  The  cubic  form  of  SiC,  3C-  or  (3-SiC,  is  face-centered 
cubic  with  alternating  layers  of  carbon  and  silicon  atoms  in  the  [111]  direction.  However,  this  polytype  is 
difficult  to  crystallize  and  was  not  used  in  this  research.  Rather,  two  forms  of  hexagonal  SiC  were  used 
exclusively.  Most  work  was  conducted  with  6H-SiC,  which  has  alternating  layers  of  silicon  and  carbon 
in  the  sequence  ABCACBABCACB. . .  along  the  c-axis  of  the  hexagonal  cell.  Limited  work  was  also 
conducted  with  the  4H-SiC  polytype,  which  has  alternating  layers  of  silicon  and  carbon  in  the  sequence 
ABAC  ABAC. . .  along  the  c-axis  of  the  hexagonal  cell. 

The  etching  reaction  proposed  above  does  not  explicitly  indicate  that  Pt  must  be  deposited  onto 
the  SiC  surface  in  order  for  electroless  etching  to  proceed.  However,  Pt  catalyzes  the  reduction  of  the 
oxidant  (in  our  case,  S2082")  and  acts  as  a  local  cathode  site  for  hole  injection  into  the  semiconductor  (Fig. 
2.1).  In  addition  to  the  Pt  catalyst,  UV  illumination  is  also  necessary  during  electroless  etching.  The 
energy  of  the  holes  available  from  S2082  reduction  is  not  sufficient  to  accomplish  direct  injection  of  holes 
into  the  valence  band  of  SiC.46,47  Fig.  2.2  shows  the  energy  of  the  S20827S042'  couple  on  the  absolute 
scale  referenced  to  the  work  function  of  Pt.46'47  Clearly  holes  generated  through  this  oxidation  process  do 
not  possess  sufficient  energy  to  access  the  SiC  valence  band  directly. 
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2.2  Experimental  Methods 

2.2.1  Etching  Techniques 

2.2.1. 1  Electroless  Etching  Technique 

Silicon  carbide  wafers  were  obtained  from  (Tree  and  included  a  6H  n-type  wafer  (p  =  0.086  Q 
cm),  a  4H  n-type  wafer  (p  =  0.067  Q  cm),  and  a  4H  p-type  wafer  (4.4  x  1017  cm'3).  Both  n-doped  wafers 
were  cut  on  axis  to  expose  the  (0001)  Si  plane;  the  p-type  wafer  was  cut  8°  off  the  c-axis.  The  two  n- 
doped  wafers  were  polished  on  the  silicon  side  only;  the  p-type  wafer  was  double  polished.  Initial 
cleaning  of  the  wafers  consisted  of  consecutive  5-min  dips  in  water,  acetone,  and  methanol  before  sputter 
deposition  of  the  Pt  catalyst.  Pt  was  sputter  coated  onto  the  SiC  wafers  using  a  Desk  II  TSC  Turbo 
Sputter  Coater  (Denton  Vacuum)  through  a  physical  mask  to  produce  small  (500  pm  in  diameter,  1.5  mm 
center- to-center)  circles  of  10  nm  mass  thickness. 

Initially,  wet  etching  was  attempted  with  CH30H:HF:H202  (1:2:1  by  volume  of  pure  methanol, 
49%  HF,  and  30%  H202),  since  this  etchant  is  known  to  etch  both  Si  and  GaN  electrolessly.21'24 
However,  this  etchant  produced  only  a  small  degree  of  pitting  in  SiC  (Fig  2.3).  Previous  work  conducted 
in  our  laboratory23  has  shown  that  careful  precleaning  of  the  GaN  surface  is  required  in  order  to  produce 
porous  morphology  consistently.  Therefore,  a  more  aggressive  precleaning  (nitric  acid)  of  the  SiC 
samples  was  conducted;  however,  the  same  sparse  pitting  resulted  from  exposure  to  the  CH30H:HF:H202 
etchant.  Next,  an  etchant  composed  of  0.1M  K2S208:0.08M  KOH  was  attempted,  since  this  etchant  had 
been  shown  to  electrolessly  etch  GaN.16  This  etchant  produced  SiC  with  a  rougher  surface  (Fig.  2.4); 
however,  the  morphology  was  not  porous  and  did  not  extend  beyond  the  top  surface  of  the  SiC  wafer. 

The  working  hypothesis  at  this  point  in  time  was  that  the  etching  reaction  needed  to  be  initiated, 
since  it  was  known  that  the  CH30H:HF:H202  etchant  can  produce  PSiC  porous  when  SiC  is  used  as  a 
substrate  below  an  etched  GaN  layer  (Fig.  2.5).48  Therefore,  the  surface  of  the  SiC  wafer  was  roughened 
with  a  scribe  prior  to  immersion  in  CH30H:HF:H202  in  order  to  produce  initiation  sites  for  etching. 
However,  SiC  prepared  in  this  manner  did  not  exhibit  porosity  either. 
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Figure  2.1.  Schematic  diagram  of  the  reduction  of  persulfate  ion  accompanied  by  hole  injection  into  the 
SiC  wafer  at  the  Pt  catalyst  deposited  on  the  surface. 


S20827S042 


Figure  2.2.46'47  Scale  showing  the  reduction  potential  of  the  persulfate  ion  in  relation  to  the  conduction 
and  valence  bands  of  crystalline  SiC  on  the  absolute  electrochemical  energy  scale. 
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Figure  2.3.  Plan  view  SEM  image  of  the  relatively  sparse  etch  pits  resulting  from  electroless  etching  of 
SiC  in  CH30H:HF:H202  etchant. 


Figure  2.4.  Plan  view  SEM  image  of  the  roughened  SiC  wafer  surface  after  exposure  to  the  0.1M 
K2S208:0.08M  KOH  etchant. 
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Figure  2. 5. 48  Plan  view  SEM  image  of  PSiC  formed  when  an  epilayer  of  GaN  grown  on  a  SiC  substrate 
was  entirely  etched  away  with  the  CIROIEHFiEECC  etchant. 

In  these  initial  experiments,  even  though  a  great  majority  of  the  surface  remained  planar  and 
unetched,  extensive  porosity  in  the  immediate  vicinity  (within  500  nm)  of  micropipes  was  often  observed 
(Figs.  2.6  and  2.7),  indicating  a  possible  relationship  between  crystalline  quality  and  porous  morphology. 
One  hypothesis  is  that  the  porous  structure  in  the  vicinity  of  micropipes  is  due  to  dislocations  terminating 
at  the  free  surface,  since  these  micropipes  obviously  interfere  with  the  short-range  order  of  the  lattice 
structure. 

Some  researchers  have  noted  the  formation  of  a  “skin  layer”  on  the  surface  of  porous  SiC  during 
anodic  etching.25,32,34’37’49'50  The  top  thin  skin  layer  remains  primarily  crystalline,  with  relatively  few  pits 
extending  to  the  free  surface  of  the  wafer.  This  skin  layer  is  approximately  20-50  nm  thick,  and  it  covers 
a  fully  porous  morphology  that  extends  deeper  within  the  wafer.  Therefore,  several  of  the  electrolessly 
etched  SiC  samples  were  reactive  ion  etched  (RIE)  to  remove  any  potential  skin  layer.  RIE  was 
conducted  at  a  power  density  of  0.8  W/cm2  and  with  20  seem  flow  of  SF6  through  the  reactor;  the  total 
pressure  in  the  reactor  was  20  mtorr  (maintained  by  argon).  Etching  times  ranged  from  2  to  6  min.  These 
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Figure  2.6.  SEM  image  of  a  micropipe  on  the  SiC  wafer  surface  after  exposure  to  the  CFEOFFFlFiFECE 
etchant.  Notice  the  porous  ring  extending  approximately  1  pm  laterally  from  the  periphery  of  the 
micropipe. 


Figure  2.7.  Close-up  SEM  image  of  the  porous  microstructure  around  the  periphery  of  a  micropipe  that 
has  been  exposed  to  the  CFEOFFElFiFECE  etchant. 
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conditions  were  selected  to  remove  approximately  1 00  nm  of  material  from  the  surface  of  the  SiC 
wafer.37  51  However,  RIE  revealed  no  definable  submerged  porosity,  although  the  RIF,  results  were  hard 
to  interpret  definitively,  because  the  RIE  damage  to  crystalline  material  (Fig.  2.8)  produced  a  morphology 
very  similar  to  that  in  porous  material.  All  of  the  samples  exposed  to  RIE  looked  like  the  normal  RIE 
damage  shown  in  Fig.  2.8,  with  the  exception  of  one  sample  (Fig.  2.9).  However,  it  is  believed  that 
anomalies  in  the  reactive  ion  etching  produced  the  structure  seen  in  Fig.  2.9,  and  that  no  submerged 
porosity  existed  in  this  wafer  either. 

Oxidation  has  also  been  cited  in  the  literature  as  a  method  of  exposing  submerged  porosity.25,26 
Several  of  the  samples  in  this  work  were  therefore  oxidized  (both  before  and  after  RIE).  Oxidation  was 
carried  out  with  steam  (in  02)  flowing  at  5  L/min  within  a  tube  furnace  maintained  at  1 150  °C;  the 
oxidation  time  was  typically  1  hr.  However,  there  was  no  evidence  of  any  submerged  porosity  after 
oxidation  and  a  dip  in  HF.  Finally,  several  techniques  of  cycling  RIE,  oxidation,  and  electroless  etching 
were  conducted  in  an  attempt  to  initiate  the  etching  process.  None  of  these  techniques  were  successful  in 
producing  extensive  porosity,  however. 

Eventually,  an  etchant  composed  of  HF  and  K2S2O8  was  tested.  This  etchant  was  investigated 
simply  to  combine  the  effects  of  the  acid  attack  of  SiCF  by  HF  with  K2S2O8,  which  is  a  slightly  stronger 
oxidizer  than  H2O2.  This  etchant  proved  the  most  successful  and  was  used  for  the  remainder  of  the 
electroless  etching  work.  S2O82"  ( E°  =  -2.05  V)  has  a  very  similar  standard  reduction  potential  to  H202  (E° 
=  -1.776  V),  and  one  would  not  expect  them  to  behave  much  differently  as  etchants  on  thermodynamic 
grounds.  We  are  still  not  certain  what  gives  rise  to  the  significantly  stronger  etch  behavior  of  S2082" 
relative  to  H2O2. 

Once  S2O82"  was  identified  as  a  suitable  oxidant,  two  different  concentrations  of  wet  etchants 
were  tested:  (a)  concentrated  HF  etchant,  consisting  of  14M  HF  (semiconductor  grade,  Transene  Inc.)  and 
0.075M  K2S208  (Fisher  Scientific);  and  (b)  dilute  HF  etchant,  consisting  of  0.50M  HF  and  0.15M  K2S208 
(the  solubility  limit  of  K2S2O8  at  room  temperature).  Tests  were  done  with  methanol  in  the  etchant  to 
lower  surface  tension  and  prevent  the  formation  of  bubbles  on  the  SiC  surface  during  etching.  However, 
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Figure  2.8.  Plan  view  SEM  image  of  a  previously-unetched  crystalline  SiC  wafer  that  has  been  etched  via 
reactive  ion  etching  using  the  conditions  described  in  the  text. 


Figure  2.9.  Plan  view  SEM  image  of  SiC  that  has  been  exposed  to  the  CE^OEEHFiEECE  etchant  for  60 
min,  followed  by  reactive  ion  etching  using  the  conditions  described  in  the  text.  The  wafer  had  shown 
only  sparse  etch  pits  on  the  surface  after  electroless  etching. 
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methanol  did  not  help  the  etching  process;  in  fact,  SiC  did  not  etch  with  substantial  methanol  added  to  the 
etchant. 

Etch  times  varied  from  1 0  to  60  min.  Ultraviolet  irradiation  was  obtained  from  a  Hg  arc-lamp 
(Oriel  Instruments  Model  6281).  The  arc-lamp  emission  supplied  for  the  duration  of  the  etching  at  a 
distance  of  10  cm  was  short-pass  filtered  to  expose  the  wafer  to  32  mW/cm2  at  X  <  360  nm.  At  the 
conclusion  of  etching,  a  rinse  of  5  min  each  in  water  and  methanol  (2x)  was  conducted  before  air  drying 
the  porous  samples. 

2.2. 1.2  Anodic  Etching  Technique 

The  same  wafers  and  pre-cleaning  procedures  were  used  for  anodic  etching  of  SiC,  with  the 
exception  that  the  Pt  pattern,  deposited  onto  the  SiC  wafers  prior  to  electroless  etching,  was  eliminated. 
Rather,  to  prepare  the  SiC  pieces  to  undergo  anodic  etching,  a  corner  of  the  wafer  (~2  x  2  mm)  was 
coated  with  20  nm  mass  thickness  of  Pt  for  electrical  contact  with  the  etching  probe. 

The  SiC  was  placed  into  a  Teflon  electrochemical  cell,  along  with  a  Pt  wire  counter  electrode  and 
a  saturated  calomel  reference  electrode  (Fig.  2.10).  Dilute  HF  (5%)  was  employed  as  the  etching 
solution.  The  SiC  sample  to  be  etched  was  illuminated  by  UV  light  (15  mW/cm2  at  X  <  365  nm)  for  the 
duration  of  the  etching  (usually  40  min).  Upon  the  conclusion  of  etching,  a  rinse  of  5  min  each  in  water 
and  methanol  (2x)  was  conducted  before  air  drying  the  porous  samples. 

2.2.2  Scanning  Electron  Microscopy 

Scanning  electron  microscope  (SEM)  images  were  collected  using  a  high-resolution  field 
emission  SEM  (Hitachi  S4700)  using  accelerating  voltages  of  10-15  keV. 

2.2.3  BET  Isotherms 

Surface  area  measurements  were  conducted  in  order  to  further  characterize  the  porous  structure, 
as  well  as  evaluate  the  suitability  of  the  porous  layers  to  serve  as  chemical  sensors.  One  method  of 
measuring  the  surface  area  of  porous  materials  uses  the  method  developed  by  Brunauer,  Emmet,  and 
Teller.  BET  isotherms  are  collected  by  adsorbing  nitrogen  to  the  surface  of  the  sample  material.  Surface 
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area  measurements  were  obtained  at  liquid  nitrogen  temperature  using  an  ASAP  2400  apparatus 
(Micromeritics).  Each  sample  was  outgassed  for  24  hours  at  150  °C  prior  to  analysis. 
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Figure  2.10.  Diagram  of  the  etching  apparatus  used  to  etch  SiC  wafers  anodically. 
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2.3  Results  and  Discussion 


2.3.1  Electrolessly  Etched  SiC 

2.3.1. 1  n-doped  6H-SiC 

SEM  analysis  indicates  that  etch  times  ranging  from  10  to  60  min  yield  a  progression  of  PSiC 
morphologies.  In  general,  for  both  the  concentrated  and  dilute  HF  etchants,  etch  times  less  than 
approximately  30  min  yield  a  ridged  porous  structure,  as  shown  in  Fig.  2.11.  The  porous  morphology  is 
homogeneous  across  the  entire  wafer  surface,  although  occasional  pockets  of  varying  degrees  of  porosity 
are  observed.  The  cause  of  the  pockets  of  relatively  unetched  regions  is  unknown  at  this  time.  These 
isolated  patches  of  lightly  etched  SiC  do  not  occur  consistently  at  a  specific  distance  away  from  the  Pt 
circles.  The  relatively  unetched  regions  could  result  from  regions  of  low  dislocation  density  at  the  surface 
of  the  wafer.  Or,  it  is  possible  that  the  wet  etchant  simply  did  not  have  good  access  to  these  small  regions 
of  the  wafer,  possibly  due  to  the  formation  of  bubbles  during  the  etching  process. 

Often,  wafers  etched  for  these  relatively  short  amounts  of  time  were  actually  more  porous  directly 
between  the  Pt  dots,  rather  than  at  the  Pt  dot  edge.  This  is  an  unexpected  result  given  the  theory  of  hole 
injection  occurring  at  the  Pt  cathode.  It  seems  unlikely  that  hole  injection  occurs  at  the  Pt,  and  then  the 
holes  are  transported  macroscopic  (mm)  distances  away  to  produce  etching.  Flowever,  Pt  is  necessary  for 
etching,  as  several  experiments  were  conducted  without  Pt,  and  no  PSiC  was  ever  observed  in  the  absence 
of  Pt. 

Longer  etch  times  (~60  min)  yield  different  morphologies  depending  upon  the  etchant  used. 

When  using  the  concentrated  HF  etchant,  porosity  restricted  to  the  near-surface  region  is  obtained  (Fig. 
2.12).  However,  when  using  the  dilute  HF  etchant  at  longer  etch  times  a  spongelike  morphology  results 
(Fig.  2.13).  One  key  difference  between  the  morphologies  obtained  at  short  and  long  times  appears  in  the 
lateral  etching.  At  shorter  etch  times  the  etching  occurs  inhomogeneously  in  the  plane,  but  is  largely 
directed  perpendicular  to  the  original  surface.  At  longer  times  the  etching  develops  more  lateral  structure 
with  etching  proceeding  to  undercut  some  of  the  structures  confined  to  the  surface. 
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Figure  2.1 1.  Plan  view  SEM  image  of  PSiC  prepared  from  n-type  6Fl-SiC  by  electroless  etching  in 
concentrated  F1F  etchant  for  20  min. 


Figure  2.12.  Plan  view  SEM  image  of  PSiC  prepared  from  n-type  6Fl-SiC  by  electroless  etching  in 
concentrated  F1F  etchant  for  60  min. 
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Figure  2.13.  Plan  view  SEM  image  of  PSiC  prepared  from  n-type  6Fl-SiC  by  electroless  etching  in  dilute 
F1F  etchant  for  60  min. 

Numerous  porous  wafers  were  cleaved  and  examined  in  cross-section  to  determine  the  internal 
depth  of  the  porous  structure.  Despite  the  depth  of  the  porous  layer  reasonably  clean  cleaves  could  be 
obtained.  As  shown  in  Fig.  2.14,  a  uniformly  porous  layer  can  be  formed  deep  into  the  bulk  of  the  wafer, 
extending  >100  pm  in  especially  long  etches.  Furthermore,  the  porous  structure  extends  all  the  way  to 
the  free  surface  of  the  wafer,  without  the  formation  of  a  skin  layer  reported  in  anodic  etching  experiments 
(Fig  2. 15).  All  of  the  porous  materials  had  the  pattern  of  networked  spherical  or  triangular 
pores/nanocrystallites  as  shown  in  Fig.  2.15;  cylindrical  nanocrystallites/pores  were  never  formed.  For 
the  deeply  porous  layers,  average  nanocrystallite  sizes  on  the  order  of  25-50  nm  diameter  were  obtained. 

The  porous  SiC  to  crystalline  SiC  interface  took  on  several  forms  when  examined  in  cross- 
section;  occasionally  multiple  forms  of  porosity  termination  were  witnessed  in  just  one  piece  of  porous 
SiC.  First,  relatively  infrequently  porosity  stops  abruptly  (at  least  to  the  viewer  of  the  exposed  fracture 
plane)  (Fig.  2.16).  Second,  the  porosity  can  end  at  a  front  interspersed  with  relatively  large  pits  (Fig. 
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2. 17).  Finally,  the  porous  structure  can  just  gradually  become  smoother  until  the  broken  cross  section 
becomes  planar  (Fig.  2.18). 

The  formation  of  deep  porous  structure  is  striking,  because  reactants  (HF)  need  to  be  transported 
to  the  evolving  etch  sites,  and  products  need  to  be  transported  away.  The  postulated  formation  of  CCT  as 
an  oxidation  product  undoubtedly  impacts  the  transport  process.  More  importantly  the  interconnected 
nature  of  the  deep  porous  structures  provides  multiples  paths  for  the  transport  of  reactive  species  to  the 
active  etch  front  and  for  transport  of  products  away.  In  PGaN  where  the  pores  are  largely  columnar,  deep 
structures  are  difficult  to  obtain.  Thus,  it  is  quite  likely  that  the  lateral  connections  between  pores 
facilitate  transport  and  make  the  attainment  of  deep  porosity  possible. 

Another  conspicuous  feature  of  the  electroless  etching  of  PSiC  is  the  lateral  extent  of  etching. 
PSiC  is  formed  adjacent  to,  and  as  much  as  millimeters  away  from,  the  edges  of  the  Pt  coated  areas.  Fig. 
2. 19  illustrates  the  difference  in  etching  behavior  between  the  region  covered  by  Pt  and  the  areas  adjacent 
to  it.  In  this  sample,  after  formation  of  PSiC,  the  Pt  was  removed  with  aqua  regia,  and  a  sharp  interface 
(±50  nm)  is  revealed.  In  the  region  adjacent  to  the  Pt  porosity  exists  both  at  the  surface  and  below  the 
surface.  In  the  region  covered  by  Pt  during  etching,  no  surface  porosity  is  observed,  although  the 
modulated  light  and  dark  areas  in  the  SEM  micrograph  are  a  clear  sign  that  some  of  the  subsurface  pores 
extend  laterally  under  the  edge  of  the  Pt  coated  region.  These  results  indicate  that  the  Pt,  in  addition  to 
serving  as  a  catalyst,  can  be  used  as  a  mask  to  direct  the  formation  of  surface  porosity  patterns.  Other 
experiments  were  conducted  in  which  a  continuous  10  nm  thick  sheet  of  Pt  was  deposited  onto  the  back 
(unpolished)  side  of  the  SiC  wafer  prior  to  etching.  This  protocol  resulted  in  porosity  on  the  front 
(polished)  side  of  the  wafer,  but  the  porosity  was  restricted  to  a  shallow  surface  layer. 

Unfortunately,  the  etch  rate  (depth  of  etch)  was  inconsistent  from  sample  to  sample.  As  shown  in 
Fig.  2.14,  deep  porous  structures  could  be  formed  in  some  samples.  Flowever,  in  many  samples, 
etching/porosity  was  limited  to  the  top  50  nm  at  the  surface,  regardless  of  how  long  the  sample  was 
exposed  to  the  etchant  (Fig.  2.20).  For  this  reason,  using  currently  available  commercial  SiC  wafers, 
electroless  etching  does  not  result  in  the  reproducible  morphology  necessary  to  mass  produce  PSiC  films. 
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Figure  2.14.  Cross-sectional  SEM  image  of PSiC  prepared  from  n-type  61 1-SiC  by  electroless  etching  in 
dilute  F1F  etchant  for  60  min.  Note  the  uniformity  of  the  etched  morphology  throughout  the  depth  of  the 
etched  region. 


Figure  2.15.  Cross-sectional  SEM  image  of  PSiC  prepared  from  n-type  6Fl-SiC  by  electroless  etching  in 
dilute  F1F  etchant  for  60  min. 
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Figure  2.16.  Cross-sectional  SEM  image  of  the  abrupt  termination  of  porosity  at  the  PSiC/SiC  interface. 


MB45Bfly2  15.0kV  10.4mm  xIOOk  SE(U)  5/10/02  13:35 _ 500nm 


Figure  2.17.  Cross-sectional  SEM  image  of  the  PSiC/SiC  interface  that  is  marked  by  the  presence  of 
large  etch  pits. 
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Figure  2.18.  Cross-sectional  SEM  image  of  a  gradual  transition  between  PSiC  and  unetched  SiC;  no 
definable  interface  of  porosity  is  discemable. 


Figure  2.19.  Plan  view  SEM  image  of  PSiC  prepared  from  n-type  6Fl-SiC  by  electroless  etching  in  dilute 
F1F  etchant  for  60  min.  The  image  is  obtained  near  the  edge  of  a  Pt-covered  region,  but  after  removal  of 
the  Pt. 
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Figure  2.20.  Cross-sectional  SEM  image  of  PSiC  prepared  from  n-type  6Fl-SiC  by  electroless  etching  in 
dilute  F1F  etchant  for  60  min.  Despite  the  long  etch  time,  only  a  surface  layer  of  pores  are  observed. 

Fligh  resolution  X-ray  diffraction  (F1XRD)  studies  were  conducted  in  order  to  determine  if  some 
feature  of  crystallographic  quality  (particularly  defects)  could  be  responsible  for  the  inconsistencies  in 
etch  depth.  F1XRD  studies  can  reveal  changes  in  lattice  constant,  defect  density,  and  impurities  (>  1%) 
for  the  typical  X-ray  sampling  depth  of  ~25  pm.  The  results  of  the  FIXRD  study  of  the  porous  samples 
were  inconclusive.  The  unetched  6Fl-SiC  wafer  was  found  to  have  higher  defect  densities  at  the  edges  of 
the  wafer,  which  is  a  normal  result  of  crystal  growth.  Upon  electroless  etching,  plan  view  SEM  images 
indicated  that  samples  from  the  edges  of  the  SiC  wafer  did  become  more  porous  than  samples  from  the 
central  regions  of  the  SiC  wafer.  Flowever,  upon  cleaving  the  samples,  cross-section  analysis  showed  that 
none  of  the  pieces  examined  in  this  test  series  exhibited  submerged  porosity. 

The  FIXRD  test  series  also  showed  a  change  in  the  lattice  constant  of  the  SiC  after  etching.  The 
porous,  etched  material  exhibited  a  variation  of  c-direction  lattice  spacing  between  the  value  for  the 
unetched  material  (15.08  A)  and  a  0.06%  contraction  of  the  lattice  parameter  (Fig.  2.21).  The  change  in 
lattice  parameter  is  similar  to  the  behavior  of  etched  PSi,  where  lattice  parameter  also  changes  upon 


etching.  Interestingly,  though,  in  PSi  the  lattice  parameter  expands  upon  the  formation  of  pores,52"56  while 
our  results  indicate  a  small  lattice  contraction  in  PSiC.  This  contraction  of  the  lattice  parameter  leads  to 
tensile  stress  within  the  porous  layer.  This  residual  stress  can  affect  the  Raman  spectra  of  these  porous 
layers,  and  so  this  effect  must  be  accounted  for  when  analyzing  the  Raman  data  presented  in  Chapter  3. 

2.3. 1.2  n-doped  4H-SiC 

The  4H  polytype  of  SiC  yields  a  distinctly  different  morphology  than  6H  material  when  exposed 
to  the  same  electroless  etching  conditions.  The  resulting  morphologies  are  independent  of  the 
concentration  of  etchants  used.  For  short  (~20  min)  etch  times  a  ridged  porous  structure  similar  to  the  6H 
polytype  is  obtained  in  4H  wafers  (similar  to  that  shown  in  Fig.  2. 1 1).  For  longer  etch  times  (~60  min), 
different  pieces  of  the  4H  wafer  yield  different  morphologies,  from  a  cellular  surface  exhibiting 
nanofmgers  (Fig.  2.22)  to  highly  porous  spongelike  structures  (Fig.  2.23).  One  possible  explanation  for 
the  variation  in  morphology  observed  when  etching  4H  material  under  nominally  identical  conditions 
arises  from  the  observation  that  the  4H  wafer  studied  contained  at  least  some  regions  of  6H  material  as 
determined  by  Raman  mapping  experiments  (see  Chapter  3).  It  is  reasonable  to  hypothesize  that 
variations  in  etched  morphology  between  different  4H-SiC  samples  arise  from  variations  in  the  polytype 
distributions  of  the  nominal  4H  starting  material.  The  principal  difference  between  electroless  etching  of 
4H  and  6H  wafers,  however,  is  the  consistent  confinement  of  the  porous  structure  to  the  near-surface 
region  in  4H  material.  At  most  a  surface  layer  of  pores  a  few  hundred  nm  deep  is  obtained  in  4H  material 
(Fig.  2.24). 

Occasionally  samples  that  contained  both  the  4H  and  6H  polytypes  at  the  surface  were  etched.  It 
was  found  that  at  the  completion  of  etching,  only  the  4H  regions  had  etched  and  become  porous;  the  6H 
regions  remained  planar  and  unetched  (Fig.  2.25).  This  is  an  unexpected  result,  given  the  ability  of  pure 
6H-SiC  to  etch  much  more  completely  than  the  apparently  “pure”  4H-SiC  pieces.  In  many  regions  the 
4H/6H  boundary  was  marked  with  relatively  large  micropipes  (Fig.  2.26).  In  addition,  in  cross-section, 
the  4H  polytype  was  detected  to  undercut  the  6H  polytype  within  the  wafer  (Fig.  2.27). 
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Figure  2.21.  Fligh  resolution  X-ray  diffraction  rocking  curve  of  unetched  6FI-SiC  (dashed  trace),  a  PSiC 
sample  with  only  a  surface  layer  of  pores  (dotted  trace),  and  a  PSiC  sample  that  had  a  porous  layer  greater 
than  100  pm  thick  (solid  trace). 


Figure  2.22.  Plan  view  SEM  image  of  PSiC  prepared  from  n-type  4FI-SiC  by  electroless  etching  in  dilute 
F1F  etchant  for  60  min. 
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Figure  2.23.  Plan  view  SEM  image  of  PSiC  prepared  from  n-type  4Fl-SiC  by  electroless  etching  in  dilute 
HF  etchant  for  60  min. 


Figure  2.24.  Cross-sectional  SEM  image  of  PSiC  prepared  from  n-type  4FFSiC  by  electroless  etching  in 
dilute  F1F  etchant  for  60  min. 
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Figure  2.25.  Plan  view  SEM  image  of  the  interface  between  4H-SiC  (etched,  bottom)  and  6H-SiC 
(planar,  top)  showing  the  crystallographic  dependence  of  the  etching  when  both  polytypes  are  present  in 
the  same  planar  region. 


Figure  2.26.  Plan  view  SEM  image  of  micropipes  at  the  interface  between  4Fl-SiC  (etched,  left)  and  6F1- 
SiC  (planar,  right). 
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Figure  2.27.  Cross-sectional  SEM  view  of  the  4H-SiC  (right,  etched)  undercutting  the  6Fl-SiC  (left, 
planar)  within  the  bulk  of  the  wafer. 

Like  the  6Fl-SiC,  porous  4Fl-SiC  often  showed  interesting  morphologies  in  the  direct  vicinity  of 
micropipes  that  terminated  on  the  free  surface.  One  example  is  the  highly  porous  structure  surrounding  a 
micropipe,  which  then  has  an  abrupt  transition  to  a  less  porous  structure  (Figs.  2.28  and  2.29).  This  once 
again  points  to  the  disruption  of  the  lattice  structure  (which  is  certainly  disturbed  near  a  micropipe) 
playing  a  key  role  in  the  formation  of  porous  structures  in  SiC. 

2.3. 1.3  p-doped  4H-SiC 

The  p-doped  4FFSiC  wafer  was  exposed  to  both  the  dilute  and  concentrated  F1F  etchant  for  a 
variety  of  times,  and  etching  was  attempted  on  both  the  polished  Si  side  as  well  as  the  polished  C  side. 
Flowever,  no  etching  of  the  wafer  was  ever  witnessed.  The  etchant  only  accentuated  the  scratches  and 
introduced  a  few  pits  on  the  surface  (Figs.  2.30  and  2.31). 

2.3. 1.4  Semi-insulating  6H-SiC 

A  small  piece  of  semi-insulating  6FFSiC  was  obtained  from  a  collaborator  (Dr.  Steve  Saddow) 
and  exposed  to  the  dilute  F1F  electroless  etchant.  Etching  of  semi-insulating  material  is  particularly 
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attractive,  since  anodic  etching  techniques  cannot  be  used  with  a  semi-insulating  wafer.  However,  after 
exposure  to  the  electroless  etchants,  only  limited  pitting  was  observed,  and  no  porous  morphology  formed 
(Fig  2.32). 

2.3. 1.5  General  Conclusions:  Electroless  Etching  and  Morphology 

Analysis  of  the  band  diagram  at  the  interface  between  crystalline  SiC  and  Pt  metal  provides 
context  in  which  to  interpret  the  relationship  of  electroless  etching  and  morphology  in  n-doped  SiC,  and 
why  semi-insulating  and  p-doped  wafers  do  not  etch  under  the  electroless  etching  conditions.  This 
behavior  can  be  understood  with  reference  to  the  band  diagram  in  Fig.  2.33.  The  left  hand  side  of  the 
figure  illustrates  the  commonly  encountered  case  of  Schottky  barrier  formation  at  a  metal/n-SiC  interface. 
The  barrier  is  formed  by  the  energy  difference  between  the  Fermi  level  in  the  metal  and  the  first  available 
state  in  SiC,  which  is  located  at  the  conduction  band  edge.  Similarly  there  is  a  barrier  for  hole  injection  in 
p-type  SiC,  as  shown  on  the  right  hand  side  of  Fig.  2.33.  Clearly  band  bending  occurs  in  a  direction 
opposite  to  that  necessary  for  etching  for  both  the  semi-insulating  and  p-doped  SiC,  and  this  explains  why 
electroless  etchants  were  unable  to  make  these  wafers  porous.  The  magnitude  of  the  barrier  for  hole 
injection  could  be  lessened  by  utilizing  a  metal  with  a  larger  work  function  (difference  between  the 
vacuum  level,  E0,  and  the  Fermi  level).  Unfortunately,  Pt  already  has  one  of  the  largest  work  functions 
available,  and  it  is  clearly  desirable  due  to  its  catalytic  activity. 

When  using  HFifUSjOs  etchant,  the  porous  morphology  obtained  is  highly  dependent  upon  the 
crystallographic  quality  of  the  starting  SiC  wafer.  Defects  such  as  dislocations  and  micropipes,  in 
particular,  become  highly  decorated  during  etching,  and  current  technology  still  produces  wafers  of 
relatively  high  (~106  cm'2)  screw  dislocation  density.2’50  Another  example  of  the  role  of  defects  is 
presented  in  Fig.  2.34.  This  cross-sectional  SEM  image  is  obtained  from  a  wafer  that  has  been  cleaved 
through  a  micropipe,  which  is  seen  on  the  right  side  of  the  image,  in  contrast  to  the  bulk  of  the  wafer  on 
the  left  side  of  the  image.  On  the  internal  micropipe  wall,  rows  of  small  (~50  nm)  etch  pits  are  clearly 
aligned  perpendicular  to  the  longitudinal  axis  of  the  micropipe.  The  alignment  of  these  pits  suggests  that 
they  are  due  to  dislocations,  since  dislocations  often  align  relative  to  a  crystalline  axis  in  order  to  reduce 
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Figure  2.28.  Plan  view  SEM  image  of  a  large  defect  and  the  region  of  increased  porosity  (lighter  regions) 
in  the  immediate  vicinity  of  the  defect. 


Figure  2.29.  Plan  view  SEM  image  of  the  wafer  in  Fig.  2.28.  This  close-up  image  shows  the  sharp 
interface  between  the  relatively  weakly  etched  regions  (dark  in  Fig.  2.28)  next  to  the  more  strongly  etched 
regions  (light  in  Fig.  2.28). 
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Figure  2.30.  Plan  view  SEM  image  of  p-doped  4Fl-SiC  (polished  silicon  side)  that  has  been  exposed  to 
the  dilute  F1F  etchant  for  60  min. 


Figure  2.31.  Plan  view  SEM  image  of  p-doped  4FI-SiC  (polished  carbon  side)  that  has  been  exposed  to 
the  dilute  F1F  etchant  for  60  min. 
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Figure  2.32.  Plan  view  SEM  image  of  semi-insulating  6H-SiC  that  has  been  exposed  to  the  dilute  F1F 
etchant  for  60  min. 

overall  stress.  Published  results  are  inconclusive  about  the  presence  of  dislocations  near  micropipes; 
studies  show  both  the  presence  of  dislocations  near  micropipes57’58  as  well  as  dislocation- free  6F1 
morphology  preserved  even  at  the  micropipe  surface.59  Alternatively,  the  etch  pits  in  the  micropipe  could 
form  selectively  on  certain  growth  planes  that  spiral  around  the  micropipe,58’60'61  although  it  is  unknown 
why  the  etch  pits  would  align  only  on  certain  growth  planes.  Regardless  of  the  specific  mechanism 
responsible  for  the  alignment,  the  locations  of  the  etch  pits  are  not  random  and  are  inherently  linked  to  the 
crystallographic  nature  of  the  SiC  substrate. 

2.3.2  Anodically  Etched  SiC 

Numerous  pieces  of  the  6Fl-SiC  wafer  were  anodically  etched,  in  order  to  study  the  effect  of 
etching  parameters  such  as  the  etchant  concentration,  etch  voltage,  and  etch  current.  Representative 
examples  of  different  properties  were  observed  from  the  etch  conditions  listed  in  Table  2.1.  All  of  the 
pieces  listed  in  Table  2.1  were  etched  for  40  min. 
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Figure  2.34.  Cross-sectional  SEM  image  of  a  PSiC  sample  cleaved  through  a  micropipe  (right  side)  and 
the  surrounding  bulk  wafer  (left  side).  Etch  pits  are  observed  at  high  density  on  the  internal  micropipe 
wall,  aligned  perpendicular  to  its  longitudinal  axis.  The  pits  are  much  more  sparse  in  the  bulk  portion  of 
the  wafer. 


Table  2.1.  PSiC  Etching  Conditions 


Sample  # 
1 
2 

3 

4 

5 

6 


Etch  Voltage  (Y) 

3 

4 

4 

5 
4 
3 


Etch  Current  (mA/cnr) 
80 
40 
35 
35 
20 
10 


Porous  Layer  Thickness  (urn) 

3 

4 
2 
7 
10 
3 


The  top  surface  of  each  anodically  etched  PSiC  wafer  is  composed  of  nanofingers  of  SiC 
arranged  in  random  fashion  (Fig.  2.35).  In  cross-section,  the  PSiC  layers  consistently  show  a  propensity 
for  triangular  pitting.  In  most  instances,  the  triangular  etch  pits  do  not  follow  a  regular  pattern,  and, 
therefore,  a  highly  branched  structure  results  (Figs.  2.36  and  2.37).  The  PSiC  nanocrystallites  have  an 
average  diameter  of  approximately  25  nm,  regardless  of  the  forward  bias  (3-5  V)  or  current  density  (10- 
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80  mA/cm2)  used  to  etch  each  wafer.  The  morphology  is  striking  in  that  the  size  of  the  pores  seen  in 
cross-section  is  relatively  monodisperse. 

In  one  PSiC  sample  (sample  #4),  the  triangular  etch  pits  were  observed  to  orient  regularly  upon 
planes  parallel  to  the  top  free  surface  of  the  wafer  (Fig.  2.38).  This  pitting  structure  is  reminiscent  of  the 
triangular  etch  pits  seen  within  the  micropipe  in  PSiC  produced  via  electroless  etching  (Fig.  2.34). 

Sample  #4  was  prepared  from  the  same  wafer  as  the  other  PSiC  samples  discussed  in  this  work,  so  the 
cause  of  this  particular  morphology  is  unknown. 

2.4  BET  Surface  Area  Analysis 

The  high  surface  area  characteristic  of  porous  materials  is  one  of  the  drivers  for  their  application 
in  sensing  platforms.  Surface  areas  for  several  PSiC  samples  formed  via  electroless  etching  were 
measured  via  BET  isotherms  to  be  -150  m2/cm3,  which  is  comparable  to  the  surface  areas  obtained  in 
anodically  etched  PSi.62  All  of  the  anodically  etched  pieces  in  this  work  were  limited  to  10  pm  deep  or 
less.  Considering  that  all  etched  pieces  were  -5x5  mm  in  size  or  less,  no  single  piece  of  anodically 
etched  material  had  a  substantial  amount  of  total  surface  area  for  accurate  adsoiption  measurements.  For 
this  reason,  accurate  surface  area  analysis  via  BET  isotherms  was  never  obtained  for  any  of  the  anodically 
etched  pieces  produced  in  this  work. 

2.5  Conclusions 

Electroless  etching  of  SiC  in  HF  and  S2082'  is  competent  to  form  PSiC,  although  reproducible 
formation  of  deep  porous  structures  remains  an  unsolved  issue.  The  spatial  extent  of  porosity  in  PSiC  is 
dependent  on  the  polytype  and  on  the  crystallographic  perfection  of  the  wafer.  Furthermore,  for  6H-SiC 
the  resulting  morphology  depends  on  the  etchant  concentrations.  Surface  areas  on  the  order  of  150 
m2/cm3  can  be  produced  easily,  and  in  6H-SiC  an  extensively  interconnected  network  of  pores  can  be 
achieved  to  significant  depths  (>  100  pm).  These  factors  argue  for  the  efficient  transport  of  reactants  and 
products  within  the  interconnected  network  and  imply  that  these  materials  are  suitable  for  use  as  chemical 
and  physical  sensors. 
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Figure  2.35.  Plan  view  SEM  image  of  PSiC  prepared  from  n-type  6Fl-SiC  by  anodic  etching  for  40  min. 


200nm 


Figure  2.36.  Cross-sectional  SEM  image  of  PSiC  prepared  from  n-type  6Fl-SiC  by  anodic  etching  for  40 
min. 
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Figure  2.37 .  Cross-sectional  SEM  image  of  PSiC  prepared  from  n-type  6Fl-SiC  by  anodic  etching  for  40 
min. 
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Figure  2.38.  Cross-sectional  SEM  image  of  PSiC  prepared  from  n-type  6FI-SiC  by  anodic  etching  for  40 
min. 
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Anodic  etching  of  SiC  yielded  morphologies  similar  to  electroless  etching.  However,  greater 
control  of  the  depth  of  porous  etching  is  possible  with  anodic  etching.  For  this  reason,  anodic  etching 
was  used  in  the  production  of  PSiC  samples  tested  for  their  sensing  capabilities  (Chapter  4). 
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CHAPTER  3 


SPECTROSCOPIC  CHARACTERIZATION 


3.1  Introduction 

In  addition  to  sensors  and  substrates  for  epitaxial  growth,  another  application  that  has  received 
widespread  attention  is  the  possibility  of  using  porous  semiconductors  as  light  emitting  diodes  (LEDs). 
Research  into  porous  semiconductors  grew  exponentially  after  the  seminal  discovery  by  Canham  that 
porous  silicon  is  capable  of  light  emission  at  photon  energies  in  excess  of  the  bulk  bandgap.1  The  blue- 
shifted  emission  from  PSi  is  apparently  tunable  via  manipulation  of  the  etching  conditions.2"5  The 
intensity  of  the  light  emission  from  PSi  is  often  reported  to  be  as  much  as  500  times  more  intense  than 
emission  from  crystalline  Si.6’7 

Quantum  confinement  effects  are  one  plausible  explanation  for  the  blue-shifted  emission.  In  an 
indirect  bandgap  semiconductor,  bandgap  transitions  are  dipole  forbidden  due  to  translational  symmetry. 
However,  when  a  small  nanocrystallite  is  formed,  a  supeiposition  of  eigenstates  with  a  distribution  of  k 
values  is  necessary  to  physically  localize  the  wavefunction  inside  the  nanocrystal.  In  turn,  physical 
localization  breaks  the  translational  symmetry,  allowing  for  the  mixing  of  eigenstates  with  different  k 
values.  As  a  result,  the  combination  of  a  relaxation  of  the  momentum  selection  rules  as  well  as  distortions 
of  the  local  band  structure  can  allow  for  direct  bandgap  transitions  in  nanocrystallites  of  indirect  gap  Si. 
However,  researchers  have  encountered  difficulties  in  producing  a  PSi  LED.  Particularly,  emission  from 
the  shorter  (blue)  wavelengths  has  proven  to  be  unstable  and  still  relatively  weak. 

The  successes  with  silicon  naturally  led  investigators  to  question  whether  comparably  enhanced 
properties  could  be  coaxed  from  porous  analogs  of  wide  bandgap  materials,  such  as  silicon  carbide  and 
gallium  nitride.  Because  silicon  carbide  is  an  indirect  bandgap  material,  the  luminescence  of  SiC  LEDs  is 
very  low  (quantum  efficiency  only  ~10"4),  an  efficiency  unsuitable  in  applications  such  as  LED  displays. 
However,  if  PSiC  could  be  developed  to  exhibit  properties  similar  to  those  of  porous  silicon,  a  bright  UV 
LED  could  be  realized.8"11 
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In  addition  to  exploring  the  possibility  of  producing  UV  LEDs,  spectroscopic  studies  provide  a 
tool  to  characterize  the  porous  nanostructure.  Cathodoluminescence  and  photoluminescence  spectra  of 
PSiC  were  collected  in  this  work  in  order  to  look  for  a  possible  blue-shift  in  emission  due  either  to 
quantum  confinement  effects  or  to  the  formation  of  new  emissive  surface  states  upon  etching.  The  exact 
peak  location  for  the  bandgap  of  SiC  shifts  slightly  from  wafer  to  wafer  due  to  differences  in  doping;  for 
the  samples  examined  in  this  work,  the  peak  cathodoluminescence  (CL)  bandgap  emission  occurred  at 
426  nm  (6H-SiC)  or  392  nm  (4H-SiC). 

Few  observations  of  CL  emission  from  PSiC  layers  have  been  reported,  and  reported  observations 
of  UV  emission  have  been  inconsistent.  Shor  et.  al.12'13  reported  a  UV  CL  emission  peak  at  335  nm, 
which  was  assigned  to  the  native  oxide  or  possibly  quantum  confinement  effects.  Danishevskii  et.  al.14 
reported  evidence  of  increased  CL  emission  centered  at  468  nm,  but  no  UV  emission. 

PL  emission  from  crystalline  SiC  has  been  well  documented  in  the  literature.  Depending  upon 
the  doping  concentrations,  a  broad  peak  appearing  between  580-620  nm  has  been  attributed  to  a 
combination  of  N'-B+  recombination  and  “defect”  luminescence.1517  When  aluminum  impurities  are 
present  in  the  wafer,  a  broad  peak  centered  near  470  nm  may  be  seen  and  is  assigned  to  N'-A1+ 
recombination.15 

A  variety  of  photoluminescence  results  have  been  published  for  anodically  etched  PSiC.  There 
are  several  reports  in  the  literature  of  an  increase  in  photoluminescence  intensity  in  the  460-500  nm 
region  of  the  spectrum  upon  formation  of  PSiC.15'23  Increased  emission  in  this  wavelength  range  has  been 
reported  for  both  porous  material  prepared  from  both  4H  and  6H  polytypes  of  SiC;  PL  emission  does  not 
appear  to  be  polytype  dependent  when  the  SiC  is  made  porous  by  anodic  etching.  In  contrast,  a  few 
authors  report  that  formation  of  PSiC  results  in  relatively  small  or  no  shift  in  the  PL  of  the  wide  visible 
band  attributed  to  defects  and  N'-B  recombination.  ’  No  published  research  has  indicated  PL 
emission  in  the  ultraviolet  region  (blue  of  the  bandgap  emission).  The  bandgap  emission  itself  is  difficult 
to  see  in  PL  experiments  due  to  the  indirect  nature  of  the  bandgap  in  SiC. 
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The  origin  of  the  large  blue  shift  in  the  visible  photo  luminescence  of  PSiC  with  respect  to  bulk 
SiC  is  not  yet  firmly  established.  Two  categories  of  mechanisms,  categorized  as  extrinsic  or  intrinsic  in 
nature,  have  been  invoked  to  explain  this  PL. 

There  are  many  proposed  extrinsic  changes  that  may  occur  in  PSiC  to  cause  blue-shifted  PL 
emission.  For  example,  changes  in  the  Raman  spectra  led  Danishevskii  et.  al.  to  report  the  formation  of 
3C-SiC  (the  cubic  polytype)  nanocrystallites  within  the  hexagonal  matrix  of  the  wafer.26  Parkhutik  et.  al. 
theorize  the  formation  of  Si  nanocrystallites  within  the  SiC  matrix.27 

A  different  train  of  thought  involves  the  formation  of  some  type  of  surface  hydride  species  or 
some  other  change  to  the  surface  of  the  SiC  crystal.  Researchers  who  have  witnessed  altered  PL  emission 
characteristics  upon  oxidation  or  hydrogenation  of  the  PSiC  surface  often  propose  this  argument.  Surface 
state  impurities  proposed  as  causes  of  blue-shifted  PL  include  species  such  as  hydrogenated  amorphous 
silicon  carbide.8  The  passivation  of  dangling  bonds  has  also  been  theorized  as  the  cause  of  blue  PL,  since 
dangling  bonds  can  serve  as  nonradiative  recombination  centers  (localized  states)  in  the  bandgap.28  Some 
researchers  claim  that  etching  results  in  an  increase  of  C-H  terminations  on  the  surface24'29'30  or  Si-H 
terminations  on  the  surface.31  Finally,  some  researchers  have  proposed  that  etching  causes  an  increase  of 
defects  on  the  surface,  and  these  defects  serve  as  recombination  centers  for  new  PL.  ’  ’  ’  ’ 

However,  the  blue-shifted  PL  may  also  be  due  to  an  intrinsic  change:  the  development  of 
topographically  altered  structures  which  can  support  quantum  confinement  effects.  Distortion  of  the  local 
band  structure  could  lead  to  the  relaxation  of  momentum  selection  rules,  thereby  allowing  direct 
transitions  in  the  semiconductor  material  and,  hence,  an  increase  in  emission  efficiency.  The  apparent 
tunability  of  the  wavelength  of  maximum  PL  emission  with  varying  etching  conditions  points  to  a 
quantum  confinement  effect  explanation.  However,  to  exhibit  quantum  confinement,  the  widths  of  the 
nanocrystallites  must  be  significantly  less  than  the  dimensions  of  the  free  exciton.  The  free  exciton  is  ~50 
A  in  Si;3  the  free  exciton  size  for  SiC  was  not  located  in  the  literature,  but  reports  in  the  literature  indicate 
that  spherical  nanocrystallites  in  PSiC  must  also  have  diameters  <50  A  to  show  significant  increases  in 
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bandgap  emission.  ’  ’  ’  These  small  dimensions  were  not  observed  in  some  of  the  samples  that 
displayed  blue-shifted  PL  in  the  literature.21 

Evidence  for  a  quantum  confinement  effect  in  PSiC  can  be  revealed  by  Raman  scattering 
measurements,  as  well.  Formation  of  a  porous  nanostructure  leads  to  broadening  and  shifting  of  Raman 
bands,  because  the  translational  lattice  symmetry  is  relaxed,  and  coupling  of  phonon  states  with  different 
k  values  occurs.  Previous  work2'37'42  has  shown  that  phonon  confinement  effects  cause  the  longitudinal 
optic  (LO)  mode  of  PSi  to  broaden  asymmetrically  to  the  low-frequency  side  and  red-shift  as 
nanostructures  are  formed.  Since  the  morphology  investigations  reported  in  Chapter  2  clearly  show 
evidence  of  nanostructure  formation,  Raman  investigations  were  conducted  in  an  attempt  to  quantify  the 
degree  of  porosity  non-destructively.  However,  Raman  data  must  be  interpreted  carefully,  as  other 
factors  such  as  changes  in  the  doping  levels43'44  and  residual  stress  in  the  crystal45'48  can  also  lead  to 
distortion  of  the  Raman  bands.  HXRD  studies,  detailed  in  Chapter  2,  provide  evidence  that  the  porous 
layers  formed  in  this  work  are  under  tensile  stress,  and  this  stress  can  distort  and  shift  the  Raman  bands,  if 
the  stress  is  great  enough. 

3.2  Cathodoluminescence  Spectra 

Cathodoluminescence  spectra  were  obtained  at  beam  energies  of  1  to  20  keV  with  a  Zeiss  DSM 
960  SEM  utilizing  an  Oxford  Instruments  MonoCL  monochromator  and  photomultiplier  tube  detector 
(Amherst  Scientific). 

3.2.1  Electrolessly  Etched  SiC 

CL  spectra  obtained  from  porous  material  produced  from  both  6H-SiC  and  4H-SiC  do  not  show 
any  blue-shift  of  the  bandgap  luminescence  nor  any  other  changes  in  the  visible  portions  of  the  spectra 
(Figs.  3.1  and  3.2).  6H-SiC  spectra  were  dominated  by  the  bandgap  emission  at  426  nm  and  a  broad  peak 
near  600  nm  due  to  donor-acceptor  recombination.  4H-SiC  spectra  exhibited  a  sharp  bandgap  peak  at  392 
nm  and  a  broad  donor-acceptor  recombination  peak  centered  near  480  nm. 
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Figure  3.1.  Cathodoluminescence  spectra  of  unetched  6H-SiC  (dashed  line)  and  PSiC  prepared  from  n- 
type  6H-SiC  via  electroless  etching  in  dilute  F1F  etchant  for  60  min  (solid  line). 


Figure  3.2.  Cathodoluminescence  spectra  of  unetched  4Fl-SiC  (dashed  line)  and  PSiC  prepared  from  n- 
type  4Fl-SiC  via  electroless  etching  in  dilute  F1F  etchant  for  60  min  (solid  line). 
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3.2.2  Anodically  Etched  SiC 

CL  spectra  of  anodically  etched  PSiC  samples  produced  in  this  work  consistently  show  peak 
emission  blue  of  the  bandgap  emission  at  426  nm  (Fig.  3.3).  CL  emission  spectra  exhibited  wavelengths 
of  maximum  emission  which  varied  between  355  and  400  nm.  Higher  etching  current  densities  yielded 
bluer  emission,  in  agreement  with  results  published  elsewhere.13  The  increase  in  UV  CL  intensity  seems 
to  be  accompanied  by  loss  of  emission  intensity  at  590  nm  (the  peak  widely  attributed  to  donor-acceptor 
(N'-B+)  recombination  and  defects).  In  addition,  the  UV  peak  can  become  so  intense  that  it  entirely 
masks  any  bandgap  CL  emission,  cf.  spectra  of  samples  #4  and  #5  in  Fig.  3.3. 

Special  care  was  taken  to  examine  the  UV  emission  evident  in  Fig.  3.3,  since  the  observation  of 
significant  emission  intensity  at  300  nm  is  surprising.  In  particular  spectra  were  collected  several  times 
with  a  variety  of  slit  widths  between  the  sample  and  monochromator  and  the  monochromator  and 
photodetector.  Even  for  the  smallest  collection  slit  widths  (10  pm),  the  UV  emission  remained 
measurable  (although  weak)  at  300  nm. 

CL  spectra  were  obtained  at  1,  5,  10,  15,  and  20  keV  in  order  to  observe  any  difference  in  the 
spectra  as  a  function  of  depth  within  the  wafer.  Utilizing  the  work  of  Kanaya  and  Okayama,49  50  we  can 
predict  that  accelerating  voltages  of  1-20  keV  correspond  to  total  sampling  depths  of  0.24-3.6  pm,  with 
maximum  emission  coming  from  depths  of  0.007-1.1  pm  within  the  porous  layer.5152  However,  no  shifts 
in  the  CL  spectra  were  observed  for  the  different  voltages;  the  weaker  excitation  voltages  resulted  simply 
in  weaker  signals.  This  indicates  that  the  UV  CL  emission  is  not  restricted  to  the  top  surface  layer  of  the 
porous  network. 

3.3  Photoluminescence  Spectra 

Photoluminescence  spectra  were  excited  at  A  =  325.0  nm  (HeCd  laser).  A  single  grating  (150 
groove/mm)  was  used  in  the  spectrograph  stage  for  wide  spectral  coverage,  and  the  filter  stage  of  the 
spectrometer  was  not  used. 
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Many  reports  in  the  literature  quantify  increases  in  PL  emission  intensity  (up  to  500  times  greater 
than  the  emission  from  crystalline  SiC)  upon  the  formation  of  PSiC.  However,  it  is  likely  that  much  of 
this  increased  intensity  is  due  to  increased  outcoupling  efficiency  attributable  to  the  rough  surface  formed 
by  etching,  and  any  inherent  increases  in  PL  intensity  due  to  the  porous  network  are  difficult  to  isolate.  In 
fact,  in  this  work  for  the  325  nm  laser,  the  unetched  PL  data  was  actually  obtained  from  the  back  side  of 
the  wafer,  since  it  was  not  polished  and  therefore  yielded  a  stronger  signal.  No  difference  was  observed 
between  PL  emission  from  the  front  and  back  surfaces  except  for  the  increase  in  intensity  from  the  back 
surface.  In  this  work,  no  comparisons  of  PL  emission  intensity  will  be  made  between  samples. 

3.3.1  Electrolessly  Etched  SiC 

In  contrast  to  PL  results  in  the  literature  for  anodically  etched  PSiC,  the  electrolessly  etched  PSiC 
produced  in  this  work  never  shows  maximum  PL  emission  in  the  460-500  nm  range.  However,  all 
spectra  of  PSiC  show  increased  PL  emission  below  500  nm.  In  some  samples  (Fig.  3.4,  top)  this 
increased  emission  is  pronounced  with  an  unresolved  emission  band  peaking  in  the  440-460  nm  region. 
Changes  in  the  PL  of  the  broad  visible  band  are  observed,  too.  As  shown  in  Fig.  3.4  the  wavelength  of 
maximum  emission  in  PL  blue-shifts  8-20  nm  in  the  porous  material. 

3.3.2  Anodically  Etched  SiC 

The  PL  spectra  collected  from  anodically  etched  PSiC  shows  wide  variation  in  the  wavelength  of 
maximum  emission  (Fig.  3.5).  However,  in  contrast  to  the  CL  spectra,  maximum  PL  emission  is  always 
sub-bandgap  in  energy.  Maximum  PL  emission  occurs  anywhere  from  550  nm  (near  the  N'-ET 
recombination  peak)  to  430  nm  (near  bandgap  emission).  Low  power  density  (15  mW/mnr)  was  utilized 
to  collect  the  spectra,  ensuring  that  the  high  energy  peaks  really  are  due  to  PL  and  not  bandgap  emission. 
No  apparent  relation  between  anodization  voltage/current  and  resulting  PL  emission  was  observed,  in 
contrast  to  some  results  in  the  literature.18,23 

Because  the  emission  in  this  spectral  region  is  below  the  bandgap  energy,  the  most  likely  source 
of  luminescence  is  band-tail  states,22  the  blue  shift  upon  etching  signifying  a  shift  in  the  distribution  of 
such  states.  When  the  power  of  the  excitation  laser  was  decreased,  the  maximum  PL  emission  from 
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Figure  3.3.  Cathodoluminescence  spectra  of  unetched  6Fl-SiC  and  three  different  pieces  of PSiC  prepared 
from  n-type  6Fl-SiC  by  anodic  etching  for  40  min.  Sample  numbers  correspond  to  the  notation  given  in 
Table  2. 1  of  Chapter  2. 
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Figure  3.4.  Photoluminescence  spectra  of  unetched  6Fl-SiC  (dashed  line)  and  two  pieces  of  PSiC 
prepared  from  n-type  6Fl-SiC  via  electroless  etching  in  dilute  F1F  etchant  for  60  min. 
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Figure  3.5.  Photo  luminescence  spectra  of  unetched  6H-SiC  and  three  different  samples  of  PSiC  prepared 
from  n-type  6H-SiC  via  anodic  etching  for  40  min.  Sample  numbers  correspond  to  the  notation  given  in 
Table  2. 1  of  Chapter  2. 

sample  #3  was  found  to  blue  shift.  This  strongly  suggests  that  the  maximum  PL  emission  for  this  piece  of 
PSiC  arises  from  donor-acceptor  recombination.53,54  However,  the  wavelength  of  maximum  emission 
from  sample  #4  does  not  shift  as  the  power  of  the  excitation  laser  beam  is  varied.  One  possible 
conclusion  is  that  the  spectrum  from  sample  #4  is  dominated  by  defect  emission. 

The  PL  data  collected  with  325  nm  excitation  give  broad  emission,  which  tail  into  the  UV,  but 
they  do  not  show  the  UV-dominated  emission  seen  in  the  CL  spectra.  The  differences  in  the  PL  and  CL 
spectra  are  probably  not  a  result  of  different  probe  depths,  because  the  CL  and  PL  sample  similar  volumes 
(3.6  pm  deep  in  CL  vs.  4.3  pm  deep  in  PL).  This  suggests  that  the  high  energy  electrons  and  photons 
used  in  CL  and  PL  respectively  access  different  emissive  states.  For  example,  the  UV  emission  witnessed 
in  CL  may  be  some  type  of  defect  state  that  is  accessed  by  the  20  kV  electrons  utilized  in  CL,  but  the  state 
cannot  be  accessed  by  the  4  eV  photons  used  to  collect  the  PL  spectra. 
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3.4  Raman  Spectra 

Raman  spectra  were  excited  at  three  different  wavelengths.  Initially,  spectra  were  excited  with 
sub-bandgap  radiation  (Ar+  laser,  A  =  457.9  nm)  and  collected  on  a  single  monochromator  (Spex  500M, 
1200  groove/mm  grating)  equipped  with  a  liquid  N2  cooled  charge  coupled  device  (CCD)  camera 
(Photometries  Ltd.).  A  holographic  notch  filter  (Kaiser)  was  used  to  remove  scattered  laser  line  from  the 
Raman  scattered  radiation  (Fig.  3.6).  This  wavelength  was  used  to  map  the  polytype  of  the  4H-SiC  wafer 
as  well  as  initially  characterize  the  porous  layers. 

Later,  Raman  spectra  were  excited  with  a  HeCd  laser  operating  at  441.6  nm  as  well  as  a  different 
HeCd  laser  emitting  at  325.0  nm.  The  sub-bandgap  441.6  nm  excitation  samples  the  entire  wafer 
thickness,55  so  data  from  these  experiments  reflect  the  properties  of  the  PSiC  film  as  well  as  the 
crystalline  SiC  substrate.  Laser  excitation  at  325.0  nm  (above  bandgap  energy)  was  utilized  in  order  to 
probe  the  near- surface  region  of  the  wafer  preferentially,  as  this  wavelength  has  a  probe  depth  of  ~4.3 
pm.56  These  Raman  spectra  were  collected  with  a  triple  spectrometer  (Spex  1877,  0.6  m)  with  two  2400 
groove/mm  gratings  in  the  filter  stage,  a  2400  groove/mm  grating  in  the  spectrograph  stage,  and  a  UV- 
sensitive  CCD  camera  (Roper  Scientific)  (Fig.  3.7).  In  SiC,  the  Ai  LO  mode  at  966  cm1  is  a  polar  mode 
in  C6v  symmetry  and,  thus,  couples  very  strongly  to  the  macroscopic  electric  field  accompanying  the  LO 
phonon,  distorting  its  intensity  relative  to  the  TO  modes  when  above-bandgap  excitation  is  used.47’57,58  It 
is  the  LO  mode,  then,  that  was  primarily  monitored  for  changes  indicative  of  quantum  confinement 
effects. 

3.4.1  Electrolessly  Etched  SiC 

Initially  spectra  were  excited  at  457.9  nm,  and  this  wavelength  provides  a  strong  signal  from  the 
PSiC  wafer.  The  spectra  are  dominated  by  the  expected  optic  branch  phonons  allowed  in  backscattering 
along  the  z-axis  (A!  LO  and  E2  TO)  (Figs.  3.8  and  3.9).  Because  the  laser  energy  is  sub-bandgap,  the 
signal  originates  from  the  crystalline  SiC  substrate  as  well  as  the  much  thinner  PSiC  layer,  making  the 
spectra  of  the  etched  and  unetched  samples  indistinguishable. 
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Figure  3.6.  Experimental  arrangement  for  visible  Raman  collection  at  457.9  nm. 
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Figure  3.7.  Experimental  arrangement  for  UV  Raman  and  PL  collection. 
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Figure  3.8.  Unpolarized  Raman  spectra  of  etched  (top,  solid)  and  unetched  (bottom,  dashed)  6H-SiC  in 
the  backscattering  configuration  excited  at  457.9  nm. 
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Figure  3.9.  Unpolarized  Raman  spectra  of  etched  (top,  solid)  and  unetched  (bottom,  dashed)  4H-SiC  in 
the  backscattering  configuration  excited  at  457.9  nm. 
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The  457.9  nm  excitation  was  useful  for  determining  the  polytype  in  the  nominally  4H-SiC  wafer. 
CL  measurements  had  shown  that  some  portions  of  the  4H-SiC  wafer  exhibited  bandgap  values  indicative 
of  the  6H  poly  type.  Therefore,  Raman  spectra  were  collected  at  75  different  points  over  the  surface  of  the 
2”  wafer  in  order  to  map  the  different  polytypes  within  the  wafer.  Raman  mapping  showed  that  a 
majority  of  the  wafer  was  4H,  although  significant  portions  of  6H  material  existed  (Fig.  3.10).  The 
interfaces  between  the  414  and  614  regions  often  showed  characteristics  of  each  polytype  when  examined 
via  Raman  spectroscopy  (Fig.  3.11).  The  614  wafer  was  also  mapped  and  found  to  be  all  614  in  polytype. 

It  is  interesting  that  the  Raman  modes  in  Fig.  3.11  show  peak  locations  between  that  of  pure  414 
and  pure  614  character,  given  the  definite  interface  seen  between  the  two  polytypes  in  cross-sectional 
SEM  (Fig  2.27).  The  Raman  results  clearly  indicate  a  strained  interface  that  would  not  be  apparent  from 
the  SEM  images  alone. 

Subsequent  to  the  457.9  nm  Raman  studies,  an  improved  optical  system  was  constructed  within 
our  laboratory.  Initially,  Raman  spectra  were  again  collected  with  sub-bandgap  excitation  (441.6  nm, 
HeCd  laser)  in  order  to  analyze  the  spectra  with  improved  resolution.  However,  once  again  no 
broadening  of  the  modes  was  witnessed  for  the  porous  SiC  samples. 

To  circumvent  problems  with  penetration  depth,  Raman  data  were  then  collected  with  325.0  nm 
excitation,  a  photon  energy  significantly  above  the  300  K  bandgap  energy.  This  above-bandgap  laser  line 
probes  a  much  smaller  volume  of  the  PSiC  at  the  expense  of  lower  signal  strength.  However,  even  with 
325.0  nm  excitation,  the  LO  mode  was  not  found  to  shift  or  broaden  with  the  development  of  porosity. 
The  lack  of  change  in  the  LO  lineshape  with  etching  means  that  either  phonon  confinement  effects  are 
inherently  small  in  electrolessly  etched  PSiC,  or  the  fraction  of  unit  cells  sampled  for  which  confinement 
effects  are  important  is  small.  Additionally  these  results  point  to  the  extreme  stability  of  PSiC,  and 
indicate  that  the  porous  layer  does  not  introduce  significant  strain. 

Unfortunately,  while  collecting  Raman  data  with  the  325  nm  laser  only,  a  “band”  always  appears 
near  854  cm"1.  The  cause  of  this  band  is  unknown,  since  it  does  not  equate  to  any  mode  in  our  polytypes 
of  SiC  (or  a  multiple  of  a  mode).  It  also  does  not  correspond  to  any  expected  unwanted  emission  from  the 
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Figure  3.10.  Polytype  map  of  the  2”  diameter  wafer  that  was  purchased  as  pure  4Fl-SiC.  The  blue 
regions  were  found  to  be  4Fl-SiC;  the  white  regions  were  found  to  be  61 1-SiC. 


Raman  Shift  (cm'1) 

Figure  3.11.  Unpolarized  Raman  spectra  of  pure  6Fl-SiC  (top),  pure  414-SiC  (middle),  and  a  wafer  that 
has  both  polytypes  present  (bottom).  The  spectra  were  excited  at  457.9  nm  and  collected  in  the 
backscattering  configuration. 
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HeCd  laser,  and  the  use  of  two  prisms  in  the  optical  system  did  not  remove  this  line.  The  854  cm"1  band 
could  originate  from  the  15R  polytype  of  SiC,  but  this  explanation  seems  improbable,  since  no  other 
evidence  for  the  15R  polytype  was  observed.  The  854  cm'1  band  is  most  likely  a  spectral  artifact  from  the 
1200  groove/mm  gratings,  and  therefore  unfortunately  appears  in  all  325  nm  Raman  spectra  presented  in 
this  work. 

3.4.2  Anodically  Etched  SiC 

Raman  spectra  collected  with  441.6  nm  excitation  (sub-bandgap)  were  unchanged  between  the 
crystalline  SiC  substrate  and  the  anodically  etched  PSiC  layers.  However,  Raman  spectra  collected  with 
325  nm  excitation  (above  bandgap)  showed  broadening  of  the  longitudinal  optic  (LO)  mode  to  lower 
wavenumbers;  the  degree  of  broadening  varied  greatly  between  samples  (Fig.  3.12).  The  peak  in  the  LO 
Raman  band  in  PSiC  showed  little  (<  3  cm'1)  or  no  shifting  to  lower  wavenumbers. 

Similar  broadening  has  been  reported  by  other  authors  in  PSiC18  26  and  attributed  to  confinement 
effects.  One  model  used  to  describe  Raman  lineshapes38,39  is 


The  model  in  Eq.  3.1  assumes  spherical  nanocrystallites  of  diameter  L.  In  Eq.  3.1,  co  is  the  Raman 
frequency,  q  is  the  Brillouin  zone  dimension  (2n la,  where  a  =  1.51 17  nm  for  6H-SiC),  and  T0  is  the  line 
width  of  the  crystalline  material.  co(q)  is  the  phonon  dispersion  relation,  for  which  we  used 

f  \ 

co2 (q)=  657, 500  +  27 0,000 cox  ^  .  (3.2) 

This  phonon  dispersion  relation  was  chosen  after  fitting  the  phonon  dispersion  curve  for  all  of  the  SiC 
polytypes  to  a  best-fit  function  (Fig.  3. 13). 47,58-64 

A  Mathematica  program  was  written  by  Dr.  Lindesay  at  Howard  University  in  order  to  fit  the 
model  given  in  Eq.  3.1  to  the  experimental  data  (see  Appendix  A).  Fitting  Eq.  3.1  to  our  Raman  spectra 
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predicts  nanocrystallite  diameters  ranging  from  6-25  nm.  However,  the  theoretical  line  shapes  do  not  fit 
the  experimental  data  well  (Fig.  3.14).  Eq.  3.1  predicts  that  broadening  of  the  LO  mode  should  be 
accompanied  by  a  significant  shift  in  the  peak  location  of  the  LO  mode  to  lower  wavenumbers.  However, 
as  seen  in  Fig.  3.12,  significant  broadening  of  the  LO  mode  is  accompanied  by  only  relatively  small  shifts 
to  lower  wavenumbers. 

Therefore,  a  bimodal  fit  was  implemented  in  an  attempt  to  account  for  the  sampling  of  some  large 
crystallites  (effectively  bulk  SiC)  as  well  as  nanocrystallites  of  some  small,  average  diameter.  Again,  the 
experimental  LO  mode  is  too  broad  to  be  effectively  represented  by  only  a  bimodal  distribution  of 
nanocrystalline  sizes  (Fig.  3.15).  In  actuality,  of  course,  the  PSiC  layers  are  probably  composed  of  a 
continuum  of  nanocrystallite  sizes,  and  the  Raman  spectra  would  be  best  fit  by  a  wide  distribution  of 
sizes,  if,  in  fact,  phonon  confinement  is  the  cause  of  the  LO  mode  broadening.  However,  even  when 
including  nanocrystallite  dimensions  on  the  order  of  just  a  few  atoms  in  diameter,  the  quantum 
confinement  equation  cannot  properly  model  the  wide  broadening  seen  in  the  experimental  Raman 
spectra. 

Another  possible  explanation  for  the  broadening  of  the  Raman  LO  mode  in  polar  materials,  such 
as  PSiC,  is  the  emergence  of  surface  Frohlich  modes.  Literature  reports  of  PSiC  LO  mode  broadening  26 
very  similar  to  the  results  found  in  this  work  have  been  attributed  to  the  longitudinal-transverse  splitting 
of  a  Frohlich  mode.65  Frohlich-type  modes  have  been  reported  in  porous  GaN;57  in  addition,  LO  mode 
broadening  attributed  to  surface  phonons  has  been  found  in  GaP  microcrystals66  and  porous  GaP.67 
Finally,  infrared  reflectance  spectra  of  PSiC  films  also  show  indications  of  surface  phonon  modes.68 

The  Raman  spectra  in  this  work  were  fit  to  a  bimodal  Gaussian  distribution  consisting  of  the  bulk 
SiC  LO  phonon  plus  a  peak  due  to  the  surface  Frohlich  modes.  Following  the  model  of  Sarua  et.  al.,65  the 
broadening  of  the  LO  mode  in  PSiC  was  assumed  to  be  due  to  the  combinatorial  effect  of  extraordinary 
LO  modes,  longitudinal  Frohlich  modes,  and  transverse  Frohlich  modes.  The  intensities  from  these  three 
types  of  modes  were  weighted  for  the  orientation  of  the  incident  laser  beam  used  in  this  work.  The  results 
of  this  model  must  be  used  with  caution;  in  our  case,  the  model  yields  information  of  a  qualitative  nature 
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Figure  3.12.  Unpolarized  Raman  spectra  (collected  at  325.0  nm)  of  unetched  6Fl-SiC  (top)  and  two 
pieces  of  PSiC  prepared  via  anodic  etching  for  40  min.  The  spectra  were  excited  at  325.0  nm  and 
collected  in  the  backscattering  configuration. 


Figure  3.13.  Dispersion  curve  in  the  I  L  direction  for  SiC  showing  the  experimentally  determined  LO 
modes  found  in  the  literature  (points)  and  the  best  fit  line  in  accordance  with  Eq.  3.2. 
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Figure  3.14.  Quantum  confinement  model  fit  (dashed)  compared  to  the  experimentally-determined  LO 
mode  (solid)  for  an  anodically  etched  piece  of  614-SiC  (sample  #5). 
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Figure  3.15.  Experimentally-determined  LO  mode  (solid)  for  an  anodically  etched  piece  of  6Fl-SiC 
(sample  #5)  compared  to  the  theoretical  curve  combining  a  Gaussian  curve  for  crystalline  6Fl-SiC 
combined  with  a  quantum  confinement  peak  (dashed). 


65 


at  best.  The  properties  of  Frohlich  modes  strongly  depend  upon  the  morphology  (symmetry,  structural 
geometry,  effective  size,  etc.)  and  on  the  properties  of  the  constituents.  The  model  assumes  the  porous 
material  is  columnar  in  structure;  however,  as  our  SEM  images  show,  our  porous  material  is  not  as 
columnar  as,  for  example,  porous  GaN.  In  addition,  this  model  uses  an  effective  medium  approximation, 
which  is  not  strictly  valid  for  our  spectra,  since  the  wavelength  of  the  laser  is  only  moderately  greater  than 
the  dimensions  of  the  porous  material.  As  a  result,  we  can  expect  that  our  Raman  spectra  include  the 
effects  of  strong  diffuse  scattering,  accompanied  by  a  broadening  of  spectral  lines  and  violation  of  the 
symmetry  selection  rules. 

Following  the  model  of  Sarua  et.  al.,65  the  total  Frohlich  mode  was  assumed  to  be  composed  of 
three  distinct  components:  a  longitudinal  Frohlich  character,  a  transverse  Frohlich  character,  and 
extraordinary  LO  modes.  The  longitudinal  Frohlich  and  extraordinary  LO  modes  have  symmetry  along 
the  c-axis  of  the  crystal;  the  transverse  Frohlich  mode  is  in  the  plane  of  the  crystal.  Therefore,  first  the 
angle  of  incidence  of  our  excitation  beam  had  to  be  measured  to  deduce  the  relative  contribution  of  each 
mode;  it  was  assumed  that  the  three  modes  listed  above  would  have  equal  intensity  contributions  to  the 
total  Frohlich  mode  before  weighting  for  orientation.  For  our  geometry  the  longitudinal  Frohlich  and 
extraordinary  LO  modes  were  given  weights  of  42.9%  each,  and  the  transverse  Frohlich  mode  was  given 
a  weight  of  14.2%. 

The  dielectric  constant  of  the  semiconducting  skeleton  was  calculated  for  700  cm'1  <  co  <  1000 
cm'1  according  to 


£(6))  =  £y 


(  2  2  ^ 
1  m LO  ~  ® TO 


V 


2  2 

COjq  —  co  j 


(3.3) 


£»  was  assigned  a  value  of  8.422  for  our  excitation  wavelength;56  coi0  and  co ro  are  the  experimentally 
determined  longitudinal  optic  mode  and  transverse  optic  mode,  respectively.  Eq.  3.3  assumes  vanishingly 
small  phonon  damping.  The  plot  of  Eq.  3.3  for  anodically  etched  sample  #5  is  shown  in  Fig.  3.16.  For 
this  sample,  the  TO  phonon  at  788.9  cm"1  and  the  LO  phonon  at  964.9  cm'1  can  be  seen,  represented  by  a 
pole  (TO  phonon)  and  a  zero  (LO  phonon). 
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Next,  etT,  the  dielectric  function  component  parallel  to  the  pores,  was  calculated  according  to 

£ti(co)  =  c£(co)  +  (l-c)£0.  (3.4) 

c  is  the  fraction  of  crystalline  skeleton  remaining  after  etching.  The  use  of  Eq.  3.4  to  determine  an 
appropriate  value  for  c  will  be  discussed  below.  s0  is  the  dielectric  constant  of  air  (1.00),  since  air  fills  the 
pores  during  our  measurements.  Eq.  3.4  yields  a  plot  given  in  Fig.  3.17.  In  plots  of  s(eo),  the  poles 
correspond  to  modes  of  transverse  character,  and  the  zeros  correspond  to  longitudinal  modes.  Therefore, 
in  Fig.  3.17,  one  sees  a  TO  mode  at  788.9  cm'1  and  a  zero  at  948  cm'1  which  for  this  plot  corresponds  to 
extraordinary  LO  phonons. 

Finally,  £_l,  the  electric  field  component  perpendicular  to  the  pore  direction,  was  calculated  via 


sL(a>)  =  s(a>) 


(2-c)s0  +  cs(  co) 
csQ  +  (2-c)s(  <x>) 


Eq.  3.5  yields  two  poles  and  two  zeros,  as  shown  in  Fig.  3.18.  The  pole  at  788.9  cm'1  corresponds  to  the 
TO  phonon  mode;  the  zero  at  923  cm'1  corresponds  to  the  longitudinal  Frohlich  excitation;  the  pole  at  958 
cm"1  is  the  transverse  Frohlich  excitation;  and  finally,  the  zero  at  964.9  cm'1  is  the  FO  phonon  mode. 

The  experimental  Raman  spectrum  was  fit  with  two  Gaussian  peaks,  one  Gaussian  representing 
the  crystalline  SiC  and  the  other  Gaussian  peak  the  cumulative  Frohlich  mode  contribution.  The 
wavenumber  for  the  Gaussian  peak  representing  the  cumulative  Frohlich  mode  contribution  was  noted. 
Then,  the  value  for  c  in  Eq.  3.4  and  Eq.  3.5  was  chosen  such  that  the  weighted  contributions  from  the 
extraordinary  FO  phonon  mode  peak  (Eq.  3.4),  the  longitudinal  Frohlich  excitation  peak  (Eq.  3.5),  and 
the  transverse  Frohlich  excitation  peak  (Eq.  3.5)  added  together  matched  the  Gaussian  peak  value  found 
to  represent  the  cumulative  Frohlich  mode  contribution. 

An  appropriate  standard  deviation  was  determined  for  the  three  individual  Frohlich  mode 
Gaussian  peaks.  In  Fig.  3.19,  a  was  assigned  a  value  of  50  to  properly  reflect  the  broad  Frohlich  mode  of 


the  experimental  data.  The  three  individual  Gaussian  peak  contributions  were  added  together  and 
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Figure  3.16.  Plot  of  Eq.  3.3  for  sample  #5. 


Figure  3.17.  Plot  of  Eq.  3.4  for  sample  #5. 
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Figure  3.18.  Plot  of  Eq.  3.5  for  sample  #5. 
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Figure  3.19.  Experimentally-determined  LO  mode  for  a  sample  of  anodically-etched  PSiC  (solid) 
compared  to  the  theoretical  curve  composed  of  Frohlich  mode  interactions  combined  with  a  Gaussian  for 
crystalline  6FI-SiC  (dashed).  Data  is  for  sample  #5. 
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combined  with  a  Gaussian  peak  for  the  crystalline  SiC  to  form  a  calculated  Raman  spectrum  for  our  PSiC 
samples  (Fig.  3.19).  Frohlich  mode  fitting  parameters  are  shown  in  Table  3.1. 

In  general,  good  fits  were  obtained  from  the  bimodal  Gaussian  distribution,  and  Frohlich  modes 
appear  to  be  the  leading  cause  of  the  LO  mode  broadening  witnessed  in  this  work.  However,  for  the  best 
Gaussian  fits,  standard  deviations  of  50  or  greater  had  to  be  used  for  the  Frohlich  modes. 

Table  3.1.  Frohlich  Mode  Fitting  Parameters 
Sample  #  Frohlich  Mode  Peak  (cm~)  Width  (  cm-) 


1 

949.9 

71 

2 

934.9 

56 

3 

955.4 

79 

4 

939.4 

45 

5 

941.9 

57 

6 

938.7 

51 

3.5  Conclusions 

PSiC  produced  via  the  electroless  etching  technique  does  not  show  significant  changes  in  the 
Raman  spectrum  or  cathodoluminescent  emission.  However,  spectroscopic  evidence  of  the 
morphological  changes  in  the  near-surface  region  is  observed  in  the  blue-shift  of  sub-bandgap 
photoluminescent  emission.  Taken  as  a  whole,  these  observations  do  not  bode  well  as  far  as  producing 
PSiC  LEDs. 

The  spectroscopic  behavior  of  the  anodically  etched  PSiC  differed  greatly  from  the  electroless 
material.  Correlation  between  the  Raman,  PL,  and  CL  behavior  is  evident  in  the  anodically  etched 
material.  The  PSiC  wafers  that  yield  the  greatest  broadening  of  the  Raman  LO  mode  also  yield  the 
largest  blue  shifts  in  the  peak  PL  emission.  This  is  similar  to  what  has  been  reported  for  porous  GaP.69  In 
addition,  these  pieces  emit  the  most  intense  CL  signal,  but  the  peak  CL  in  these  pieces  is  relatively  red 
compared  to  the  other  porous  samples.  The  increase  in  UV  CL  emission  appears  to  be  accompanied  by  a 
loss  of  emission  in  the  peak  at  590  nm  (donor-acceptor  recombination),  indicating  a  possible  shift  by 
electrons  from  one  form  of  excited  energy  state  to  another. 


70 


A  few  results  in  this  work  are  consistent  with  a  quantum  confinement  explanation  for  the  trends 
in  CL/PL/Raman  behavior.  First,  CL  peaks  that  are  bluer  than  bandgap  indicate  the  possibility  of 
quantum  wells.  Also,  the  variability  in  the  wavelength  of  maximum  CL  emission  points  to  a  tunability 
that  could  be  dependent  upon  nanocrystallite  size.  In  addition,  the  broadest  Raman  LO  modes  (possibly 
an  indication  of  phonon  quantum  confinement)  correspond  to  the  samples  that  are  dominated  by  UV  CL 
emission,  with  a  loss  of  bandgap  and  donor-acceptor  recombination  peaks  in  the  CL  spectra. 

However,  other  results  in  this  work  argue  against  quantum  confinement.  First,  in  all  samples,  the 
Raman  LO  mode  is  still  dominated  by  a  peak  that  is  unmoved  from  its  position  in  crystalline  SiC.  This 
indicates  that  the  great  majority  of  the  PSiC  layer  is  composed  of  crystallites  too  large  to  exhibit  a 
quantum  confinement  effect.  Also,  fitting  of  the  Raman  spectra  indicate  a  Frohlich  mode  effect  that 
cannot  be  represented  well  by  mathematical  descriptions  of  phonon  quantum  confinement.  Next,  the 
Raman  spectra  fits  yielding  the  smallest  nanocrystallite  sizes  actually  correspond  to  the  pieces  with  the 
smallest  blue  shifts  in  CL  emission.  In  addition,  PL  spectra  never  indicate  a  UV  peak  as  was  seen  in  CL 
measurements.  Finally,  SEM  images  do  not  seem  to  indicate  the  formation  of  numerous  nanocrystallites 
that  are  actually  small  enough  to  produce  quantum  confinement  effects.  Taken  as  a  whole,  the 
spectroscopic  results  (especially  the  fitting  of  Raman  broadening  with  Frohlich  mode  models)  indicate 
that  features  of  the  anodically  etched  porous  surface  (such  as  surface  effects/species  or  defects)  are  most 
likely  responsible  for  the  UV  cathodoluminescent  emission  witnessed  in  this  work. 

Unlike  PSi  the  Raman,  PL,  and  CL  spectra  of  PSiC  are  stable  indefinitely  upon  storage  in  the 
laboratory  ambient,  suggesting  electrolessly  etched  PSiC  is  more  robust  than  PSi  and  therefore  more 
suitable  for  applications  in  hostile  environments.  For  this  reason,  the  chemical  and  temperature  sensing 
capabilities  of  PSiC  were  examined  in  the  remainder  of  this  work. 
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CHAPTER  4 


CHEMICAL  AND  TEMPERATURE  SENSING 


4.1  Introduction 

One  primary  method  of  monitoring  an  engine’s  performance  is  by  measuring  the  percentage  of  H2 
in  the  exhaust  flow.  Therefore,  in  recent  years,  NASA  Glenn  Research  Center  has  led  an  effort  to 
develop  sensors  that  can  detect  H2  concentrations  at  high  temperatures.  However,  the  demands  upon  the 
sensors  are  far  greater  than  current  capabilities.  The  sensor  must  be  able  to  operate  within  a  wide  range  of 
temperatures.  It  must  have  a  low  detection  limit  (<  100  ppm  H2),  yet  exhibit  a  wide  dynamic  range  so 
that  the  sensor  can  also  operate  in  environments  approaching  100%  H2.  The  sensor  must  have  both  a 
quick  response  time  as  well  as  a  rapid  return  to  baseline  signal  upon  removing  H2  from  the  gas  flow.  The 
sensor  must  be  functional  in  both  inert  and  oxygen-rich  environments.  Finally,  the  sensor  must  show 
specificity  and  be  able  to  quantify  the  H2  concentration  amongst  an  environment  of  other  (possibly 
caustic)  gases. 

The  large  specific  surface  area  coupled  with  deep  porosity  and  the  need  for  efficient  transport 
mechanisms  to  shuttle  reactants  and  products  between  the  solid- liquid  interface  and  the  locus  of  etching 
bode  well  for  the  use  of  PSiC  in  high  sensitivity  gas-sensing  applications.  The  large  surface  area 
correlates  directly  with  sensitivity,  and  the  observation  of  a  highly  interconnected  network  of  pores  in  the 
bulk  argues  for  efficient  molecular  transport  to  the  interior  of  the  crystal.  In  this  research,  control  of  the 
etch  depth  was  never  established  using  electroless  etching,  so  PSiC  films  of  with  predictable  porous  layer 
depths  were  never  achieved  via  electroless  etching.  Therefore,  for  the  sensing  studies  that  follow, 
anodically  etched  PSiC  samples  were  used  in  all  experiments,  because  porous  layers  of  desired  thickness, 
as  well  as  a  more  planar  porous-nonporous  interface,  could  be  produced. 

There  are  several  methods  of  actually  measuring  changes  in  the  environment  (  ambient  gas 
identities,  gas  concentrations,  temperature,  pressure)  with  porous  materials.  Ellipsometry  has  been  used 
to  detect  changes  in  the  average  refractive  index  of  porous  layers  due  to  analyte  adsorption.1 4  A 
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reduction  in  photoluminescence  intensity  has  also  been  utilized  to  detect  molecular  constituents.5"10  In 
addition,  interferometry  has  been  used  to  detect  different  chemicals  and  biological  compounds  with  PSi 
sensors.11'18  Finally,  many  studies  have  been  conducted  measuring  changes  in  conductivity  to  sense 
chemicals  in  porous  silicon,  "  crystalline  silicon  carbide,"  ’  and  porous  silicon  carbide."  ’ 

Interferometry  and  conductivity  measurements  were  chosen  as  the  modes  of  detection  in  this  research, 
since  work  with  PSi  reported  in  the  literature  has  shown  these  detection  methods  to  exhibit  promising 
sensitivity,  fast  recovery  to  baseline,  high  dynamic  ranges,  small  minimum  detection  limits,  reversibility, 
fast  time  response,  and  limited  selectivity. 

Experiments  were  designed  to  take  advantage  of  the  inherent  capabilities  and  advantages  of  SiC 
described  in  Chapter  1.  Many  applications  exist  for  a  sensor  that  can  detect  certain  chemicals  in  a  hostile 
(high  temperature,  caustic)  environment.  Applications  that  could  benefit  from  improved  sensing  ability 
exist  in  the  automotive  industry,  wafer  processing,  and  in  monitoring  materials  aging  in  nuclear  reactors. 
The  experiments  detailed  below  were  chosen  to  mimic  conditions  that  may  be  encountered  within  a  liquid 
fuel  rocket  engine,  although  simulating  a  rocket  nozzle  environment  in  a  laboratory  is  obviously  very 
difficult.  Neither  the  extreme  temperatures  nor  unstable  exhaust  species  (such  as  atomic  hydrogen, 
methyl  radicals,  etc.)  within  a  rocket  nozzle  can  be  recreated  easily. 

In  a  rocket  engine,  both  the  concentrations  of  the  exhaust  gases  as  well  as  the  temperature  of  the 
exhaust  flow  can  be  expected  to  fluctuate.  Any  one  method  of  detection  in  itself  will  probably  be 
insufficient  to  deconvolute  the  confounded  effects  of  concentration  and  temperature  changes.  In  addition, 
due  to  selectivity  and  cross-sensitivity,  individual  sensors  are  limited  in  the  amount  of  information  that 
they  can  provide  in  environments  that  contain  multiple  chemical  species.  Therefore,  in  this  research,  both 
interferometry  and  conductivity  measurements  were  attempted  with  the  goal  of  simultaneously  detecting 
rocket  exhaust  species  with  molecular  specificity  as  well  as  the  chamber  temperature.  In  the  future,  for 
any  operational  detection  system,  it  is  envisioned  that  a  large  array  of  different  porous  sensors  would  have 
to  be  constructed  and  probed  in  more  than  one  way  in  order  to  simultaneously  identify  the  various  gas 
species  present,  their  concentrations,  and  the  temperature  of  the  gas  flow. 
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4.2  Interferometry 
4.2.1  Background 

Interferometry  relies  upon  constructive  and  destructive  interference  of  reflections  from  two 
surfaces.  If  a  PSiC  layer  can  be  formed  that  has  a  small  feature  size  compared  to  the  detection  optical 
wavelength,  thin  film  interference  will  occur  from  light  reflected  at  the  air/PSiC  and  PSiC/SiC  interfaces. 
The  resulting  spectrum  will  show  Fabry-Perot  interference  fringes.  The  maxima  in  the  fringes  occur  at 
wavelengths  satisfying  the  relation: 

m-  A  =  2-n-l .  (4.1) 

In  Eq.  4.1,  m  is  the  spectral  order  (an  integer),  X  is  the  wavelength  of  the  fringe  maximum,  n  is  the 
average  refractive  index  of  the  film,  and  /  is  the  thickness  of  the  porous  film.  When  a  gas  adsorbs  to  the 
walls  of  the  porous  film,  the  average  refractive  index  of  the  PSiC  layer  changes,  resulting  in  a  shift  of  the 
Fabry-Perot  fringes  (cf  Fig.  1.1).  The  change  in  n  can  be  attributed  to  two  different  effects.  First,  n  is 
expected  to  increase  due  to  a  higher  refractive  index  species  displacing  air  from  the  nanopores. 
Simultaneously,  the  adsorption/absorption  of  the  analyte  may  change  the  distribution  of  carriers  within 
the  semiconducting  skeleton,  thereby  altering  n  for  the  porous  layer.  Utilizing  Eq.  4.1,  the  shift  in  the 
fringe  maximum  can  be  used  to  quantify  the  amount  of  analyte  in  the  gas  stream  flowing  past  the  porous 
layer,  as  long  as  some  of  the  gas  adsorbs  to  the  walls  of  the  nanopores  in  the  porous  semiconductor. 

The  precise  refractive  index  of  any  given  PSiC  layer  is  a  function  of  the  detailed  structure  and 
porosity  and  so  varies  from  sample-to-sample,  even  for  samples  prepared  under  nominally  identical 
conditions.  Initially,  an  approximate  value  of  the  refractive  index  is  needed,  though,  in  order  to  form 
appropriate  porous  layer  thicknesses  for  the  interferometry  studies.  Petrova-Koch  et.  al.  have  reported  n 
=  2.2±0. 1  for  their  PSiC  samples.10  From  ellipsometry  measurements,  Choyke  has  found  that  the  best  fit 
for  n  is  simply  linear  between  0%  and  100%  porosity  at  632.8  nm  (n= 2.63  to  «=1).31  These  values  were 
used  to  estimate  an  optimum  initial  PSiC  layer  thickness  of  approximately  3  pm  for  interferometry 
studies. 
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Because  the  exact  thickness  of  an  etched  porous  layer  is  not  known  immediately  after  etching 
(until  the  sample  is  cleaved  and  examined  via  SEM),  the  change  in  the  porous  films  due  to  gas  adsoiption 
is  usually  quantified  in  terms  of  nl,  the  effective  optical  thickness  (EOT).  At  the  conclusion  of  the 
interferometry  and  conductivity  studies,  porous  layers  can  be  cleaved  to  measure  the  exact  porous  layer 
thickness  via  SEM.  Once  /  is  known,  a  value  for  the  average  refractive  index  of  the  film  in  nitrogen  and 
the  different  adsorbate  gases  can  be  deduced. 

In  the  literature,  interferometry  of  PSi  layers  has  been  utilized  to  detect  species  such  as 
hydrocarbons,  alcohols,  HF,  CO,  H20,  02,  and  biological  compounds.  In  this  research,  adsorbates  were 
chosen  that  are  anticipated  to  be  found  in  the  exhaust  of  liquid  fuel  rocket  engines.  Unfortunately,  many 
of  the  gases  in  rocket  engines  have  indices  of  refraction  very  close  to  air/nitrogen  (Table  4.1).  As  a  result, 
for  these  species,  changes  in  n  for  the  porous  layer  must  be  due  solely  to  changes  in  the  carrier 
concentration  distribution  in  the  semiconductor  skeleton,  since  the  refractive  index  of  the  pores  will  not 
change  significantly  when  the  adsorbate  has  nearly  the  same  refractive  index  as  air/nitrogen. 

Furthermore,  a  large  dipole  moment  is  required  in  order  to  perturb  the  charge  distribution  in  the 
semiconducting  skeleton  and  be  detected. 

Table  4. 1  Index  of  Refraction  for  Several  Gases  of  Interest 


Species 

SiC 

Air 

Nitrogen 

H2 

CO 

co2 

Ethanol 

Kerosene 

Acetone 

Trimethylamine 


Index  of  Refraction  Dipole  Moment  (D) 

2.6  at  632.8  nm; 

2.7  at  467  nm 


1.0002926 

1.000300 

0 

1.000132 

0 

1.000340 

0.112 

1.00045 

0 

1.36 

1.69 

1.442 

~0 

1.357 

2.88 

1.400 

0.612 
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4.2.2  Experimental 

The  PSiC  samples  were  placed  in  a  custom  made  flow  cell  that  allowed  for  flow  of  nitrogen 
diluent  and  various  adsorbate  gases  through  the  cell  at  controlled  flow  rates  and  concentrations. 

Individual  mass  flow  controllers  controlled  the  nitrogen  and  adsorbate  gases.  Adsorbate  vapors  were 
generated  by  bubbling  a  stream  of  nitrogen  through  analytes  existing  as  liquids  at  300  K,  while  pure  CO 
was  obtained  from  a  compressed  gas  cylinder.  The  total  flow  rate  was  varied  between  20- 1 000  seem. 
Simultaneously  the  PSiC  sample  could  be  heated  from  25  °C  to  600  °C. 

Reflection  spectra  were  obtained  using  the  experimental  arrangement  shown  in  Fig.  4.1.  The 
beam  from  a  tungsten  lamp  was  collimated  and  then  focused  to  a  small  spot  (500  pm  diameter)  on  the 
PSiC  surface.  The  reflected  light  entered  a  spectrometer  (Jobin  Yvon  CP200,  entrance  slit  of  50  pm). 

The  light  was  detected  by  a  liquid  nitrogen  cooled  CCD  camera  (Photometries  CH210). 

It  was  soon  discovered  that  only  a  small  range  of  porous  layer  thicknesses  would  serve  these 
experiments.  If  the  porous  layer  became  too  thin,  then  the  interferometric  fringes  were  spaced  too  far 
apart,  allowing  only  one  or  two  fringes  to  be  monitored  within  the  250  nm  window  of  the 
monochromator.  However,  if  the  porous  layers  were  made  too  thick,  scattering  losses  limited  the 
intensity  of  light  reaching  the  porous-nonporous  interface  that  was  back-reflected  out  to  the  detector.  So, 
for  the  thick  porous  layers,  the  collected  reflections  simply  mimicked  the  source  spectrum,  and  no  fringes 
were  observed.  Optimum  porous  layers  were  found  to  have  thicknesses  of  3-4  pm.  In  addition,  relatively 
low  current  densities  during  etching  yielded  the  best  success  in  producing  Fabry-Perot  fringes,  possibly 
because  the  pore  size  was  small  enough  to  retain  optical  reflectivity  of  the  PSiC  surface. 

CO  is  a  common  rocket  fuel  exhaust  gas,  and  therefore  it  was  examined  first,  even  though  success 
with  CO  was  not  anticipated  due  to  its  low  refractive  index  and  dipole  moment.  Other  rocket  fuel  exhaust 
gases  with  similarly  low  refractive  indices  (such  as  hydrogen,  oxygen,  carbon  dioxide,  and  small 
hydrocarbons)  were  not  investigated.  In  addition  to  CO,  two  other  common  rocket  fuels  (kerosene  and 
ethanol,  a  representative  alcohol)  were  initially  examined  via  interferometry.  Finally,  due  to  the  lack  of 
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success  with  these  vapors  (denoted  below),  acetone  and  trimethylamine  were  bubbled  through  the  test 
cell,  too. 

4.2.3  Results  and  Discussion 

Several  samples  of  PSiC  were  produced  via  anodic  etching  that  yielded  well-resolved  Fabry -Perot 
fringes  in  their  reflection  spectra.  The  reflection  spectrum  for  the  PSiC  layer  in  pure  nitrogen  was 
compared  to  the  spectra  of  the  layer  in  different  analytes,  as  shown  for  example  in  Fig.  4.2. 

Unfortunately,  none  of  the  chemicals  used  in  this  work  could  ever  be  detected  with  PSiC  layers 
via  interferometry.  PSiC  layers  with  a  native  oxide  layer,  as  well  as  samples  that  had  been  F1F  dipped  or 
RCA  cleaned  (Table  4.2),  were  examined  in  these  studies.  CO  (0-17%),  ethanol  (0-7%),  acetone  (0- 
17%),  kerosene  (0-0.25%),  and  trimethylamine  (0-40%)  were  all  investigated,  but  no  shifts  of  the  Fabry- 
Perot  peaks  were  ever  witnessed  at  room  temperature  or  temperatures  up  to  500  °C. 

Table  4.2  RCA  Cleaning  Protocol 

1.  Dip  in  1 : 1 :5  (vol)  NH4OH  :  30%  H202 :  H20  at  70  °C  for  5  min 

2.  Rinse  in  DI  FFO 

3.  Dip  in  5  wt%  F1F  for  5  min 

4.  Rinse  in  DI  H20 

5.  Dip  in  1:1:6  (vol)  37  wt%  HC1 :  H202 :  H20  at  70  °C  for  5  min 

6.  Rinse  in  DI  H20 

4.2.4  Conclusions 

Interferometry  does  not  appear  to  be  a  viable  detection  method  for  the  analytes  of  interest  with 
PSiC  layers.  Even  for  the  relatively  wide  range  of  analyte  properties  examined  in  this  work,  no  shifts  of 
the  Fabry-Perot  fringes  were  ever  witnessed. 
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Figure  4.1.  Experimental  set-up  for  interferometry  collection.  The  “analyte”  was  either  ethanol,  acetone,  kerosene,  or  trimethylamine. 
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Figure  4.2.  Fabry-Perot  fringes  from  a  PSiC  layer  exposed  to  0-17%  acetone.  Notice  the  lack  of  any 
shifting  of  the  Fabry-Perot  fringes,  as  was  observed  for  all  of  the  analyte  species  investigated  in  this  work. 


4.3  Conductivity  Measurements 
4.3.1  Background 

When  a  gaseous  species  absorbs  or  adsorbs  to  a  porous  semiconductor  surface,  the  conductivity 
of  the  porous  region  often  changes  (usually  an  increase  in  conductivity  is  measured).  There  are  several 
possible  explanations  for  the  increased  conductivity,  and  the  explanations  depend  upon  the  identity  of  the 
adsorbate  involved.  In  some  cases,  adsorbates  can  redistribute  charge  in  the  nanocrystallites,  reducing 
charge  carrier  depletion  at  the  free  surface.  In  addition,  at  high  enough  temperatures,  the  conductivity 
may  increase  due  to  absoiption  of  the  adsorbate  into  the  crystal  lattice,  altering  the  lattice  composition 
and,  thus,  the  conductivity.  Adsorbed  species  such  as  hydrogen32  can  also  passivate  dangling  bonds 
(which  can  act  as  charge  traps)  on  the  SiC  surface.  Next,  for  certain  vapors  (at  high  enough 
concentrations),  the  adsorbate  can  condense  into  a  liquid  in  the  pores  and  carry  current  in  addition  to  the 
current  within  the  SiC  matrix.  If  the  adsorbate  does  not  condense,  it  can  at  least  lower  barriers  for  carrier 
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hopping  between  nanocrystallites  that  are  in  close  proximity.  Finally,  the  relationship  between 
conductivity  and  applied  voltage  is  a  function  of  the  adsorbate  concentration  as  well  as  its  dipole  moment, 
so  these  dependencies,  which  are  molecularly  specific,  work  together  to  determine  the  output  signal. 

There  are  several  reports  in  the  literature  detailing  hydrogen  sensing  via  SiC  devices.  All  of  these 
devices  exploit  the  Fl-induced  change  in  barrier  height  at  the  metal-SiC  interface  which  in  turn  shifts  the 
I/V  characteristic,  thereby  changing  the  measured  conductivity.  Published  results  denote  a  wide  range  of 
concentrations  that  can  be  sensed,  reaching  from  the  single  digit  ppm  to  tens  of  percent  of  hydrogen  gas, 
although  no  single  structure  covers  the  entire  composition  range.33'35  In  general,  these  devices  exhibit 
rapid  response  to  the  presence  of  hydrogen  in  the  gas  flow,  usually  responding  on  the  order  of  a  few 
seconds. 

However,  a  variety  of  technological  hurdles  remain  in  order  to  develop  a  viable  hydrogen  sensor 
using  SiC.  Most  of  the  devices  in  the  literature  work  in  a  fairly  tight  range  of  temperatures.  Lower  limits 
on  operating  temperature  are  dictated  by  the  need  for  hydrogen  to  dissociate  on  the  activated  metal,  be 
absorbed,  and  diffuse  through  the  metal  gate.  Upper  temperature  limits  are  often  dictated  by  the 
formation  of  oxides  or  other  deleterious  effects  to  the  activated  metal  gate  material.  Few  devices  are 
reported  to  work  at  room  temperature;35  some  reports  have  shown  response  at  100  °C.36  Upper  limits  of 
operability  are  typically  in  the  500-600  °C  range.  Obviously,  a  wide  range  of  operating  temperatures  is 
desirable. 

The  crystallographic  quality  of  the  metal  gate  material  (usually  Pt  or  Pd)  is  also  an  important 
parameter  governing  the  response  of  hydrogen  sensing  devices.37"44  Both  sensitivity  and  selectivity  are 
highly  dependent  upon  the  morphology  of  the  gate  metal.  Hydrogen  and  some  hydrogen-containing 
compounds  (small  hydrocarbons)  can  usually  be  sensed  regardless  of  the  thickness  and  morphology  of  the 
gate  metal.  However,  whether  other  species,  such  as  CO,  NO2,  and  O2,  contribute  to  the  sensing  signal  is 
very  dependent  upon  the  gate  morphology  and  thickness.  These  chemicals  require  near-direct  contact 
with  the  SiC  surface  in  order  to  be  detected,  so  the  degree  of  coverage  and  “cracks”  in  the  gate  metal 
determine  the  sensitivity  to  these  non-hydrogen  containing  compounds.  Therefore,  exactly  which 
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chemicals  are  detected,  as  well  as  their  individual  sensitivities  by  the  device,  are  all  affected  by  the 
thickness  and  quality  of  the  metal  gate  layer.45,46  Furthermore,  Pd  gates  may  be  restricted  to  relatively 
low  temperatures  (<  400  °C)  due  to  the  formation  of  palladium  oxides  at  higher  temperatures.34  In 
addition,  prolonged  heating  at  425  °C  of  Pd/SiC  Schottky  diode  structures  have  been  shown  to  change  the 
sensor  properties  and  decrease  sensor  response,  partially  due  to  the  formation  of  palladium  silicides  PdxSi 
(x=l,  2,  3,  4)  at  the  interfacial  region.39,40  Finally,  prolonged  exposure  to  elevated  temperatures  can  lead 
to  the  formation  of  PdO  and  SiOx  compounds  on  the  top  surface  of  the  gate  metal.  Even  a  thin  layer  of 
these  compounds  appears  to  have  significant  impact  on  the  catalytic  properties  of  the  surface,  forming  a 
surface  passivation  layer  that  reduces  the  sensitivity  to  gas  adsorption.39 

Just  as  the  processing  conditions  of  the  gate  metal  can  determine  sensor  properties,  selectivity  and 
sensitivity  can  similarly  be  dictated  by  how  the  gate  oxide  is  processed  in  MISiC  devices.  Differences  in 
the  deposition  of  Sn02  gate  films  have  been  shown  to  affect  the  sensor  response,  even  if  the  geometry  of 
the  film  is  the  same.41  Particularly,  tight  control  of  the  thickness  and  stoichiometry  of  the  Sn02  film  are 
required.  Furthermore,  drift  in  the  properties  of  Sn02  films  with  long-term  heating  due  to  grain  boundary 
annealing  has  been  noted.47  The  performance  of  capacitors  utilizing  Si02  insulating  layers  have  shown 
similar  dependence  on  processing  conditions,  especially  annealing.43 

Many  of  the  devices  in  the  literature  that  can  detect  small  (ppm)  concentrations  of  hydrogen 
saturate  at  relatively  low  concentrations  of  hydrogen,  too.  Kim  et.  al.36,48  report  a  SiC  Schottky  diode  that 
can  detect  less  than  160  ppm  H2,  yet  the  device  saturates  at  approximately  600  ppm  H2  (Fig.  4.3).  Similar 
saturation  levels  have  been  found  with  alternative,  thicker  gate  insulators.44  An  attractive  possibility  is 
that  porous  devices  may  exhibit  a  much  larger  range  of  detectable  hydrogen  concentrations. 

The  metal  gates  used  in  the  Schottky  diode  devices  behave  quite  differently  depending  on 
whether  the  gas  mixture  is  oxygen-rich  or  oxygen-lean.33,37,46,49  Therefore,  these  devices  often  yield 
different  results  for  the  same  concentration  of  hydrogen,  depending  on  whether  the  amount  of  oxygen  in 
the  total  flow  is  greater  or  less  than  the  amount  of  hydrogen  being  detected.  Careful  monitoring  of  the  02 
concentration  appears  to  be  a  necessary  component  of  using  Schottky  diode  sensors  to  detect  H2. 
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The  devices  detailed  in  the  literature  have  relatively  slow,  and  sometimes  incomplete,  recovery  to 
baseline  signal  conditions.  For  very  limited  (<  0.1%)  hydrogen  exposure,  Schottky  diode  devices  have 
been  shown  to  recover  to  baseline  conditions  within  a  few  tens  of  seconds.16  After  further  hydrogen 
exposure,  hydrogen  and  oxygen  must  be  balanced  in  order  to  prevent  baseline  drift.  If  the  device  was 
originally  exposed  to  any  oxygen,  then  the  sensors  can  be  “rezeroed”  via  periodic  exposure  to  oxygen 
between  hydrogen  sensing  applications  (Fig.  4.4).45  49  Alternatively,  a  device  can  be  outgassed  before  use 
in  a  pure  nitrogen  environment  in  order  to  desorb  all  oxygen  from  the  device,  then  baseline  drift  is 
avoided  as  long  as  the  device  is  never  exposed  to  oxygen.  Initial  outgassing  of  the  device  to  remove 
oxygen  can  take  hours,  though.45  Even  after  outgassing  all  oxygen  from  the  device,  signal  drift  can  be 
observed  as  a  result  of  incomplete  hydrogen  desoiption  from  the  device  (Fig.  4.5).5U  In  some  cases  (when 
the  exposure  to  high  temperature  is  long  enough),  the  drift  in  the  baseline  is  unrecoverable.  Diodes  made 
of  Pd  on  SiC  have  shown  background  drift  due  in  part  to  reactions  between  the  Pd  and  the  SiC;40’42  some 
method  of  stabilizing  the  Pd/SiC  interface  must  be  found  in  order  to  avoid  this  type  of  baseline  signal 
drift.  Like  the  Schottky  devices,  MISiC  devices  have  shown  similar  problems  of  baseline  drift.14  These 
problems  currently  prohibit  NASA  from  utilizing  these  sensors  in  the  numerous  uses  in  which  the  sensor 
must  be  able  to  survive  exposure  to  100%  hydrogen  without  damage  or  change  in  calibration.51 

Finally,  devices  that  detect  hydrogen  often  have  similar  responses  to  other  gases.  Methane,36 
propane,37’45'49  butane,45  ethylene,42’45  propylene,42’49  CF4,35  NH3,46  CO, 45,46,49  and  NCE45  have  all  been 
detected  by  solid  state  SiC  devices.  Obviously,  in  a  mixture  of  gases,  all  of  the  individual  contributions 
of  these  components  would  have  to  be  isolated  and  measured  separately.  Therefore,  the  additive  response 
of  single  devices  does  not  bode  well  for  applications  in  which  a  single  gas  of  interest  in  a  mixture  needs  to 
be  monitored.  Rather,  fabrication  of  an  array  of  uniquely  responding  sensors  is  the  most  likely  solution. 

In  order  to  measure  changes  in  conductivity  within  porous  layers,  ohmic  contacts  to  the  porous 
morphology  need  to  be  formed.  Ohmic  contacts  are  much  easier  to  make  when  the  SiC  wafer  is  heavily 
n-doped.  Because  our  initial  wafers  had  relatively  low  doping  levels  (less  than  3  x  1018  cm'3),  a  new 
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Figure  4.3. 36  Example  of  data  from  a  PSiC  Schottky  diode  showing  low  detection  limits  for  H2  coupled 
with  relatively  low  saturation  points. 


Time  [min] 

Figure  4.4.46  Data  showing  drift  in  the  baseline  conditions  for  a  Pt-TaSi-SiC  Schottky  diode  sensor 
operating  at  500  °C.  The  sensor  drift  is  believed  to  be  due  to  slow  desorption  of  surface  oxygen  ions  from 
the  sensor  during  the  exposure  to  the  N2  background  atmosphere. 
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Figure  4. 5. 51  Response  of  a  Pd/SiC  Schottky  sensor  to  H2  at  400  °C.  Notice  the  rising  baseline  condition 
with  repeated  H2  exposure. 

wafer  was  ordered  with  a  higher  n-doping  level  (8  x  1018  cm'3).  The  conductivity  tests  detailed  below 
were  conducted  with  this  highly  n-doped  6H-SiC  wafer. 

The  literature  details  numerous  studies  addressing  which  contact  materials  form  the  most  stable, 
ohmic  contacts  to  SiC.52"65  However,  the  most  simple  and  relatively  reliable  ohmic  contacts  are  formed 
by  depositing  nickel  on  the  SiC  surface.66'70  These  contacts  have  shown  excellent  thermal  stability  up  to 
500  °C.52'66  For  nickel  to  form  an  ohmic  contact  with  the  SiC,  nickel  silicide  (Ni2Si)  must  form. 
Therefore,  an  annealing  step  at  temperatures  of  about  800-900  °C  or  greater  is  necessary.  However, 
because  the  silicon  in  the  SiC  is  consumed  in  the  formation  of  the  silicide,  excess  carbon  accumulates  at 
the  contact/SiC  interface.  This  carbon  then  has  a  tendency  to  migrate  to  the  free  surface  (and  form 
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graphite  nanocrystallites)  and  leave  microvoids  at  the  contact/SiC  interface.71"75  These  microvoids  do  not 
appear  to  be  detrimental  to  contact  performance  at  temperatures  below  500  °C  for  hundreds  of  hours.66 

Another  promising  method  of  producing  ohmic  contacts  to  SiC  is  detailed  in  a  series  of  papers  by 
Pecz  et.  al.71"74  They  deposit  alternating  layers  of  Ni  and  Si  in  the  correct  proportion  to  form  Ni2Si,  and 
therefore  no  silicon  in  the  SiC  wafer  must  be  consumed  to  form  the  ohmic  contact.  Pecz  and  coworkers 
have  found  one  effective  architecture  to  be  Si(15  nm)/Ni(16.4  nm)/Si(30  nm)/Ni  (32.8  nm)/Si(30 
nm)/Ni(32.8  nm)/SiC.  Upon  annealing,  small  specific  contact  resistances  (~1  x  10’5  Q  cm2)  along  with 
excellent  adhesion  and  relatively  low  microvoid  volumes  were  observed. 

Of  course,  most  of  the  work  reporting  ohmic  contacts  to  SiC  has  utilized  only  crystalline  SiC; 
very  little  work  has  been  performed  to  establish  ohmic  contacts  to  a  porous  layer.  Shields  et.  al.28’29  have 
produced  contacts  with  PSiC  that  display  fairly  ohmic  behavior;  Wu  et.  al.76  have  also  made  electrical 
contact  to  n-doped  and  p-doped  PSiC,  although  the  quality  of  the  contact  was  not  evaluated  during  their 
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photodetecting  studies.  Slightly  more  work  has  been  performed  investigating  ohmic  contacts  to  PSi,  ’ 
but  the  nature  of  ohmic  contacts  to  Si  is  much  different  than  contacts  to  SiC,  and  the  electrical  properties 
are  quite  different. 

Specific  contact  resistance  is  a  property  used  to  characterize  a  contact  as  ohmic  or  Schottky. 
“Ohmic”  contacts  are  generally  expected  to  have  specific  contact  resistances  of  ~1  x  10"5  Q  cm2  or  less. 
There  are  several  methods  of  measuring  the  specific  contact  resistance.  KuphaTs  modified  4-point  probe 
technique  uses  structures  which  are  easy  to  fabricate  but  is  relatively  inaccurate,  as  this  method  does  not 
account  for  current  crowding  at  the  contact  edge.80  The  circular  transfer  length  method  (TLM)  requires 
tedious  mathematical  analysis  that  can  be  inaccurate.  The  method  of  Cox  and  Strack81  is  unique  in  that  it 
vertically  measures  the  resistance,  requiring  four  different  sized  circular  contact  pads  on  the  top  surface  of 
the  wafer  and  a  sheet  of  ohmic  contact  on  the  bottom  surface  of  the  wafer.  The  linear  TLM  was  used  in 
this  work,  though,  as  this  method  allows  for  easy  extraction  of  relatively  small  values  of  lateral  contact 
resistance.  This  geometry  also  limits  the  resistance  measurements  to  the  porous  layers. 
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TLM  involves  measuring  the  total  resistance  between  multiple  contacts  with  unequal  spacing 
between  the  contacts.  Analysis  of  the  data  yields  the  sheet  resistance,  the  contact  resistance,  and  the 
specific  contact  resistance.  In  this  work,  TLM  was  first  used  to  ensure  that  ohmic  contacts  could  be 
deposited  on  our  crystalline  SiC  wafers.  These  results  were  then  compared  to  contacts  deposited  on 
porous  layers  to  assess  changes  in  the  contact  properties  and  the  suitability  of  porous  contacts  for  sensing 
studies. 

An  excellent  review  of  TLM  data  interpretation  is  presented  by  Schroder.82  In  general  the  total 
resistance  ( RT )  is  measured  for  various  contact  spacings,  d,  between  the  different  contact  pads.  RT  is  then 
plotted  as  a  function  of  d.  Analysis  of  the  data  reveals  three  pieces  of  data.  First,  the  slope  of  the  plot 
yields  the  sheet  resistance  ( ps ),  as 


A  (Rt)Z 
A  (d) 


(4.2) 


where  Z  is  the  contact  width.  Next,  the  contact  resistance  Rc  can  be  calculated  from  the  y-intercept,  as 


R 


C 


(4.3) 


Finally,  the  specific  contact  resistance  pc  can  be  calculated  from  the  x-intercept  using 


(4.4) 


in  which  Lr  is  the  transfer  length, 

d  =  2Lt.  (4.5) 

As  in  the  crystalline  SiC  hydrogen  sensors  discussed  above,  several  options  exist  for  the  exact 
architecture  of  PSiC  sensors.  However,  two  basic  unique  sensing  configurations  can  be  used  (Fig.  4.6). 
First,  a  capacitor  type  configuration  can  be  constructed  in  which  current  flows  vertically  through  the  top 
porous  (insulating)  layer  and  then  the  crystalline  substrate  underneath.  A  grid  contact  can  be  formed  on 
top  of  the  PSiC  layer,  and  a  solid  contact  on  the  substrate’s  bottom.  The  grid  contact  pattern  is  beneficial 
because  it  produces  a  relatively  large  and  even  contact  over  the  entire  surface  area,  while  allowing  for 
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enough  free  surface  area  for  analyte  adsoiption.  Sensing  is  then  accomplished  by  measuring  the 
magnitude  of  current  flow  as  a  function  of  the  applied  bias. 

Conversely,  a  horizontal  (FET)  configuration  can  be  used.  In  this  configuration,  the  porous  layer 
becomes  the  “gate.”  The  current  through  the  device  is  measured  as  a  function  of  applied  potential 
between  the  “source”  and  “drain”  contacts.  Flowever,  the  effectiveness  of  this  architecture  for  a  bulk  SiC 
wafer  is  suspect,  given  the  crystalline  SiC  substrate  could  short-circuit  the  device  and  provide  a  path  for 
the  current  to  bypass  the  porous  region. 

This  research  first  investigated  the  formation  of  ohmic  contacts  to  crystalline  and  then  porous 
SiC.  After  reliable  contacts  could  be  made  to  PSiC  with  predictable  specific  resistances,  tests  were 
conducted  to  characterize  the  suitability  of  porous  layers  to  measure  different  chemical  environments  via 
conductivity. 


contact 


contact 


(b) 


Figure  4.6.  Two  basic  conductivity  sensing  configurations  to  measure  the  presence  of  adsorbate  gases 
within  a  porous  layer:  (a)  a  capacitive  device,  and  (b)  a  FET  configuration. 


4.3.2  Experimental 

4.3.2. 1  Specific  Contact  Resistance  Measurements 

Initially,  TLM  investigations  were  conducted  to  ensure  ohmic  contact  could  be  made  to  both  the 
crystalline  as  well  as  the  porous  SiC  layers.  Therefore,  wafer  samples  suitable  for  TLM  specific  contact 
resistance  measurements  were  fabricated.  Many  variations  of  sample  preparation  were  studied  (as 
discussed  in  section  4.3.3),  but  the  general  procedure  is  as  follows.  Before  processing,  the  SiC  sample 
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was  dipped  in  HF  and  rinsed  in  water  in  order  to  remove  any  native  oxide.  The  wafer  was  then  ICP-RIE 
etched  in  SF6  (25  °C,  1.5  mtorr  total  pressure,  8.0  seem  SF6,  2.0  seem  02,  -200  V  DC  with  a  4”  diameter 
plate,  800  W  coil  power,  3  min  etch  time)  to  remove  the  top  500  nm  of  the  wafer  and  provide  a  clean  SiC 
surface.  The  sample  was  coated  with  5214E  photoresist  (applied  at  5000  rpm  for  30  sec,  resulting  in  ~1.5 
pm  thick  photoresist)  and  subsequently  baked  at  1 10  °C  for  35  sec.  The  wafer  was  then  patterned  via  UV 
lithography  to  form  patterns  of  mesas  100  pm  wide  by  650  pm  long.  After  UV  exposure  (approximately 
400  watts  for  20  sec),  the  wafer  was  developed  (AZ327MIF  or  AZ300MIF  for  1  min),  rinsed  in  water, 
and  blown  dry.  The  wafer  was  next  ICP-RIE  etched  again  to  form  mesas  in  the  SiC  approximately  1  pm 
tall.  These  mesas  were  needed  to  isolate  the  TLM  pads  for  conductivity  measurements.  Results  in  the 
literature83  have  utilized  mesas  about  2.5  pm  tall  for  TLM  measurements  of  bulk  6H-SiC.  Electrical 
isolation  of  our  mesas  was  ensured  by  measuring  large  resistances  (off  scale,  in  excess  of  1  GO)  between 
pads  on  neighboring  mesas. 

After  formation  of  the  mesas,  the  remaining  photoresist  on  top  of  the  mesas  was  removed  with 
acetone  followed  by  cleaning  in  isopropyl  alcohol  and  water  and  finally  blown  dry.  Next,  the  pattern  for 
the  TLM  pads  had  to  be  applied.  The  wafer  was  coated  again  with  photoresist  as  described  above.  The 
TLM  pad  pattern  was  applied  via  UV  lithography  as  described  above,  except  that  now  the  exposure  time 
was  only  8  sec,  as  an  inverse  pattern  of  the  mask  had  to  be  formed.  In  this  work,  eight  identically  sized 
pads  (100  pm  wide  by  75  pm  long)  were  deposited  in  a  row,  the  pads  separated  by  2,  3,  4,  5,  6,  8,  and  10 
pm  respectively.  The  photoresist  was  cured  at  125  °C  for  1  min.  Next,  the  entire  wafer  was  exposed  to 
ultraviolet  irradiation  (total  dose  of  375  mJ/cm2  at  365  nm,  or  6.8  mW/crn2,  125  sec  exposure  time).  The 
photoresist  was  developed  (AZ327MIF  or  AZ300MIF  for  ~3  min)  and  rinsed  in  water.  Oxygen  plasma 
treatment  (02  at  0.8  W/in2  for  30  sec,  500  mtorr  total  pressure)  ensued  in  order  to  clean  any  residual  resist 
from  the  wafer  surface. 

Several  different  metallization  schemes  were  investigated.  Nickel,  titanium,  and  nichrome  (80% 
nickel  and  20%  chromium  by  weight)  contact  material  were  deposited  via  electron  beam  evaporation  at 
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pressures  below  8  x  10~7  torr.  Ni,  Ti,  and  NiCr  were  applied  at  deposition  rates  of  approximately  0.2 
nm/sec.  Alternatively,  contacts  made  of  alternating  Ni/Si  layers  or  Ni/Al  layers  were  deposited  via  argon 
sputtering  at  pressures  below  3  x  1  O'6  torr.  Gold  was  thermally  evaporated  onto  the  contacts  at  pressures 
below  8  x  10"7  torr  at  a  rate  of  approximately  0.2  nm/sec.  After  metal  deposition,  lift  off  of  excess  metal 
occurred  in  an  acetone  dip  for  1  min. 

The  annealing  time  and  temperature  has  proven  to  be  a  critical  parameter  in  the  Schottky-to- 
ohmic  contact  transition,80  and  therefore  various  annealing  conditions  were  investigated  in  the  TLM 
portion  of  this  work.  Annealing  temperatures  between  800-1100  °C  and  times  from  0-20  min  were 
performed  in  pure  dry  nitrogen. 

4.3.2.2  Fabrication  and  Testing  of  Devices 

All  of  the  sensors  fabricated  for  this  test  series  were  anodically  etched  in  a  similar  manner.  The 
wafers  were  etched  for  60  min  with  a  current  of  4.5  mA.  This  consistently  yielded  a  PSiC  film 
approximately  4  pm  thick  in  each  sensor  (as  verified  by  SEM  measurements  at  the  conclusion  of  sensing 
studies). 

During  lift  off  of  the  excess  metal,  a  short  duration  sonication  was  often  needed  in  order  to  aid  in 
the  removal  of  the  excess  metal.  This  sonication  did  not  damage  the  contacts  made  to  crystalline  SiC. 
However,  it  was  found  that  sonication  would  often  fracture  the  porous  SiC  film.  Therefore,  sonication 
was  not  used  with  porous  wafers,  and  the  success  rate  of  forming  isolated  TLM  pads  on  porous  SiC  was 
low  (~2%).  As  a  result,  when  making  actual  devices,  a  new  (thicker)  photoresist  was  tested  and  used. 

The  SiC  sensor  sample  measuring  5x5  mm  was  first  dipped  in  acetone,  methanol,  HF,  and  then 
water  to  initially  clean  the  wafer.  The  sample  was  coated  with  AZ4330  photoresist  (applied  at  4000  rpm 
for  30  sec,  resulting  in  ~5  pm  thick  photoresist)  and  subsequently  baked  at  100  °C  for  1  min.  The  sample 
was  then  patterned  with  a  mask  4.8  x  4.8  mm  in  size,  containing  24  squares  800  pm  on  a  side  (Fig.  4.7). 
The  mask  was  illuminated  at  a  total  dose  of  400  mJ/cm2  at  365  nm.  Next,  the  pattern  was  developed  for 
approximately  90  sec  using  AZ400K  developer  diluted  1 :4  in  DI  water.  Oxygen  plasma  treatment 
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ensued,  followed  by  metal  deposition.  Finally,  liftoff  in  acetone  was  used  (without  any  sonication)  to 
remove  the  excess  metal. 


When  making  devices,  bonding  wire  to  the  contacts  was  necessary.  Originally,  a  Flybrid  Model 
572  Thermosonic  Gold  Wire  Bonder  with  0.001”  diameter  gold  wire  was  used.  Flowever,  bonding  wire 
applied  with  this  bonder  often  would  not  adhere  well  to  the  device,  and  often  large  craters  formed  in  the 
contacts  (Fig.  4.8).  Therefore,  0.01”  diameter  gold  wire  was  bonded  to  the  sensing  devices  with  silver 
conductive  epoxy. 

The  devices  were  placed  into  a  steel  flow  cell  for  conductivity  tests.  The  porous  device  and  a  K- 
type  thermocouple  were  suspended  1  mm  from  a  strip  of  nichrome  that  was  resistively  heated  to  the 
desired  temperature.  Temperature  was  manipulated  during  the  sensing  studies  in  one  of  two  ways.  First, 

Porous  SiC  Metal  Grid 


Figure  4.7.  Masking  pattern  used  for  sensor  fabrication.  The  black  areas  represent  the  metal  contact  (200 
pm  wide  around  the  periphery  and  100  pm  wide  in  the  interior). 
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Figure  4.8.  SEM  images  of  three  bonding  sites  performed  with  the  Hybrid  Model  572  Thermosonic  Gold 
Wire  Bonder.  Notice  the  craters  formed  in  the  gold  TLM  pads.  Because  these  craters  could  not  easily  be 
eliminated,  and  the  0.001”  diameter  bonding  gold  wire  was  not  robust  to  manual  handing,  manual 
bonding  with  conductive  epoxy  and  0.01”  diameter  gold  wire  was  chosen  to  make  the  sensing  devices. 


in  constant  temperature  tests,  the  analyte  gas  was  adjusted,  and  then  the  temperature  of  the  nichrome 
heating  strip  was  adjusted  to  keep  the  temperature  of  the  porous  device  and  thermocouple  constant. 
Alternatively,  in  constant  applied  heat  tests,  once  an  original  temperature  was  attained,  further 
manipulations  of  the  analyte  gas  flow  were  conducted  with  no  change  in  the  resistance  through  the 
nichrome  heating  strip,  allowing  the  temperature  of  the  porous  device  and  thermocouple  to  rise  during  the 
test  as  the  analyte  concentration  increased.  In  this  mode  of  operation  the  suspended  sensor  experienced 
both  conductive  and  convective  heating  from  the  gas  flow. 

Several  variations  of  data  collection  were  attempted.  The  data  presented  below  were  collected  by 
sweeping  the  applied  voltage  across  a  set  range  and  measuring  the  current  at  each  voltage  step.  Then,  the 
slope  of  the  I/V  plot  was  used  to  determine  the  resistance  of  the  porous  sensor  for  each  gas  flow  and 
temperature  condition.  Several  different  voltage  ranges  were  tested,  including  -2500  to  2500  mV  (50  mV 
steps,  40  sec  collection  time),  -500  to  500  mV  (25  mV  steps,  8  sec  collection  time),  and  -5  to  5  mV  (1  mV 
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steps,  4  sec  collection  time).  Results  between  these  different  voltage  ranges  were  indistinguishable; 
therefore,  for  most  of  the  data  below,  sweeps  from  -5  to  5  mV  were  used  since  this  consumed  the  least 
amount  of  power  (as  would  be  desirable  in  a  functioning  sensor).  Of  course,  data  collection  time  could  be 
reduced  if  one  single  voltage  could  be  applied  and  the  measurement  made  from  that  single  resistance. 
However,  it  was  found  that  single  voltage  measurements  were  considerably  more  noisy  and  inconclusive 
than  the  -5  to  5  mV  voltage  sweep.  In  addition,  AC  measurements  of  the  resistance  were  made  with  an 
AC  lock-in  amplifier.  In  general,  the  data  trends  and  magnitudes  were  identical  to  those  found  via  the  DC 
experiments.  However,  since  once  again  the  AC  measurements  were  single  point  measurements,  the  data 
was  considerably  noisier.  As  a  result,  we  can  conclude  that  a  majority  of  the  scatter  in  the  data  for  a 
single  sensor  is  due  to  thermal  noise,  and  1/f  noise  is  probably  negligible  in  these  DC  voltage  sweeps. 

Initial  studies  were  performed  to  quantify  the  combinatorial  effects  on  conductivity  by  an 
adsorbate  gas  simultaneous  with  temperature  changes.  Adsorbates  tested  in  the  conductometric 
experiments  include  N2,  H2,  ethanol,  and  acetone.  The  H2,  ethanol,  and  acetone  were  each  diluted  in  N2  to 
the  desired  concentration.  The  total  flow  rate  was  100  seem.  The  wafer  temperature  was  varied  between 
25-500  °C. 

4.3.3  Results  and  Discussion 

4.3.3. 1  Contact  Formation 

4.3.3. 1.1  Crystalline  SiC 

Initially,  pure  Ni  contacts  (150  nm  thick)  were  studied  via  TLM.  As  expected,  the  unannealed 
contacts  were  Schottky  (Fig.  4.9).  After  annealing  at  800-1100  °C  for  30  sec- 10  min,  none  of  the  nickel 
contacts  adhered  well  to  the  crystalline  SiC  surface,  as  the  TLM  sampling  probes  easily  caused  the  nickel 
to  crumble  and  delaminate  from  the  surface.  Therefore,  no  current  would  flow  between  neighboring 
annealed  contact  pads. 

Next,  150  nm  thick  Ni  pads  were  capped  with  50  nm  of  gold  before  annealing.  The  gold  cap 
layer  was  applied  to  prevent  any  oxidation  of  the  nickel  contact.  During  TLM  measurements,  the  gold 
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cap  layer  eliminated  the  crumbling  of  the  contact,  but  delamination  of  the  annealed  contacts  remained  an 
issue. 

Adhesion  problems  have  been  successfully  alleviated  in  the  literature  by  depositing  a  thin  Ti  layer 
first  to  the  substrate.  An  added  benefit  of  applying  Ti  to  Ni  contacts  on  SiC  is  that  the  Ti  consumes  some 
of  the  excess  carbon  to  form  TiC  while  maintaining  good  ohmic  contact  behavior. 70,8,  84  Therefore,  the 
next  set  of  contacts  were  made  of  50  A  Ti,  1500  A  Ni,  and  then  a  gold  cap  layer  500  A  thick.  These 
contacts  were  structurally  sound  and  did  not  delaminate  from  the  SiC  wafer.  The  unannealed  contacts 
showed  good  adhesion  to  the  underlying  crystalline  SiC  substrate,  although  some  discontinuities  existed 
in  the  Ni  layer  (Fig.  4.10).  Annealing  at  900  °C  for  1  min  or  more  resulted  in  improved  contact 
homogeneity  (Fig.  4.1 1).  The  unannealed  contacts  were  Schottky,  with  results  very  similar  to  the  pure 
nickel  contact  results  shown  in  Fig.  4.9.  After  testing  various  annealing  conditions  (800-1100  °C  for  0-20 
min),  optimum  conductivity  was  found  to  occur  after  annealing  at  900  °C  for  1 0  min.  Annealing  times  in 
excess  of  10  min  result  in  no  reduction  of  the  contact’s  resistance.  An  example  of  contact  behavior  as  a 
function  of  annealing  time  is  shown  in  Fig.  4. 12.  10  min  annealing  times  at  900  °C  resulted  in  ohmic 
contacts  forming  -15%  of  the  time.  Figs.  4.13  and  4.14  present  data  collected  from  TLM  pads  for  these 
contacts.  This  data  yields  a  specific  contact  resistance  of  1.3  x  10'5  Q  cm2  for  the  crystalline  SiC  contacts. 

Because  the  Ti/Ni/Au  contacts  were  ohmic  only  ~15%  of  the  time,  several  other  metallization 
schemes  were  investigated  as  described  below.  However,  none  of  the  contacts  described  below  exhibited 
a  greater  rate  of  success  than  the  Ti/Ni/Au  contacts,  and  eventually  all  devices  were  made  with  Ti/Ni/Au 
or  Ti/NiCr(20%  Cr)/Au  contacts. 

Following  the  work  of  Pecz  et.  al.71"74  described  in  section  4.3.1,  alternating  Ni  and  Si  layers  were 
deposited  via  argon  sputtering  in  the  stoichiometric  ratio  to  form  Ni2Si.  Annealing  at  1000  °C  for  3-10 
min  yielded  contacts  that  were  ohmic  -100%  of  the  time,  so  these  contacts  were  very  promising. 

However,  SEM  analysis  of  the  contact  morphology  indicated  these  contacts  would  be  problematic  to  cap 
with  gold.  Gold  capping  is  necessary  to  protect  the  contact  from  oxidation  as  well  as  allow  for  bonding. 
Argon  sputter  deposition  resulted  in  poor  adhesion  of  the  unannealed  contact  materials  to  a  rough  surface, 
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Figure  4.9.  I/V  curve  of  an  unannealed  150  nm  thick  nickel  contact  to  crystalline  SiC,  showing  the 
expected  Schottky  character. 


Figure  4.10.  Cross-sectional  SEM  image  of  an  unannealed  contact  composed  of  50  A  Ti,  1500  A  Ni,  and 
500  A  Au.  The  substrate  is  crystalline  SiC. 
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Figure  4.1 1.  Cross-sectional  SEM  image  of  a  contact  composed  of  50  A  Ti,  1500  A  Ni,  and  500  A  Au, 
annealed  at  900  °C  for  1  min.  The  substrate  is  crystalline  SiC. 

as  shown  in  Fig.  4.15.  Annealing  at  900-1000  °C  resulted  in  improved  adhesion  to  the  SiC  substrate; 
however,  the  contact  material  formed  large  voids  during  annealing.  In  Fig.  4.16,  after  just  20  sec  of 
annealing,  voids  can  be  seen  to  form  in  the  Ni  and  Si  layers  reacting  to  form  Ni2Si.  In  other  samples, 
large  sheet-like  voids  were  formed  during  annealing  (Fig.  4.17).  Even  longer  annealing  times  resulted  in 
contacts  composed  of  many  nanowires  and  large  pockets  (Fig.  4.18  and  Fig.  4.19).  When  gold  cap  layers 
were  thermally  evaporated  onto  the  Ni/Si  layers  prior  to  annealing  at  1000  °C,  the  gold  caps  failed  with 
only  short  annealing  times.  Within  a  few  minutes  of  annealing,  the  gold  cap  layers  would  become 
discontinuous  (Fig.  4.20).  When  the  contacts  were  annealed  for  10  min  or  more,  the  change  in  contact 
morphology  caused  nearly  all  of  the  gold  cap  layers  to  form  small  pools  or  completely  delaminate  from 
the  contact.  Alternatively,  a  protocol  was  investigated  in  which  the  Ni/Si  layers  were  annealed  to  make 
ohmic  contacts  and  then  the  gold  cap  layer  was  deposited.  Elowever,  no  better  success  rate  was  achieved. 
A  significant  (25%  or  more)  portion  of  the  gold  cap  layers  would  delaminate  from  the  contacts  during  lift 
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—  1  min 
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—  no  anneal 
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3  min 

—  5  min 
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Figure  4.12.  I/V  curves  for  contacts  composed  of  50  A  Ti,  1500  A  Ni,  and  500  A  Au  annealed  at  900  °C 
for  the  times  indicated.  The  substrate  is  crystalline  SiC.  (a)  is  shown  in  full  scale,  while  (b)  has  the 
current  scale  enlarged  to  show  detail.  Anneal  times  in  excess  of  1 0  min  result  in  no  reduction  of  the 
contact  resistance. 
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(a) 


(b) 


Figure  4.13.  I/V  data  collected  from  a  row  of  contact  pads  composed  of  50  A  Ti,  1500  A  Ni,  and  500  A 
Au  annealed  at  900  °C  for  10  min.  The  substrate  is  crystalline  SiC.  (a)  is  shown  in  full  scale,  while  (b) 
has  the  current  scale  enlarged  to  show  detail. 
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Pad  Spacing,  d  (microns) 


Figure  4.14.  TLM  data  collected  from  a  row  of  contact  pads  composed  of  50  A  Ti,  1500  A  Ni,  and  500  A 
Au  annealed  at  900  °C  for  10  min.  The  substrate  is  crystalline  SiC.  Analysis  of  the  data  according  to 
Eqs.  4. 2-4. 5  yields  a  specific  contact  resistance  of  1.3  x  10'5  Q  cm2  for  these  contacts. 


Figure  4.15.  Cross-sectional  SEM  image  of  alternating  Ni  and  Si  layers  argon  sputtered  onto  a  rough  SiC 
crystalline  surface.  Notice  the  poor  adhesion  to  the  crystalline  SiC  substrate. 
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Figure  4.16.  Cross-sectional  SEM  image  of  a  contact  composed  of  argon  sputtered  Si(  1 5  nm)/Ni(16.4 
nm)/Si(30  nm)/Ni(32.8  nm)/Si(30  nm)/Ni(32.8  nm)/SiC  that  has  been  capped  with  25  nm  Au  and  then 
annealed  for  20  sec  at  900  °C.  Notice  the  voids  forming  during  the  annealing  reaction.  The  substrate  is 
crystalline  SiC. 


Figure  4.17.  Cross-sectional  SEM  image  of  a  contact  composed  of  argon  sputtered  Si(  1 5  nm)/Ni(16.4 
nm)/Si(30  nm)/Ni(32.8  nm)/Si(30  nm)/Ni(32.8  nm)/SiC  that  has  been  annealed  for  3  min  at  1000  °C  and 
subsequently  capped  with  250  nm  Au.  Reaction  between  the  alternating  Ni  and  Si  layers  resulted  in  large 
voids  and  poor  contact  adhesion.  The  substrate  is  crystalline  SiC. 
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Figure  4.18.  Plan  view  SEM  image  of  a  contact  composed  of  argon  sputtered  Si(  1 5  nm)/Ni(16.4 
nm)/Si(30  nm)/Ni(32.8  nm)/Si(30  nm)/Ni(32.8  nm)/SiC  that  has  been  annealed  for  10  min  at  1000  °C. 
Reaction  between  the  alternating  Ni  and  Si  layers  resulted  in  cratering  of  the  top  surface  as  well  as  the 
formation  of  nanowires.  The  substrate  is  crystalline  SiC. 


Figure  4.19.  Cross-sectional  SEM  image  of  a  contact  composed  of  argon  sputtered  Si(  1 5  nm)/Ni(16.4 
nm)/Si(30  nm)/Ni(32.8  nm)/Si(30  nm)/Ni(32.8  nm)/SiC  that  has  been  annealed  for  10  min  at  1000  °C. 
Reaction  between  the  alternating  Ni  and  Si  layers  resulted  in  the  formation  of  nanowires.  The  substrate  is 
crystalline  SiC. 
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Figure  4.20.  Plan  view  SEM  image  of  a  contact  composed  of  argon  sputtered  Si(  1 5  nm)/Ni(16.4 
nm)/Si(30  nm)/Ni(32.8  nm)/Si(30  nm)/Ni(32.8  nm)/SiC  that  has  been  capped  with  250  nm  gold  and  then 
annealed  for  2.5  min  at  1000  °C.  Reaction  between  the  alternating  Ni  and  Si  layers  resulted  in  gaps 
forming  in  the  gold  cap  layer.  The  substrate  is  crystalline  SiC. 

off.  In  addition,  many  of  the  remaining  contacts  were  Schottky  in  nature  (even  though  the  gold  was 
deposited  over  an  ohmic  Ni/Si  contact),  pointing  to  poor  Au  adhesion  to  the  Ni/Si  contact. 

Nichrome  has  also  been  suggested  as  a  possible  contact  material,  particularly  showing  excellent 
adhesion  to  the  gold  cap  layer  and  improved  bonding  characteristics.63  Chromium  also  can  form  a  stable 
carbide  (CrTT),  providing  a  benefit  similar  to  the  presence  of  titanium.70  Nichrome  contacts  investigated 
in  this  work  (1500  A  NiCr,  500  A  Au)  were  ohmic  ~15%  of  the  time  (similar  to  the  Ti/Ni/Au  contacts), 
and  showed  specific  contact  resistances  comparable  to  the  Ti/Ni  contacts.  Contact  conductivity  as  a 
function  of  annealing  time  is  shown  in  Fig.  4.21 ;  results  are  similar  to  those  found  for  the  Ti/Ni  contacts. 

Finally,  aluminum  was  investigated  in  the  contacts  to  prevent  the  formation  of  microvoids  by 
reacting  with  excess  carbon  to  form  AI4C3.60  Contacts  composed  of  200  A  A1  and  1250  A  Ni  were 
deposited  via  argon  sputtering,  followed  by  500  A  Au  deposited  via  thermal  evaporation.  The  argon 
sputtering  resulted  in  large  grains  of  Ni  being  deposited  on  the  crystalline  SiC  substrate  (Fig.  4.22). 
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Subsequent  annealing  resulted  in  a  granular  contact  (Fig.  4.23).  These  contacts  were  Schottky  when 
unannealed  and  entirely  prevented  any  current  flow  after  annealing  for  1  min  or  more  at  900  °C. 

4.3.3. 1.2  Porous  SiC 

Upon  selection  of  Ti/Ni/Au  and  Ti/NiCr/Au  as  the  best  contact  materials,  ohmic  contact  to  porous 
SiC  layers  was  investigated.  As  detailed  in  section  4. 3.2.2,  the  success  rate  of  isolated  TLM  pad 
formation  on  porous  SiC  was  low,  due  to  the  inability  to  use  sonication  to  aid  lift  off  of  excess  metal.  For 
this  reason,  TLM  studies  on  porous  SiC  were  only  conducted  for  the  Ti/Ni/Au  contacts. 

The  contact  was  found  to  have  good  adhesion  to  the  porous  substrate.  The  unannealed  contact 
fills  the  top  layer  of  pores  well  (Fig.  4.24).  Upon  annealing,  the  contact  takes  on  a  roughened  surface  and 
continues  to  fill  the  top  layer  of  the  porous  morphology  (Fig.  4.25).  In  fact,  the  contact  adheres  to  the 
PSiC  so  strongly  that  sonication  first  leads  to  fracture  within  the  porous  film,  while  the  contact  retains  its 
adhesion  to  the  top  of  the  porous  layer  (Fig.  4.26).  Annealing  at  900-1000  °C  does  not  appear  to  change 
the  morphology  of  the  PSiC  film. 

When  making  porous  SiC  TLM  structures,  a  complete  row  of  ohmic  contact  pads  was  never 
fabricated.  Rather,  some  rows  would  have  several  good  pads  of  one  spacing,  another  row  would  have  a 
couple  of  good  pads  of  another  spacing,  etc.  Therefore,  specific  contact  resistance  data  for  the  porous 
samples  was  determined  by  average  resistance  values  from  numerous  different  rows  in  order  to 
accumulate  enough  data.  Raw  data  for  the  resistances  found  with  porous  SiC  are  given  in  Table  4.3. 
Using  the  data  presented  in  Table  4.3,  the  specific  contact  resistance  of  the  contacts  to  porous  SiC  is 
found  to  be  5.5  x  10  4  Q  cm2  (Fig.  4.27). 

4.3.3.2  Device  Behavior 

Crystalline  devices  displayed  total  resistances  across  the  wafer  always  near  1  -3  Q.  Of  the  porous 
devices  tested,  total  resistances  showed  greater  variance.  When  the  contact  was  formed  from  Ti/Ni/Au, 
measured  total  resistances  included  3  Q,  40  Q,  55  Q,  161  Q ,  and  185  Q.  The  one  device  formed  from  a 
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Figure  4.21.  I/V  curves  for  contacts  composed  of  1500  A  NiCr  and  500  A  Au  annealed  at  900  °C  for  the 
times  indicated.  The  substrate  is  crystalline  SiC.  (a)  is  shown  in  full  scale,  while  (b)  has  the  current  scale 
enlarged  to  show  detail. 
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Figure  4.22.  Cross-sectional  SEM  image  of  an  unannealed  contact  composed  of  argon  sputtered  A1  (200 
A)  and  Ni  (1250  A)  followed  by  500  A  Au  deposited  via  thermal  evaporation.  The  substrate  is  crystalline 
SiC, 


Figure  4.23.  Cross-sectional  SEM  image  of  a  contact  composed  of  argon  sputtered  A1  (200  A)  and  Ni 
(1250  A)  followed  by  500  A  Au  deposited  via  thermal  evaporation.  The  contact  has  been  annealed  at  900 
°C  for  1  min.  The  substrate  is  crystalline  SiC. 
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Figure  4.24.  Cross-sectional  SEM  image  of  an  unannealed  contact  composed  of  50  A  Ti,  1500  A  Ni,  and 
500  A  Au.  The  substrate  is  porous  SiC. 


Figure  4.25.  Cross-sectional  SEM  image  of  a  contact  composed  of  50  A  Ti,  1500  A  Ni,  and  500  A  Au, 
annealed  at  900  °C  for  5  min.  The  substrate  is  porous  SiC.  Notice  the  macroscopic  adhesion  of  the  metal 
contact  to  the  porous  morphology. 
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Figure  4.26.  Cross-sectional  SEM  image  of  a  contact  composed  of  50  A  Ti,  1500  A  Ni,  and  500  A  Au 
annealed  at  900  °C  for  5  min.  The  substrate  is  porous  SiC,  and  the  sample  has  been  sonicated  to  aid  in 
metal  lift  off.  Notice  how  sonication  has  lead  to  fracture  of  the  porous  layer,  but  the  contact  remains 
adhered  to  the  porous  morphology. 


Pad  Spacing,  d  (microns) 

Figure  4.27.  TLM  data  collected  from  a  row  of  contact  pads  composed  of  50  A  Ti,  1500  A  Ni,  and  500  A 
Au  annealed  at  900  °C  for  10  min.  The  substrate  is  porous  SiC.  Analysis  of  the  data  according  to  Eqs. 

4. 2-4. 5  yields  a  specific  contact  resistance  of  5.5  x  10'4  F2  cm2  for  these  contacts. 
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Table  4.3.  Raw  Data  for  the  Resistance  of  Different  Pad  Spacings 
During  TLM  Studies  of  Ti/Ni/Au  Contacts  on  Porous  SiC 


Pad  Spacing  ( nm)  Measured  Resistance  (Q.) 


2  144,  165,  149,  143 

3  153,  157,  126,  159,  163 

4  125,154,184,134,115 

5  165,  143,  170,  174,  146 

6  146,  156,  187,  185,  183 

8  187,207,203,209,194 

10  203,211,171,189,204 


Ti/NiCr/Au  contact  had  a  total  resistance  of  6  Q.  Asa  source  of  comparison,  a  Schottky  (unannealed) 
Ti/Ni/Au  contact  showed  a  resistance  of  285  Q  across  a  PSiC  device  at  300  mV  applied  potential. 


When  exposed  to  a  constant  (250  seem)  flow  of  nitrogen  gas  and  heated  from  25-500  °C, 
crystalline  and  porous  SiC  devices  behaved  as  shown  in  Table  4.4. 

Table  4.4  Response  of  SiC  Sensors  to  Temperature  Increases  in  an  Ambient  of  Pure  N2 


Sample 

Contact 

Charm c  in  resistance 

25  -  250  °C 

25  -500  -C 

Crystalline  SiC 

(1)  Ohmic  Ti/Ni/Au 

(1)  0.9%  decrease 

(1)  4.8%  decrease 

(2)  Ohmic  Ti/Ni/Au 

(2)  6.3%  decrease 

(2)  12.0%  decrease 

(3)  Ohmic  Ti/Ni/Au 

(3)  0.3%  increase 

(3)  4.7%  increase 

Porous  SiC 

(1)  Ohmic  Ti/Ni/Au 

(1)  15.1%  decrease 

(1)  45.5%  decrease 

(2)  Ohmic  Ti/Ni/Au 

(2)  29.4%  decrease 

(2)  42.7%  decrease 

(3)  Ohmic  Ti/NiCr/Au 

(3)  20.7%  decrease 

(3)  34.5%  decrease 

The  data  in  Table  4.4  indicates  that  the  SiC  wafer  is  behaving  as  an  intrinsic  semiconductor  for 
the  most  part  (except  for  the  one  crystalline  sample  that  experienced  an  increase  in  resistance  upon 
heating).  This  behavior  is  to  be  expected  if  the  donor  levels  are  relatively  shallow  (<  25  meV),  as  nearly 
all  of  the  donor  electrons  would  already  be  in  the  conduction  band  at  room  temperature  and  therefore 
intrinsic  behavior  would  result.  Previous  results  in  the  literature  support  this  possibility.85'86  Moore85 
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shows  that  the  binding  energy  of  nitrogen  donors  in  cubic  SiC  is  54  meV,  as  confirmed  by 
photoluminescence  studies.  However,  conductivity  studies  confirm  the  existence  of  another  donor  at  <  20 
meV,  which  could  arise  from  a  variety  of  sources.  First,  the  shallow  donor  energy  level  may  be  due  to 
higher  concentration  pockets  of  nitrogen  donors.  Lin87  showed  that  at  concentrations  higher  than  2.5  x 
1018  cm'3  (slightly  less  doped  than  our  wafer),  arsenic  in  silicon  (whose  binding  energy  approximates  that 
of  N  in  SiC)  showed  a  similar  reduction  in  binding  energy.  Alternatively,  the  shallow  donor  level  could 
be  due  to  the  existence  of  D'  (overcharged  nitrogen)  centers  or  donor  complexes.  Nitrogen  in  the  vicinity 
of  inversion  domain  boundaries  or  some  other  source  of  excess  electrons  could  create  the  necessary  D‘ 
centers.85 

The  data  in  Table  4.4  also  indicates  that  the  resistance  of  the  porous  devices  is  much  more 
sensitive  to  changes  in  the  temperature  of  the  surrounding  gas.  One  possible  explanation  for  this 
phenomenon  is  the  presence  of  deep  (on  the  order  of  a  few  hundred  meV  below  the  conduction  band 
edge)  traps.  It  is  entirely  plausible  that  the  crystalline  material  has  very  few  deep  traps,  but  that  the 
porous  material  has  many  more  as  a  result  of  the  anodization  process.  If  the  electrochemical  etching 
produces  an  abundance  of  deep  traps  that  are  significantly  populated  within  the  porous  material,  then  the 
promotion  of  electrons  from  these  trap  states  into  the  conduction  band  could  account  for  the  decreased 
resistance  in  the  porous  material  upon  heating.  Regardless  of  the  cause,  PSiC  has  a  much  larger 
temperature  coefficient  of  conductivity  than  crystalline  SiC,  and  this  may  eventually  open  up  new 
avenues  of  temperature  sensing  with  PSiC. 

Conductivity  changes  as  a  function  of  exposure  to  hydrogen,  acetone,  ethanol,  and  pure  nitrogen 
are  presented  in  Table  4.5  and  Table  4.6  (the  hydrogen,  acetone,  and  ethanol  are  diluted  in  N2  to  the 
concentrations  indicated).  An  example  of  some  representative  data  is  presented  in  Figs.  4.28  and  4.29. 
The  acetone,  ethanol,  and  hydrogen  tests  were  conducted  at  total  flow  rates  of  1 00  seem.  Each  individual 
resistance  measurement  has  error  bars  of  ±0.6%. 


Ill 


Voltage  (mV) 


Figure  4.28.  An  example  of  conductometric  data  collection.  This  typical  data  shows  the  change  in 
resistance  of  a  PSiC  sensor  when  exposed  from  an  increasing  concentration  of  0-20%  1 12.  The  F^  is 
diluted  in  N2  to  a  total  flow  rate  of  100  seem;  the  wafer  is  at  250  °C. 
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Figure  4.29.  An  example  of  conductometric  data  collection.  This  typical  data  shows  the  change  in 
resistance  of  a  PSiC  sensor  when  exposed  from  an  increasing  concentration  of  0-20%  142,  followed  by 
decreasing  the  concentration  from  20-0%  142.  4'he  142  is  diluted  in  N2  to  a  total  flow  rate  of  100  seem;  the 
wafer  is  at  250  °C. 
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Table  4.5.  Summary  of  Tests  in  which  the  Device  was  Maintained  at  a  Constant  Temperature; 

Data  from  Unique  Devices  is  Denoted  by  (#) 

Crystalline  SiC 


Concentration 


Gas 

range  (at%) 

Room  temp 

250  -C 

500  2C 

Acetone 

0-0.85% 

No  response 

No  response 

No  response 

Ethanol 

0-0.65% 

No  response 

No  response 

No  response 

Hydrogen 

0-1.6% 

No  response 

No  response 

No  response 

Hydrogen 

0-20% 

No  response 

No  response 

No  response 

Porous  SiC 

Acetone 

0-0.85% 

No  response 

No  response 

No  response 

Ethanol 

0-0.65% 

No  response 

No  response 

No  response 

Hydrogen 

0-1.6% 

No  response 

No  response 

No  response 

Hydrogen 

0-20% 

No  response 

(1)  No  response 

(1)  AR  =  -5.1%(145mQ) 

(2)  No  response 

(2)  AR  =  -7.0%  (321  mQ) 

(3)  AR  =  -1.2%  (50  mfi)  (3)  AR  =  -5.6%  (224  mQ) 

Table 

4.6.  Summary  of  Tests  in  which  the  Device  was  Supplied  a  Constant  Applied  Heat  (Device 
Temperature  was  Allowed  to  Change  as  the  Gas  Flow  was  Varied); 

Data  from  Unique  Devices  is  Denoted  by  (#) 

Crystalline  SiC 

Gas 

Concentration 
range  (at%) 

-250 2C 

-500  2C 

Acetone 

0-0.85% 

AR  =  +0.4%  (7  mO); 

AT  =  +4.3  °C 

AR  =  -0.7%  (1 1  mQ);  AT  =  0 

Ethanol 

0-0.65% 

AR  =  +2.4%  (42  mO); 

AT  =  +6.1  °C 

AR  =  -1.0%  (17  mQ);  AT  =  0 

Hydrogen 

0-20% 

AR  =  -1.3%  (23  mfl);  AT  =  +18.6  °C 

AR  =  -3.5%  (60  mQ); 

AT  =  +31.1  °C 

Porous  SiC 

Acetone 

0-0.85% 

No  response,  AT  =  +3.8  °C 

AR  =  -1.5%  (32  mQ);  AT  =  0 

Ethanol 

0-0.65% 

No  response,  AT  =  +3.4  °C 

No  response;  AT  =  0 
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Hydrogen  0-1.6% 


No  response,  AT  =  0 


(1)  AR  =  -0.9%  (29  mO); 
AT  =  0 


(2)  AR  =  -0.5%  (16  mQ); 

AT  =  +3.2  °C 

Hydrogen  0-20%  (1)  No  response,  AT  = +44.8  °C  (1)  AR  =  -21%  (388  mQ); 

AT  =  +28  °C 

(2)  AR  =  -1.2%  (5 1  mQ);  (2)  AR  =  -8.9%  (227  mQ); 

AT  =  +17  °C  AT  =  +26  °C 

As  seen  in  Tables  4.5  and  4.6,  none  of  the  devices  were  able  to  sense  any  analyte  flow  at  room 
temperature.  In  addition,  even  at  high  temperature,  neither  acetone  nor  ethanol  causes  changes  in  the  total 
resistance  much  greater  than  the  noise  for  the  measurements.  Of  course,  as  shown  above,  only  small 
concentrations  of  ethanol  and  acetone  have  been  investigated,  as  these  initial  experiments  were  conducted 
in  order  to  look  for  low  detection  limits  comparable  to  the  results  found  from  the  Schottky  diode  sensor 
results  in  the  literature.  Our  later  studies  focused  on  hydrogen  at  a  wide  range  of  concentrations  due  to  its 
economic  and  safety  considerations  for  NASA  operations. 

Examining  the  hydrogen  data  collected  at  constant  temperature  (Table  4.5),  the  crystalline  device 
is  unable  to  sense  H2  at  any  temperature  investigated  in  this  work.  However,  the  PSiC  sensor  begins  to 
reproducibly  detect  H2  once  the  system  is  heated  to  500  °C.  The  experiments  conducted  at  constant 
applied  heat  (Table  4.6)  verify  this  result,  although  the  constant  heat  tests  show  slightly  greater  decreases 
in  resistance,  probably  due  to  the  additive  effects  of  convective  wafer  heating  (H2  has  a  higher  heat 
capacity  than  N2)  and  changes  in  conduction. 

Given  that  20%  H2  in  the  gas  flow  yields  approximately  a  6%  decrease  in  total  resistance,  and 
resistance  measurements  are  ±0.6%,  a  linear  response  of  the  porous  sensors  would  yield  an  estimated 
minimum  detection  limit  of  2%  H2  at  500  °C.  Consistent  with  this  estimate  the  sensors  were  unable  to 
detect  1.6%  H2  at  500  °C,  yet  reliably  and  reproducibly  measured  a  decrease  in  resistance  when  exposed 
to  2.5%  H2  (the  first  data  point  in  the  0-20%  H2  tests,  as  shown  for  example  in  Fig.  4.29). 
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One  device  was  made  using  the  same  contact  recipe,  but  was  not  annealed,  therefore  yielding  a 
Schottky  contact.  For  a  300  mV  potential,  this  device  yielded  a  380%  increase  in  current  when  heated 
from  25-500  °C.  This  is  a  relatively  large  increase  in  signal  strength  when  compared  to  the  results  from 
the  ohmic  contacts  presented  above  (Table  4.4).  Flowever,  it  was  found  that  the  current  through  the 
Schottky  device  changed  (decreased)  from  experiment  to  experiment,  and  it  is  possible  that  the  heating 
cycles  were  instigating  interfacial  reactions  within  the  contact  layers.  Like  the  annealed,  ohmic  contacts, 
this  device  did  not  sense  FL  flow  at  250  °C.  The  Schottky  diode  did  experience  a  14%  increase  in  current, 
though,  when  exposed  to  20%  F^  at  500  °C.  This  result  is  greater  than  the  results  found  for  the  ohmic 
contacts  above.  The  increase  in  current  could  be  due  to  the  combined  effects  of  H2  adsoiption  to  the 
semiconducting  skeleton  surface  as  well  as  lowering  the  Schottky  barrier  at  the  contact/SiC  interface. 
Flowever,  upon  a  few  hours  of  exposure  to  500  °C,  the  Schottky  contact  was  found  to  fail  and  suffer  a 
dramatic  decrease  in  current. 

Wide  voltage  sweeps  (-2.5  to  2.5  V)  were  conducted  at  room  temperature  to  compare  the 
ohmicity  of  the  various  contacts  for  a  wide  range  of  potentials  (Fig.  4.30).  As  expected,  the  ohmic 
contacts  to  the  crystalline  devices  have  the  most  linear  and  steepest  slopes.  The  ohmic  contacts  to  the 
porous  devices  show  good  linearity,  too,  although  the  slope  is  less  than  the  crystalline  ohmic  contacts 
(indicating  the  higher  resistances  described  in  the  beginning  of  this  section).  As  expected,  the  Schottky 
(unannealed)  contact  to  the  porous  device  shows  a  dramatically  decreased  current  flow  and  non-linearity. 
4.3.4  Conclusions 

These  experiments  do  not  address  whether  hydrogen  absoiption  into  the  SiC  crystal  lattice  is  the 
root  cause  of  the  changes  in  conductivity  found  during  these  conductometric  tests.  Flowever,  studies  of 
hydrogen  passivation  of  SiC  support  this  interpretation.  When  hydrogen  absorption  into  the  wafer  is 
conducted  at  high  pressure  (10  bars)  accompanied  by  high  temperatures  (1500-1700  °C),  quenching 
results  in  hydrogen-nitrogen  complexes  (passivation),  which  decrease  the  free  carrier  concentration  of  n- 
type  SiC  wafers, 8S  because  the  hydrogen  bonds  with  the  extra  electron  on  the  donor  nitrogen  atoms.  The 
decrease  in  the  free  carrier  concentration  should  lead  to  an  increase  in  resistance,  not  the  decrease  in 
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Figure  4.30.  Large  sweeps  of  current,  demonstrating  the  properties  of  ohmic  and  Schottky  contacts  to 
both  crystalline  and  porous  SiC. 

resistance  observed  in  these  studies.  However,  studies  of  SiC  exposed  to  H2  gas  at  500  °C  (a  temperature 
closer  to  those  used  in  this  study),  show  that  hydrogen  is  present  in  the  bulk  of  the  SiC,  but  that  the 
hydrogen-nitrogen  binding  reaction  is  thermodynamically  or  kinetically  impeded.89  Furthermore,  in  n+ 
material,  such  as  the  wafers  used  in  this  study,  H'  is  expected  to  be  the  dominant  species,  because  the 
Fermi  level  is  high  and  therefore  the  probability  for  the  hydrogen  acceptor  level  to  be  occupied  is  large.89 
Finally,  it  has  been  shown  that  atomic  hydrogen  is  mobile  in  SiC  even  down  to  room  temperature.90  As  a 
result,  it  is  reasonable  to  surmise  that  ambient  molecular  hydrogen  may  be  dissociated  and  absorbed  into 
the  porous  semiconducting  skeleton,  and  that  hydrogen  may  be  free  to  serve  as  an  additional  charge 
carrier. 

This  work  indicates  that  PSiC  conductometric  sensors  cannot  compete  with  the  low  detection 
limits  realized  with  crystalline  SiC  Schottky  diode  devices.  However,  the  PSiC  sensors  may  be  able  to 
fill  a  niche  in  a  sensor  array  to  quantify  relatively  high  (>  2%)  H2  concentrations.  As  indicated  earlier  the 
hydrogen  sensors  which  respond  to  ppm-level  hydrogen  loads  saturate  at  levels  far  below  the  2-20% 
range  of  H2  concentration  to  which  the  PSiC  sensors  respond. 
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There  are  bulk  sensors  that  are  of  semiconductor  oxide  or  “hot-wire”  type.  In  these  sensors,  the 
resistance  of  the  sensor  materials  changes  as  the  sensor  surface  is  reduced.  One  architecture  reported  in 
the  literature  utilized  Sn02  sintered  in  a  bead  covering  a  Pt  wire  coil.91  Due  to  reduction  at  the  Sn02 
surface,  the  resistance  of  the  Pt  wire  changes  as  the  sensor  is  exposed  to  H2.  However,  these  devices  still 
saturate  at  relatively  low  (<  1%)  H2  concentrations. 

Therefore,  planar  thermoelectric  gas  sensors  have  also  been  investigated.92'96  Utilizing  alumina 
substrates,  these  sensors  are  fabricated  by  partially  coating  a  nickel  oxide  thick  (~700  nm)  film  with  a  thin 
(~75  nm)  film  of  Pt.  Catalytic  reactions  of  the  Pt  with  hydrogen  heats  the  platinum  coated  surface,  and 
then  the  thermoelectric  voltage  builds  up  along  the  hot  and  cold  region  of  the  oxide  film.  These  sensors 
have  been  shown  to  detect  0.025-10%  H2  in  air  at  80  °C  via  relatively  large  (>  100%)  changes  in  voltage. 
However,  temperature  control  of  these  devices  could  be  an  issue  during  implementation.  For  example, 

3%  H2  in  air  shows  only  a  0.12  °C  temperature  difference  (and  0.087  mV  signal)  at  room  temperature.92 
In  addition,  the  catalytic  activity  is  highly  dependent  upon  the  oxidation  state  and  the  degree  of 
crystallization  of  the  thermoelectric  materials.  Finally,  weak  adhesion  between  the  thermoelectric  thick 
film  and  the  substrate  has  been  a  problem. 

In  order  to  eliminate  the  processing  concerns  associated  with  the  production  of  the  nickel  oxide 
layer,  thermoelectric  gas  sensors  have  also  been  formed  by  partially  coating  a  Sio.xGco.2  thin  film  with 
Pt.97'98  This  device  exhibited  a  6-10  mV  potential  increase  when  exposed  to  10%  H2.  However,  the 
device  has  been  researched  primarily  at  temperatures  below  150  °C,  because  the  device  is  selective 
toward  H2  below  150  °C,  but  showed  response  to  several  other  gases  (methanol,  ethanol,  butane,  CO, 
methane)  at  higher  temperatures. 

In  current  solid  state  hydrogen  detection  systems,  a  hydrogen  sensitive  resistor  provides  detection 
of  higher  hydrogen  concentrations  (greater  than  approximately  1  %  H2)  while  a  Schottky  diode  quantifies 
lower  hydrogen  concentrations.99  Our  results  indicate  that  PSiC  sensors  could  fill  the  void  in  the  array  of 
sensors  for  the  relatively  high  concentrations  of  hydrogen  detection. 
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In  the  early  1990’s,  thin  film  hydrogen  resistors  were  composed  of  Pd/Ag  or  later  Pd/Ni.  The 


resistances  of  these  metal  films  have  a  response  proportional  to  the  partial  pressure  of  H2:  [  P/,  j  .  This 

dependence  is  due  to  the  sensor  detection  mechanism:  migration  of  hydrogen  into  the  bulk  of  the  metal 
changing  the  bulk  conductance  of  the  metal.100  These  sensors  generally  have  low  sensitivity  at  low 
hydrogen  concentrations  but  a  continued  response  over  a  wide  range  of  hydrogen  concentrations.  Pd/Ni 
thin  films  were  made  that  showed  ~3%  decrease  in  resistance  when  exposed  to  10%  H2.101 

Current  thin  film  hydrogen  resistors  are  usually  composed  of  palladium  chrome  or  palladium 
titanium  alloy,  as  the  alloys  show  a  greater  resistance  to  phase  transformation  during  hydrogen 
absorption.99,102  The  palladium  titanium  alloys  (71-99  at%  Pd)  have  shown  an  increase  in  resistance  of  up 
to  16%  when  exposed  to  100%  hydrogen,  and  an  increase  in  resistance  of  5%  when  exposed  to  an 
environment  of  10%  hydrogen.102  Therefore,  it  can  be  seen  that  the  results  obtained  in  this  dissertation 
show  magnitudes  of  resistance  change  slightly  less  than  current  high  concentration  hydrogen  sensing  via 
thin  metallic  films. 
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CHAPTER  5 


FUTURE  DIRECTIONS 


5.1  Electroless  Etching 

There  are  many  etching  parameters  that  can  be  manipulated  during  electroless  etching,  such  as 
the  identity  of  the  metal  catalyst,  the  etchant  composition,  etc.  However,  the  primary  factor  that 
determines  the  morphology  of  electrolessly  etched  PSiC  is  determined  before  etching  even  begins:  the 
quality  of  the  starting  SiC  wafer.  Current  SiC  production  is  simply  still  in  the  research  phase,  and  the 
growth  of  pure,  defect-controlled  SiC  is  not  available  with  today’s  technology  and  growth  techniques.  An 
example  of  current  deficiencies  is  our  4H-SiC  wafer  with  its  preponderance  of  6H  regions.  In  addition, 
dopant  concentrations  often  differ  by  more  than  an  order  of  magnitude  between  the  center  of  a  2”  wafer 
and  its  edges.  Relatively  high  concentrations  of  dopants  such  as  boron  and  nitrogen  are  present  in  all 
commercially  available  wafers  today  except  for  the  purist  (and  most  expensive)  of  growth  conditions. 

But  the  wafer  characteristic  that  appears  to  most  determine  the  resulting  PSiC  morphology  is 
defect  structure  and  concentration.  The  concentration  of  micropipes  has  decreased  dramatically  in  recent 
years,  but  even  now  micropipe  densities  on  the  order  of  100  pipes/cm2  or  less  are  difficult  to  obtain  in 
bulk  wafers.  Simple  edge  and  screw  dislocation  densities  are  much  higher  than  this  number,  and  are 
much  greater  than  the  dislocation  densities  obtained  in  the  mature  Si  production.  This  research  showed 
that  the  morphology  of  the  porous  structure  is  highly  dependent  upon  the  micropipes  and  defects  (both 
stacking  faults  and  threading  dislocations)  in  the  SiC  wafer.  Once  the  impurity  fluctuations  and  defect 
concentrations  in  crystalline  SiC  wafers  can  be  controlled  and  minimized,  electroless  etching  of  PSiC  can 
begin  in  earnest. 

5.2  Further  Characterization  of  Porous  SiC 

An  additional  spectroscopic  characterization  technique  of  PSiC  that  could  prove  useful  is  near¬ 
field  scanning  optical  microscopy  (NSOM).  NSOM  would  allow  for  simultaneous  correlation  of  light 
emission  and  scattering  properties  with  the  morphology  of  the  porosity.  Therefore,  further  correlations 
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could  be  made  in  regard  to  the  origin  of  spectral  emission  from  PSiC.  In  addition,  NSOM  would  permit 
the  collection  of  photoluminescent  data  that  is  much  more  localized  to  the  surface  region  of  the  PSiC  film. 

Two-phonon  Raman  studies  may  also  prove  useful  if  porous  SiC  layers  are  ever  developed  into 
temperature  sensors.1’3  Two-phonon  Raman  experiments  would  probe  the  local  phonon  density  of  states, 
which  directly  determines  the  thermal  conductivity  of  the  material.  Thermal  conductivity  is  obviously  an 
important  material  parameter  to  know  accurately  in  order  to  form  a  temperature  sensor.  Therefore,  the 
phonon  density  of  states  could  be  mapped  as  a  function  of  etching  conditions.  This  information  would  be 
particularly  useful  for  studies  being  performed  by  a  collaborator  at  UIUC,  in  which  the  flexing  of 
interdigitated  arrays  of  porous  SiC  cantilevers  is  being  monitored  as  a  temperature  detector. 

5.3  High  Temperature  Sensing  with  PSiC 

All  of  the  porous  sensors  in  this  work  were  anodically  etched  with  the  same  conditions  (current, 
time,  and  therefore  porous  film  thickness).  Experiments  should  be  conducted  in  which  the  etching  current 
is  varied  between  sensors,  as  varying  the  current  could  alter  both  the  porous  morphology  to  a  limited 
degree  as  well  as  possibly  the  surface  termination  (including  the  presence  of  dangling  bonds).  In 
addition,  anodically  etching  for  longer  times  in  order  to  produce  thicker  porous  films  could  lower  the 
detection  limit  of  these  sensors. 

Making  ohmic  electrical  contact  directly  to  porous  SiC  is  still  an  art  in  need  of  improvement.  Lift 
off  of  excess  metal  continued  to  be  a  problem  throughout  this  work.  There  are  a  variety  of  processing 
parameters  that  can  be  manipulated  during  contact  formation;  however,  the  single  most  useful  alteration 
would  probably  be  the  use  of  a  thick  photoresist.  Coupled  with  the  correct  developer  to  form  sharp  (“boat 
dock”)  features,  lift  off  of  excess  metal  from  porous  regions  could  be  greatly  enhanced.  When  one  of  the 
devices  that  had  poor  lift  off  was  examined  in  SEM,  it  was  found  that  the  photoresist  had  dissolved,  but 
that  the  overlying  metal  was  so  thick  that  the  metal  remained  suspended  over  an  area  that  should  have 
become  bare  (Fig.  5.1). 

The  initial  response  time  as  well  as  the  relaxation  time  of  the  sensor  to  nominal  conditions  after 
the  H2  flow  was  stopped  was  not  investigated  in  this  work.  Rather,  as  the  H2  concentration  was 
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Metal  arching  over  the  sensor  region 
(should  have  been  removed  during  liftoff) 


Former  location  of  PR 


15.0kV  10.1mm  x13.0k  6/1/04 


Metal  deposited 
directly  on  PSiC 


1  II  I  I  l  I  I  I  iT 

2.00um 


Fig.  5.1  Cross-sectional  SEM  image  of  a  contact  composed  of  50  A  Ti,  1500  A  Ni,  and  2500  A  Au.  The 
region  on  the  right  of  the  image  is  the  location  where  the  contact  metal  is  bonded  directly  to  the  PSiC  and 
should  remain  after  lift  off.  The  region  in  the  left  of  the  image  is  the  location  in  which  photoresist  existed 
so  that  lift  off  would  remove  this  metal.  The  metallization  is  too  thick  for  the  photoresist  and  has  bridged 
the  gap  where  the  photoresist  dissolved  away  during  lift  off. 


manipulated,  the  sensor  was  simply  given  an  adequate  (5  min)  time  to  equilibrate  to  the  new  gas  flow 
condition.  Transient  condition  time  responses  could  be  conducted  by  applying  a  steady  voltage  and 
monitoring  the  change  in  current  with  time.  This  would  require  the  construction  of  a  smaller  flow  cell 
that  could  allow  for  immediate  exposure  of  the  sensor  to  the  new  flow  conditions,  as  the  relatively  large 
flow  cell  used  in  this  work  was  designed  for  simultaneous  interferometry  and  conductivity  studies. 

Results  of  these  studies  could  also  elucidate  whether  changes  in  conductivity  are  due  to  adsorption  to  the 
surface  or  absoiption  into  the  lattice. 

As  noted  in  Chapter  4,  the  Schottky  diode  sensors  detailed  in  the  literature  behave  quite 
differently  depending  upon  whether  the  gas  flow  is  fuel  rich  or  fuel  lean.  No  experiments  were  conducted 
in  this  work  with  02  in  the  flow  with  the  H2.  The  presence  of  oxygen  in  the  gas  flow  may  not  change  the 
response  of  these  PSiC  sensors  considering  that  (1)  a  stable  native  oxide  already  exists  on  the  sensors 
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investigated  in  this  work,  and  (2)  these  sensors  proved  insensitive  to  low  concentrations  of  acetone  and 
ethanol.  Experiments  with  both  fuel  rich  and  fuel  lean  conditions  should  be  conducted  to  verify  this 
hypothesis. 

Furthermore,  experiments  with  higher  concentrations  of  other  hydrogen-containing  species  (such 
as  ethanol)  should  be  conducted.  This  work  only  investigated  low  concentrations  (<  1%)  of  ethanol  and 
acetone  in  order  to  compare  the  response  of  the  PSiC  sensors  to  the  Schottky  sensors  detailed  in  the 
literature.  Flowever,  as  our  future  FE  tests  showed,  the  PSiC  sensors  respond  to  only  relatively  large 
concentrations  of  ambient  hydrogen.  The  PSiC  sensors  may  also  be  able  to  detect  relatively  large 
concentrations  of  hydrocarbons. 

Buriak  and  Allen  showed  that  the  hydrosilylation  of  alkynes  and  alkenes  to  porous  Si  produced 
remarkably  stable  surfaces.4  The  PSi  surface  was  partially  covered  with  the  covalently  bound  vinyPalkyl 
groups,  and  the  original  Si-FT  porous  surface  termination  remained  to  some  degree.  Flowever,  the  Si-H 
surface  appeared  to  be  effectively  capped,  as  exposure  of  the  PSi  to  F1F  and  KOFI  did  not  result  in  oxide 
formation  or  dissolution  of  the  porous  film.  A  similar  surface  treatment  may  be  beneficial  in  producing 
PSiC  sensors.  It  may  be  possible  to  anodically  etch  a  bulk  SiC  wafer  to  form  PSiC,  immediately  followed 
by  a  surface  treatment  that  continues  to  permit  FE  adsorption  on  the  surface,  yet  is  resistant  to  oxidation. 

A  PSiC  surface  that  is  resistant  to  oxidation,  yet  still  responsive  to  FE  in  the  surrounding  environment, 
would  be  an  extremely  powerful  sensor  for  the  ever-changing  fuePoxygen  ratios  found  in  an  engine’s 
exhaust. 
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APPENDIX  A 


QUANTUM  CONFINEMENT  FIT  MATHEMATICA  CODE 

C  =  1.5117;  (*  NANOMETERS  *) 

GAMMA0=  11.80;  (*  INVERSE  CM  *) 

«  STATISTICS'DATASMOOTHING'; 

OMEGASQ[Q  J  :=  657500  +  270000  COS[PI  Q/2]; 

OMEGA[Q  J  :=  SQRT[OMEGASQ[Q]]; 

CC0[Q_,  LJ  :=  EXP[-QA2  LA2  /(4*CA2)];  (*  L  IN  UNITS  OF  NANOMETERS  *) 
INTEGRAND [Q_,  W_,  L  J  :=  4  PI  (QA2)  CC0[Q,  L]/((W  -  OMEGA[Q])A2  +  ( GAMMA0/2)A2) ; 
INTENSITY  [W_,  LJ  :=  NINTEGRATE[INTEGRAND[Q,  W,  L],  {Q,  0,  1}]; 

(*  FUNCTIONS  TO  SELECT  AND  SMOOTH  DATA  *) 

TERMFIND  [D  ATA_,  VALUE  ]  :=  ( 

FOR[JJ  =  1,  JJ  <  LENGTH[DATA],  JJ++, 

IF[DATA[[JJ]] [[  1  ]]  <  VALUE,  COUNTER  =  JJ,  BREAK] 

]; 

COUNTER 

); 

PEAKFIND[DATA  J  :=  {DATA[[(POSITION[DATA,  MAX[DATA]][[1]][[1]])]][[1]], 
DATA[[(POSITION[DATA,  MAX[DATA]][[1]][[1]])]][[2]]}; 
SMOOTHDATA[DATA_,  TERMS_,  FIRST_,  LAST  ]  :=  ( 

SMOOTH  =  MOVING  AVERAGE  [DATA,  TERMS]; 

TABLE[SMOOTH[[N]],  {N,  FIRST,  LAST}] 

); 

SMOOTHWAVE [DATA_,  TERMS_,  LOWLAM_,  HIGHLAM J  :=  ( 
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FIRST  =  TERMFIND  [DATA,  LOWLAM]; 


LAST  =  TERMFIND  [DATA,  HIGHLAM]  +  1; 

SMOOTHDATA[DATA,  TERMS,  FIRST,  LAST] 

); 

SKIMDATA[DATA_,  NUMBER  ]  :=  ( 

FCN[XJ  :=  INTERPOLATION  [DATA]  [X] ; 

LEN  =  LENGTH[DATA]; 

DX  =  (DATA[[LEN]][[1]]  -  DATA[[1]][[1]])/NUMBER; 

XX [N_]  :=  DATA[[1]][[1]]  +  N*DX; 

TABLE[{XX[N],  FCN[XX[N]]},  {N,  0,  NUMBER}] 

); 

THRESHOLD[DATA_,  OFFSET  ]  :=  MIN[PART[TRANSPOSE[DATA],  2]]  -  OFFSET; 
TREATDATA [DATA_,  OFFSET  ]  :=  ( 

THRHLD  =  THRESHOLD  [DATA,  OFFSET]; 

TABLE[{DATA[[N]][[1]],  DATA[[N]][[2]]  -  THRHLD},  {N,  1,  LENGTH[DATA]}] 

); 


(*  INTEGRATES  OVER  DATA  *) 

DATASUM[DATA_]  :=  SUM[DATA[[N]][[2]]*(DATA[[LENGTH[DATA]]][[1]]  - 
DATA  [  [  1  ]  ]  [  [  1  ]  ])/LEN  GTH  [DAT  A] ,  {N,  1,  LENGTH  [DATA]}]; 

(*  BESTLPEAK  FINDS  BEST  L  VALUE  BASED  ON  POSITION  OF  PEAK  IN  DATA  *) 
BESTLPEAK[DATA_]  :=  ( 

PEAKLAMBDA  =  PEAKFIND  [DATA]  [  [  1  ]] ; 

LOWLAM  =  DATA[[1]][[1]]; 

HIGHLAM  =  D ATA  [  [LEN GTH  [DATA]  ]][[!]]; 
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L  /.  FINDROOT[PEAKLAMBDA  == 

(W/.FINDMINIMUM[1/INTENSITY[W,  L],  {W,  LOWLAM,  HIGHLAM}][[2]]),  {L,  1,  10}] 

); 

COMPAREPEAK[DATA_]  :=  ( 

PLT1  =  LISTPLOT[DATA]; 

LVAL  =  BESTLPEAK[DATA]; 

PRINT["BEST  FITTED  L=",  LVAL]; 

PEAKDATA  =  PEAKFIND[DATA]; 

NORMDATA  =  PEAKDATA[[2]]; 

LEN  =  LENGTH[DATA]; 

NORMFIT  =  INTENSITY[W  =  PEAKDATA[[1]],  LVAL]; 

PLT2  = 

PLOT[NORMDATA*INTENSITY[W,  LVAL]/NORMFIT,  {W,  DATA[[1]][[1]], 
DATA[[LEN]][[1]]}]; 

SHOW[PLTl,  PLT2,  PLOTLABEL  ->  "COMPARISON  OF  DATA  WITH  BEST  PEAK  FIT"]; 

); 


(*  FWHM  FINDS  FULL  WIDTH  AT  HALF  MAX  FOR  INTERPOLATED  DATA  *) 
FWHM[DATA_]  :=  ( 

CHEKPTS  =  100; 

LEN  =  LENGTH[DATA]; 

START  =  DATA[[1]][[1]]; 

ENDD  =  DATA[[LEN]][[1]]; 

PEAK  =  PEAKFIND[DATA]; 

CURVE  [X]  :=  INTERPOLATION  [DATA]  [X]; 

FIRST  =  CATCH[FOR[FIRST  =  START; 
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JJ  =  0,  JJ  <  CHEKPTS,  JJ++, 


LAMFRONT  =  START  +  JJ*(ENDD  -  START)/CHEKPTS; 

IF  [CURVE  [LAMFRONT]  >  PEAK[[2]]/2,  THROW  [LAMFRONT  ]  ] ; 

]]; 

LAST  =  CATCH  [FOR[L  AST  =  ENDD; 

KK  =  0,  KK  <  CHEKPTS,  KK++, 

LAMBACK  =  ENDD  -  KK*(ENDD  -  START)/CHEKPTS; 

IF  [CURVE  [LAMB  ACK]  >  PEAK[[2]]/2,  THROW  [LAMBACK]]; 

]]; 

LAST  -  FIRST 

); 

FWHMTHEORY[DATA_,  L  J  :=  ( 

CHEKPTS  =  100; 

LEN  =  LENGTH[DATA]; 

START  =  DATA[[1]][[1]]; 

ENDD  =  DATA[[LEN]][[1]]; 

THEORY  CURVE  [X]  :=  INTENSITY[X,  L]; 

PEAK  =  1/(FINDMINIMUM[1/THE0RYCURVE[XX],  {XX,  {START,  ENDD}}][[1]]); 
FIRSTLAM  =  CATCH  [FOR[FIRST  =  START; 

JJ  =  0,  JJ  <  CHEKPTS,  JJ++, 

LAMFRONT  =  START  +  JJ*(ENDD  -  START)/CHEKPTS; 

IF  [THEORYCURVE  [LAMFRONT]  >  PEAK/2,  THROW  [LAMFRONT]]; 

]]; 

LASTLAM  =  CATCH[FOR[LAST  =  ENDD; 

KK  =  0,  KK  <  CHEKPTS,  KK++, 

LAMBACK  =  ENDD  -  KK*(ENDD  -  START)/CHEKPTS; 
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IF  [THEORY  CURVE  [LAMB  ACK]  >  PEAK/2,  THROW  [LAMBACK]]; 

]]; 

LASTLAM  -  FIRSTLAM 

); 

(*  FINDS  BEST  L  BASED  ON  FITTING  WIDTHS  OF  THEORY  AND  DATA  *) 
BESTLFWHM[DATA_]  :=  LL  /.  FINDROOT[FWHMTHEORY[DATA,  LL]  ==  FWHM[DATA],  {LL, 
{1,100}}]; 

JLMATKHFWHM[L_,  DATA  ]  :=  ( 

BSTL  =  BESTLFWHM[DATA]; 

NORMDATA  =  DATASUM[DATA]; 

LEN  =  LENGTH[DATA]; 

NORMFIT  =  NINTEGRATE[INTENSITY[WW,  BSTL],  {WW,  DATA[[1]][[1]],  DATA[[ 
LEN]][[1]]}]; 

SUM [( INTEN SIT Y [D ATA [ [N] ] [ [  1  ] ] ,  L]/NORMFIT  -  DATA[[N]][[2]]/NORMDATA)A2,  {N,  1, 
LEN}] 


(*  BESTLSQPK  FINDS  BEST  L  BASED  ON  RENORM  MIN  SQUARE  DEVIATION  *) 
BESTLSQFWHM[DATA_]  :=  ( 

BSTL  =  BESTLFWHM[DATA]; 

NORMDATA  =  DATASUM[DATA]; 

LEN  =  LENGTH[DATA]; 

NORMFIT  =  NINTEGRATE[INTENSITY[WW,  BSTL],  {WW,  DATA[[1]][[1]],  DATA[[ 
LEN]][[1]]}]; 

{BSTL,  LL  /.  FINDMINIMUM[SUM[(INTENSITY[DATA[[N]][[1]],  LL]/NORMFIT  - 
DATA[[N]][[2]]/NORMDATA)A2,  {N,  1,LEN}],  {LL,  {1,  10}}][[2]] 
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); 

BESTLSQIT[DATA_,  TIMES  ]  :=  ( 

BSTL  =  BESTLFWHM[DATA]; 

NORMDATA  =  DATASUM[DATA]; 

LEN  =  LENGTH[DATA]; 

TRIALL  =  BSTL; 

FOR[JJ  =  1,  JJ  <  TIMES,  JJ++, 

NORMFIT  =  NINTEGRATE[INTENSITY[WW,  TRIALL],  {WW,  DATA[[1]][[1]], 
DATA[[LEN]]  [[  1  ]  ] }  ] ; 

TRIALL  =  LL  /.  FINDMINIMUM[SUM[(INTENSITY[DATA[[N]][[1]],  LL]/NORMFIT  - 
DATA[[N]][[2]]/NORMDATA)A2,  {N,  1,  LEN}],  (LL,  {1,  10}}][[2]]; 

]; 

TRIALL 

); 

COMPARESQIT[DATA_,  ITERATIONS  ]  :=  ( 

PLT1  =  LISTPLOT[DATA]; 

LVL  =  BESTLSQIT[DATA,  ITERATIONS]; 

PRINT["BEST  FITTED  L=",  LVL,  "NANOMETERS"]; 

NORMDATA  =  DATASUM[DATA]; 

LEN  =  LENGTH[DATA]; 

NORMFIT  =  NINTEGRATE[INTENSITY[W,  LVL],  {W,  DATA[[1]][[1]],  DATA[[LEN]][[1]]}]; 
PLT2  = 

PLOT[NORMDATA*INTEN SITY[W,  LVL]/NORMFIT,  {W,  DATA[[1]][[1]], 

DATA[[LEN]]  [[  1  ]  ] }  ] ; 

SHOW[PLTl,  PLT2,  PLOTLABEL  ->  "COMPARISON  OF  DATA  WITH  BEST  FIT"]; 

LVL 
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(*  FITNCOMPARE  GIVES  BEST  L  =  NANOCRYSTALLITE  SIZE  FOR 
ARRAYNM  DATAPOINTS  FROM  MINIMUM  OMEGA  WMIN  TO  MAXIMUM  OMEGA  WMAX, 
GIVING  A  BACKGROUND  SIZE  ABOVE  THE  MINIMUM  DATAVALUE  OF  BACKGRND, 
SMOOTHING  THE  DATA  BY  AVERAGING  OVER  PTS2AV  POINTS,  MINIMIZING  THE 
DEVIATIONS  OVER  SKIMPTS  FINAL  NUMBER  OF  POINTS,  AND  INTERATING 
INTERATIONS  TIMES  *) 

F ITN COMPARE [ARRA YNM_,  WMIN_,  WMAX_,  BACKGRND_,  PTS2AV_,  SKIMPTS_, 
ITERATIONS  ]  := 

COMPARESQIT[SKIMDATA[TREATDATA[SMOOTHWAVE[ARRAYNM,  PTS2AV,  WMIN, 
WMAX],  BACKGRND],  SKIMPTS],  ITERATIONS]; 

TKHDATPTSO  =  READLIST['’D:\SICT7A.TXTM,  {REAL,  REAL}]; 

TKHDATPTS7  =  READLIST["D:\SICT7A.TXT",  {REAL,  REAL}]; 

BACKGRND  =  10; 

TREATDATO  =  TREATDATA[SMOOTHWAVE[TKHDATPTSO,  3,  900,  1000],  BACKGRND]; 
TREATDAT7  =  TREATDATA[SMOOTHWAVE[TKHDATPTS7,  3,  900,  1000],  BACKGRND]; 

POINTS  =15;  (*  NUMBER  OF  SKIMMED  POINTS  TO  FIT  *) 

SKIMDAT0  =  SKIMDATA[TREATDAT0,  POINTS]; 

SKIMDAT7  =  SKIMDATA[TREATDAT7,  POINTS]; 

BESTL[DATA_,  TIMES_,  SKIMPOINTS  J  :=  ( 

BSTLL  =  BESTLFWHM[DATA]; 


SKIMPTS  =  SKIMDATA[DATA,  SKIMPOINTS]; 


NORMDATA  =  DATASUM[SKIMPTS]; 

LEN  =  LENGTH[SKIMPTS]; 

TRIALL  =  BSTLL; 

FOR[JJ  =  1,  JJ  <  TIMES,  JJ++, 

NORMFIT  =  NINTEGRATE[INTENSITY[WW,  TRIALL],  {WW,  SKIMPTS[[1]][[1]], 
SKIMPTS[[LEN]][[1]]}]; 

TRIALL  =  LL  /.  FINDMINIMUM  [SUM[(INTENSITY[SKIMPTS[[N]][[1]],  LL]/NORMFIT  - 
SKIMPTS[[N]][[2]]/NORMDATA)A2,  {N,  1,  LEN}],  (LL,  {1,  10}}][[2]]; 

]; 

TRIALL 


COMPAREPLOT[DATA_,  ITERATIONS  ,  SKIMPOINTS  J  :=  ( 

PLTO  =  LISTPLOT[TREATDATO,  PLOTJOINED  ->  TRUE, 

PLOTSTYLE  ->  RGBCOLOR[l,  0,  0],  PLOTLABEL  ->  "POINTSO"]; 

PLT1  =  LISTPLOT[DATA,  PLOTLABEL  ->  "TREATED  DATAPOINTS"]; 

LVL  =  BESTL[DATA,  ITERATIONS,  SKIMPOINTS]; 

NORMDATA  =  DATASUM[DATA]; 

LEN  =  LENGTH[DATA]; 

NORMFIT  =  NINTEGRATE[INTENSITY[W,  LVL],  {W,  DATA[[1]][[1]], 

DATA[[LEN]]  [[  1  ]] }  ] ; 

PLT2  =  PLOT[NORMDATA*INTENSITY[W,  LVL]/NORMFIT,  {W,  DATA[[1]][[1]], 
DATA[[LEN]][[1]]} 

PLOTLABEL  ->  "THEORETICAL  BEST  FIT  OF  TREATED  DATAPOINTS"]; 
PRINT["BEST  FITTED  L=",  LVL,  "  NANOMETERS"]; 
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SHOW[PLTO,  PLT1,  PLT2,  PLOTLABEL  ->  "COMPARISON  OF  DATA  WITH  BEST  FIT"]; 


LVL 

); 


REN ORMINTEN SITY[DATA_,  W_,  LLJ  :=  ( 

NORMDATA  =  DATASUM[DATA]; 

LEN  =  LENGTH[DATA]; 

NORMFIT  =  NINTEGRATE [INTENSITY [WWW,  LL],  {WWW,  DATA[[1]][[1]], 
DATA[[LEN]]  [[  1  ]  ] }  ] ; 

NORMDATA*INTENSITY[W,  LL]/NORMFIT 

); 


DATANTHEORY[DATA_,  LL  J  :=  ( 

LEN  =  LENGTH[DATA]; 

TABLE[{DATA[[NNN]][[1]],DATA[[NNN]][[2]], 

RENORMINTENSITY[DATA,  DATA[[NNN]][[1]],  LL]},  {NNN,  1,  LEN}] 

); 

PLTO  =  LISTPLOT[TREATDATO,  PLOTJOINED  ->  TRUE, 

PLOTSTYLE  ->  RGBCOLOR[l,  0,  0],  PLOTLABEL  ->  "POINTSO"]; 

PLT7  =  LISTPLOT[TREATDAT7,  PLOTJOINED  ->  TRUE, 

PLOTSTYLE  ->  RGBCOLOR[0,  0,  1],  PLOTLABEL  ->  "POINTS7"]; 
SHOW[PLTO,  PLT4,  PLT7,  PLOTLABEL  ->  "ALL  POINTS"]; 
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PointsO 


compareplot[treatdat7,  iterations  =  2,  skimpoints  =  15]; 


PointsO 
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Treated  Datapoints 
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Theoretical  Best  Fit  of  Treated  Datapoints 
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Comparison  of  data  with  best  fit 
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{{901.036,24.0333,38.6076},  {902.266,  48.5333,40.8618},  {903.496,  57.3333,43.3257},  {904.725, 
41.9333,46.0247},  {905.955,52.3,  48.9876},  {907.185,  110.767,  52.2468},  {908.415,  157.167, 
55.8396},  {909.644,  166.767,  59.8081},  {910.874,  110.433,  64.2003},  {912.104,  80.3667,  69.0705}, 
{913.333,  101.367,  74.4807},  {914.563,  145.,  80.5009},  {915.793,  173.233,  87.2109},  {917.023, 
203.633,94.7006},  {918.252,  247.867,  103.072},  {919.482,300.933,  112.441},  {920.712,308.567, 
122.937},  {921.942,  320.033,  134.706},  {923.171,329.933,  147.913},  {924.401,359.867,  162.741}, 
{925.631,389.067,  179.395},  {926.86,  407.7,  198.103},  {928.09,  407.633,219.116},  {929.32,387.433, 
242.712},  {930.55,  375.333,  269.192},  {931.779,  404.7,  298.883},  {933.009,  456.733,  332.135}, 
{934.239,479.3,369.32},  {935.469,481.2,410.824},  {936.698,478.9,457.041},  {937.928,  505.533, 
508.359},  {939.158,  541.333,  565.146},  {940.387,  564.5,  627.728},  {941.617,  575.767,  696.355}, 
{942.847,  597.8,  771.163},  {944.077,  633.233,  852.118},  {945.306,  699.433,  938.951},  {946.536, 
739.567,  1031.07},  {947.766,  767.833,  1127.45},  {948.996,  797.667,  1226.49},  {950.225,  834.233, 
1325.89},  {951.455,  909.233,  1422.44},  {952.685,  1011.7,  1511.86},  {953.914,  1162.43,  1588.63}, 
{955.144,  1247.77,  1645.99},  {956.374,  1276.57,  1676.09},  {957.604,  1275.6,  1670.69},  {958.833, 
1305.63,  1622.61},  {960.063,  1384.8,  1528.31},  {961.293,  1455.7,  1391.17},  {962.523,  1515.47, 
1223.69},  {963.752,  1507.97,  1045.38},  {964.982,  1433.07,  875.731},  {966.212,  1343.77,  727.346}, 
{967.441,  1215.1,604.359},  {968.671,  1052.57,505.248},  {969.901,825.967,426.189},  {971.131, 
585.5,  363.089},  {972.36,385.933,312.405},  {973.59,  235.033,271.324},  {974.82,  148.467,  237.689}, 
{976.05,97.0667,209.867},  {977.279,  57.9667,  186.627},  {978.509,  25.0333,  167.033},  {979.739,  10., 
150.371},  {980.969,21.3,  136.089},  {982.198,46.6667,  123.757},  {983.428,58.9333,  113.036}, 
{984.658,  63.5333,  103.658},  {985.887,  36.7667,  95.4082},  {987.117,  13.6,  88.112},  {988.347,  16.5, 
81.6279},  {989.577,  36.5,  75.8392},  {990.806,  45.1,  70.6498},  {992.036,  24.8667,  65.9794},  {993.266, 
42.9,  61.7607},  {994.496,  88.2667,  57.937},  {995.725,  131.833,  54.4605},  {996.955,  136.333,  51.29}, 
{998.185,  105.033,  48.3905},  {999.414,  74.1333,  45.732},  {1000.64,  62.7,  43.2881},  {1001.87,  82.4333, 
41.0364}} 
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